AD-A206  946  I 

I  *  * 

"llWBIp  I  I  IB  ■  WIMMI 

|  ?».  SECURITY  CLASSlFlCATl<3j|j^UT«ayT.jl  ^  ^  JQg 
L  Ot  CLASS: T I  C  A  T  ION/OC)‘;^^^^^iN(l  SCHCOULk 

1 _ _ |J 

[*.  P£ R  forming  organization  report  n^M^risi 


IggflT  DOCUMENTATION  PAGE 

p  p  fpST  |  to.  RESTRICTIVE  M6 


to.  RESTRICTIVE  MARKINGS  .;  if  !  ■;  ’  '  :  j  " 

3.  0'Stribut ion/a vai lab  i  lit  y  ok  report 

Approved  for  public  release; 
distribution  is  unlimited 

5.  monitoring  organization  report  NUNtssmsT 

RFOSR  TR-  8  9-  0  406 


3^4.  NAME  0*  PtAfO  MiNG  ORGANIZATION  9b.  OrPiCc  SYMBOL  ?«.  NAM£  OP  MONITORING  ORGANIZATION 
1  UfappUcaam 

i  Princeton  University  AJOSR/NA 


fits*.  AOOPaSS  f Cl ry .  St-’t  ana  ZIP  C#UI 

s 

I  I  Dept,  of  Meet.  &  Aero  Eng 

Princeton,  NJ  08544 

la  NAM*  OK  FUNOiNG/SPONSORiNG  to.  OFFlC 

ORGANIZATION  drip? 

AF0SR/NA 

ftc.  AQQRE3C  (City.  Slati  ana  ZH  CoOai  Building  410, 
Bolling  AFB ,  DC  20332-6448 


it.  TITl€  tfnctuC:  S*cu  t y  Clarification) 

_ ‘  Charged  Slurry  Droplet  Research 

13.  PERSONAL  AL)TnOR(S| 

A . J ■  Kelly _ 

13a  TVFF  OK  RSPORT  ~  1 3B.  Tl  ME  CO  VE  R€  D 

yinal  ; _ jFROMl/]0/86  TO  30/9738 

16.  SUPPLtMkN  CARY  NOVATION 


?S.  AOOTSSS  fCily.  Hilt  ana  ZIP  Ctuia)  „  .  ,  , 

j  Building  410, 

Bolling  AFB,  DC  20332-6448 


t*>.  office  Symbol  9.  procurement  instrument  identification  number 
nr  ippucatpK  1 

K)lV  \  AFOSR  86-00 j  3 


10.  SOURCE  of  funding  not 

program  froze  c 

E  Li  ME  NT  NO.  NO. 

-  t 

61102F  2308 


16.  O.STE  of  REFOTST  lYr..  Up..  Day)  1 B,  PAC,ij  COUN 

] 989  Feb  20  118 


PROJECT  . 

task 

HO. 

no. 

2308 

a2 

WORK  U* 
NO. 


j'L  \ 

CCEATl  CODES 

1  Field 

GROUP 

SUB  OR 

L  21  / 

01 

1  21'.  / 

_ 02 

IB  SUBJSqTOTSRMJ  IConhnitt  cn  rtvjtrut  U  niramary  t M  itntify  |y  Woe P  n>nM 

Charged  Droplets  •>  Charged  Sprays 

Rayleigh  Bursting)  Droplet  Levitations 

_ Droplet  Charging  Behavior  \  Quadrupole  Mass  Spectromet 


!  Rayleigh  Bursting,  wherein  critically  charged  droplets  explosively  expel  a  number  of  m 
cron  size  "sibling"  droplets,  enhances  atomization  and  combustion  of  all  liuqid  fuels. 
Droplet  :  -face  charge  is  retainer  during  evaporation,  permitting  multiple  Rayleigh  Burst 
to  occur.  Moreover,  the  charge  is  available  for  the  deagglomeration  of  residual  pvrticul 
floes  from  slurry  droplet  evaporation.  To  fill  gaps  in  our  knowledge  of  these  processes, 
an  experimental  program  involving  the  use  of  a  charged  droplet  levitator  and  a  quadrupole 
mass  spectrometer,  high  speed  electrometer  (QMS/HSE)  has  been  undertaken  to  observe  the 
disruption  and  to  measure  quantitatively  the  debris.  A  charged  droplet  levi  ator  based  c 
a  new  video  "frame  grabber"  technology  to  image  transient  events,  is  described.  Sibling 
a  droplet  size  is  ten  microns  or  less  and  is  close  to,  if  r.ot  coincident  with,  the  predict 
|  phase  transition  in  droplet  charging  level.  The  research  effort  has  focused  on  the  expL 
I  ation  of  this  transition  and  its  implications.  -x0ctoil  droplet  charge  and  ma  s  data,  spa' 
|  ning  a  thirty  fold  range  in  the  charge  to  mass  ratio  interval  centered  on  the  predicted 
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phase  transition,  verifies  the  existence  of  discontinuous  behavior.  In  particular,  the 
limiting  droplet  surface  electric  field  (the  emission  limit)  is  found  to  differ  by  a  factor  of 
5/3  *  8%  across  the  transition,  with  the  charging  distribution  for  droplets  smaller  than  the 
transition  having  the  higher  limit  and  lower  charge  level.  Droplets  larger  than  the 
transition  have  a  relative  dielectric  constant  2.1  ±  0.1  times  that  of  the  smaller  droplet 
population  and  a  minimum  emission  limit  of  4.2  GV/m.  The  data  are  consistent  with 
condensation  of  the  surface  electrons  from  a  plasma  state  to  a  rigid  electron  "crystal'’  as 
the  droplet  size  is  increased  through  the  transition.  Electron  crystallization  is  concluded  to 
dominate  both  the  Rayleigh  Bursting  process  and  electrostatic  spraying. 
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Introduction  — 

The  behavior  of  charged  and  uncharged  sprays  are  markedly  different.  The  unique 
characteristics  of  charged  sprays,  which  provide  options  not  available  from  conventional 
atomization,  recommend  their  consideration  for  combustion  systems  involving  difficult  to 
atomize  and  burn  fuels,  such  as  the  high  solid  content  slurries  of  current  interest. 

All  of  the  distinctive  characteristics  of  charged  sprays  are  directly  traceable  to  the 
presence  of  free  charge  on  the  individual  droplet  surface.  The  fact  that  each  droplet  is 
unipolarly  charged  means  that  droplet  dispersal  will  automatically  occur,  even  under  the 
most  demanding  fluid  mechanical  conditions.  As  a  consequence,  these  sprays  are 
inherently  homogeneous  and  less  susceptible  to  the  formation  of  fuel  rich  striations. 
Moreover,  droplet  agglomeration  will  be  suppressed  or  eliminated.  Beyond  these  important 
attributes  the  presence  of  free  charge  on  the  droplet  surface  yields  other  advantages  that 
cannot  be  offered  by  uncharged  sprays.  Two  of  these  advantages  involve  the  unique 
manner  in  which  a  charged  droplet  undergoes  evaporation. 

It  is  now  well  established  [1—5]  that  droplet  surface  charge  level  will  be  uninfluenced 
by  evaporation.  In  fact,  all  evidence  shows  that  the  charge  is  tenaciously  bound  to  the 
surface  and  will  not  be  released  no  matter  how  vigorous  the  evaporation.  This  is  true  for 
all  conditions  below  the  critical  charging  level,  the  so-called  Rayleigh  Limit  [6]  of  droplet 
charging.  As  evaporation  proceeds,  and  surface  charge  is  conserved  on  a  diminished 
surface,  the  stabilizing  surface  tension  force  is  inexorably  overtaken  by  the  monotonically 
increasing,  outwardly  directed  and  repulsive  electrostatic  force.  Although  the  charge 
induced  force  makes  the  droplet  more  susceptible  to  aerodynamic  shattering  than  it 
ordinarily  would  be  [7],  even  in  the  absence  of  motion,  a  critical  condition  is  reached 
wherein  the  electrostatic  and  surface  tension  forces  come  into  balance  and  the  droplet 
becomes  unstable. 

What  little  is  known  concerning  the  subsequent  behavior  of  the  droplet  is  consistent 
with  the  explosive  expulsion  of  several  highly  charged  "sibling"  droplets  having  a  diameter 
of  the  order  of  ten  microns.  As  a  result  of  this  "Rayleigh  Bursting"  process  the  "parent" 
droplet  behaves  like  an  insitu  atomizer  that  is  capable  of  repeatedly  undergoing  cycles  of 
evaporation  and  bursting  and  thereby  seeding  the  spray  volume  with  vigorously  dispersing 


siblings.  It  should  be  noted  that  the  siblings  also  undergo  this  forced  fragmentation  further 
adding  to  the  mixing  and  dispersal  process  [lj. 

As  charged  slurry  fuel  droplets  participate  in  this  episodic  bursting,  a  point  is 
reached  where  the  free  surface  charge  will  "plate  out"  on  the  residue  of  evaporatively  dried 
particulates  and  provide  a  mechanism  for  deagglomeration.  To  the  extent  that  neutralizing 
ions  from  the  combustion  zone  fail  to  penetrate  the  outwardly  streaming  and  insulating 
vapor  doud  surrounding  the  evaporating  and  shattering  droplet,  the  free  charge  is 
potentially  capable  of  completely  dispersing  the  slurry  droplet  and  making  it  accessible  for 
combustion.  It  is  this  promise  of  enhanced  atomization  and  combustion  that  has  provided 
the  stimulus  for  the  three  year  research  program  summarized  in  this  report. 

The  Research  Effort  - 

Rayleigh  Bursting  of  charged  droplets  and  the  electrostatic  deagglomeration  of 
micron  sized  particulate  flocculi  are  fundamental  processes  that  are  inadequately 
documented,  poorly  understood  and  are  commonly  accepted  to  involve  the  rapid  dispersal 
of  small  (1  to  10/zm)  charged  entities.  This  lack  of  fundamental  knowledge  can  be  placed  in 
perspective  by  the  realization  that  while  there  is  little  experimental  or  theoretical 
information  concerning  the  episodic  Rayleigh  Bursting  phenomena,  and  then  only  for 
theologically  simple  Newtonian  [8— 11]  fluids  (no  work  has  been  undertaken  on  high  solid 
contents  liquids),  nothing  is  known  concerning  the  deagglomeration  process.  No 
quantitative  data  exist  concerning  the  most  elemental  features  of  either  of  these  processes. 
In  fact,  even  for  the  more  extensively  studied  Rayleigh  Bursting  process,  the  temporal 
development  of  the  disruption,  the  number  of  particles  generated  and  their  mass,  charge, 
velocity,  spatial  distribution  and  timing  have  yet  to  be  observed.  This  research  effort  was 
specifically  directed  toward  filling  in  gaps  in  our  knowledge  of  these  phenomena  and 
thereby  to  provide  insights  into  the  larger  concern  of  whether  charged  sprays  and  the  use  of 
electrostatic  atomization  can  benefit  slurry  fuel  combustion. 

Overview  — 

In  the  absence  of  available  diagnostic  tools  that  are  capable  of  making  the  types  of 
measurements  indicated,  the  research  program  has  been  concerned  with  the  development 
and  validation  of  the  necessary  instrumentation.  This  effort  has  involved  the  development 
of  two  new  research  instruments:  a  levitated  droplet  apparatus  for  imaging  the  disruption 
of  isolated  droplets/agglomerates  and  a  high  mass  (1014AMU)  quadrupole  mass 
spectrometer/high  speed  electrometer  (QMS/HSE)  to  measure  both  the  charge  and  the 
mass  of  the  post  disruption  debris.  These  instruments  are  designed  specifically  to  probe 
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the  behavior  of  Theologically  complex  charged  slurry  fluids,  provide  essential  insight  into 
electrostatic  particulate  deagglomeration  and  to  bolster  our  meager  knowledge  of  charged 
fluid  droplet  dynamics. 

Both  the  charged  droplet  levitator  and  the  quadrupole  mass  spectrometer/high 
speed  electrometer  have  been  designed  to  be  used  individually  or  as  a  combined 
experiment.  The  charged  droplet  levitator,  a  modern  version  of  the  Millikan  oil  droplet 
experiment  is  capable  of  suspending  single  charged  droplets  to  make  them  accessible  for 
prolonged  observation.  By  using  active  feedback  control  in  conjunction  with  a  precision 
optical  imaging  system,  the  evaporating  droplet  can  be  held  stationary  in  the  optical  field 
of  view  by  appropriate  manipulation  of  the  levitating  electrode  voltages.  Although  the 
optical  system  is  used  for  the  detailed  control  of  the  suspended  droplet,  its  primary 
function  is  to  image  the  breakup  (bursting)  and  deagglomeration  processes  as  they  occur. 
Because  these  processes  are  estimated  to  occur  on  microsecond  to  millisecond  time  scales, 
which  are  orders  of  magnitude  shorter  than  those  required  for  levitation,  it  has  been 
necessary  to  design  the  system  to  have  comparable  response  characteristics.  Consequently, 
considerable  effort  has  been  devoted  to  this  aspect  of  the  imaging  system  design.  For 
convenience,  the  initial  levitator  experiment  was  designed  to  operate  at  atmospheric 
conditions.  While  this  mode  of  operation  will  permit  imaging  of  the  disruption  process  and 
will  permit  additional  data  concerning  the  Rayleigh  Limit  to  be  collected,  it  is  only  when 
the  levitator  is  used  in  conjunction  with  the  QMS/HSE  instrument  operating  at  high 
vacuum  that  quantitative  data  concerning  the  Rayleigh  Bursting  and  deagglomeration 
debris  will  be  obtained. 

As  the  primary  instrument  for  the  measurement  of  particle  charge  and  mass,  the 
development  of  the  QMS/HSE  has  been  the  major  focus  of  this  effort.  One  of  the  most 
challenging  aspects  of  this  work  has  been  the  calibration  of  this  instrument.  Conventional 
quadrupole  mass  spectrometer  devices  designed  for  use  in  the  1  to  1000 AMU  can  be 
calibrated  using  well  characterized  chemical  species.  Despite  the  comforting  agreement 
between  the  theoretical  and  observed  behavior  of  these  instruments  there  is  no  guarantee 
that  the  device  we  fabricated  would  also  exhibit  predictable  behavior  when  dealing  with 
droplets  and  particulates  in  the  lO^AMU  range.  The  disparity  in  mass  and  in  charge  to 
mass  ratio  is  simply  too  large  to  warrant  extrapolation.  The  substantial  effort  required  to 
validate  QMS/HSE  operation  at  the  extreme  mass  range  associated  with  10  micron  sized 
particulates  has  been  well  rewarded  by  the  fundamental  insight  this  has  provided  into  the 
behavior  of  charged  droplets.  This  effort  is  summarized  after  first  a  discussion  of  the 
levitator  work  is  presented. 
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Charged  Droplet  Levitator  — 

The  charged  droplet  levitator  apparatus  depicted  in  the  schematic  of  Figure  1  has 
been  specifically  designed  to  provide  direct  quantitative  information  concerning  the 
behavior  preceding  disruption  and  a  quantitative  measure  of  the  number  of  siblings 
produced.  This  apparatus,  an  extensively  refined  variant  of  the  design  first  described  by 
Rhim,  et  al.  [11,12],  is  described  fully  in  the  research  thesis  in  Appendix  A.  This  document 
should  be  consulted  for  details  of  construction  and  operation. 

As  shown  in  the  schematic  a  Compaq  386  computer  is  used  to  control  all  functions 
of  the  droplet  levitation  and  manipulation  processes.  During  operation,  the  Compaq 
computer  issues  a  launch  command  to  form  and  to  elevate  electrostatically  a  charged 
droplet  form  the  capillary  source  that  is  centrally  located  within  the  lower  electrode.  Tap 
water,  selected  for  convenience  for  all  testing,  is  sufficiently  conductive  that  the  application 
of  an  appropriate  voltage  on  the  upper  electrode  will  induce  sufficient  free  charge  to  Dow  to 
the  fiuid  at  the  tip  of  the  grounded  capillary  to  produce  a  charged  droplet. 

Rhim  has  demonstrated  that  the  most  important  prerequisite  for  successful  droplet 
formation  and  subsequent  projection  with  the  interelectrode,  observation  area  is  the  form 
of  the  voltage  (electric  field)  profile.  The  most  reliable  profile  involves  a  rapid  rise  to  an 
elevated  voltage  which  is  maintained  for  sufficient  time  for  the  droplet  to  form  and  detach 
from  the  capillary  tip.  Once  free,  the  droplet  is  then  subjected  to  a  diminished  electric 
field  of  sufficient  magnitude  to  continue  the  upward  trajectory  but  of  insufficient  strength 
to  further  accelerate  the  droplet  and  have  it  collide  with  the  upper  electrode.  The 
computer  then  issues  a  command  for  a  sharply  defined  retarding  voltage  pulse  to  be  applied 
to  the  lower  electrode.  This  serves  to  decelerate  the  droplet  and  to  bring  it  within  the  field 
of  view  of  the  CCD  camera  optics  that  will  then  serve  to  provide  position  and  velocity 
information  for  further  control  of  the  droplet  position.  As  discussed  in  the  appendix,  the 
high  speed,  high  voltage  electrode  voltage  supplies  are  the  most  essential  element  of  this 
first  phase  of  the  levitation  process. 

Once  droplet  stabilization  has  been  achieved,  the  droplet  is  now  maintained  within 
the  field  of  view  of  the  CCD  optics.  During  this  "observation  phase"  the  full  capability  of 
a  real  time  "frame  grabber"  system  is  used.  This  function  is  performed  by  the  Pyramid 
Tracker  computer,  &  proprietory  development  of  the  GE/SRI  David  Sarnoff  Laboratories 
[13,14].  The  Pyramid  Tracker  is  specifically  designed  to  automatically  analyze  visual  data 
form  the  CCD  camera  at  any  of  six  resolution  levels  ranging  from  level  1,  the  crudest  and 
most  easily  and  hence  most  rapidly  processed  image  of  7«8  pixels,  to  the  maximum 
resolution  of  the  full  240*250  pixel  image  of  level  6.  By  automatically  adjusting  the  image 
resolution,  the  Pyramid  Tracker  can  rapidly  process  the  visual  data  to  generate  droplet 
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Figure  1:  System  Block  Diagram 
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position  control  signals.  This  is  typically  accomplished  using  low  resolution  level  1 
imaging.  Once  stabilized,  the  droplet /agglomerate  is  then  imaged  at  level  6  in  anticipation 
of  the  disruption  process.  In  this  way  it  is  possible  to  use  the  same  system  to  position  the 
particle  and  to  capture  the  disruption  process.  All  of  the  necessary  software  has  been 
written  and  all  of  the  D/A  interface  electronics  for  levitation  are  now  operational. 

In  a  yet  to  be  published  paper,  Elghazaly  and  Castle  [15]  describe  the  results  of  a 
droplet  bursting  study  in  which  the  charged  droplet  is  immobilized  or  "tethered"  at  the  end 
of  a  string.  There  are  several  reasons  that  recommended  the  use  of  this  technique  for  the 
initial  validation  testing  of  the  Pyramid  Tracker  system.  First,  the  use  of  an  immobilized 
droplet  permits  the  Pyramid  Tracker  system  to  be  tested  independent  of  the  dynamic 
levitator  unit.  Secondly,  while  the  limited  data  on  droplet  description  provided  by  this 
method  do  agree  with  Elghazaly ’s  model,  the  data  base  is  too  sparse  to  establish  the 
general  validity  of  the  model.  Finally,  if  the  Pyramid  Tracker  data  demonstrates  that  the 
technique  does  permit  a  tethered  droplet  to  Rayleigh  burst  in  the  same  manner  as  a  truly 
isolated  levitated  droplet,  it  would  be  far  simpler  to  use  with  the  QMS/HSE  than  the 
levitator,  where  the  test  droplet  has  to  be  remotely  manipulated  under  high  vacuum 
conditions. 

Accordingly  the  tethered  droplet  experiment  was  implemented  and  tests  were 
undertaken  with  the  Pyramid  Tracker.  In  these  tests  a  water  droplet  is  extruded  from  a 
300/rm  OD,  250/zm  ID  hypodermic  tube  sting,  which  served  not  only  as  a  conduit  for 
replenishment  of  the  test  fluid  but  also  as  the  charging  electrode.  In  all  instances  the 
charging  voltage  was  adjusted  to  be  within  the  range  defined  by  Elghazaly  [15]  as  being 
required  for  bursting  rather  than  spraying  or  dripping  to  occur.  As  a  result  of  these  tests,  a 
sample  of  which  is  included  as  Figure  2,  the  following  conclusions  have  been  drawn. 

1.  The  Pyramid  Tracker  computer/optics  system  is  capable  of  consistently 
sensing  droplet  breakup  and  can  track  the  siblings  as  they  are  energetically 
repelled  from  the  parent. 

2.  The  sampling  frequency  of  the  Pyramid  tracker  is  adequate  for  preliminary 
validation  testing  but  will  have  to  be  improved  by  a  factor  of  at  least  two  if 
quantitative  imaging  of  the  bursting  process  is  to  be  obtained. 

3.  The  observed  droplet  breakup  pattern  is  consistent  with  the  results  obtained 
by  Elgazhaly.  The  sampling  limitations  imposed  by  the  Pyramid  Tracker 
system  do  not  permit  a  more  detailed  validation  to  be  made  at  this  time. 
Nevertheless,  bursting  and  sibling  expulsion  do  occur  and  have  been 
observed. 
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OCCURRENCE  RATE 


This  graph  represents  data  obtained  for  50  bunts.  The  pulse 
used  in  these  bursts  had  a  magnitude  of  14kV ,  a  duration  of 
60ms  and  a  transient  time  of  12ms. 
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Figure  2:  NUMBER  OF  DROPLETS  PER  BURST 


DISTRIBUTION  OF  THE  NUMBER  OF  DROPLETS 
FORMED  IN  RZSPONCE  TO  A  PULSE 


The  inherent  simplicity  of  the  tethered  droplet  technique,  as  compared  to  the 
levitator  systems  led  to  the  decision  to  focus  on  the  use  of  this  method  rather  than  the 
levitator  with  the  QMS/HSE  instrument.  This  decision  permitted  us  to  focus  our 
attention  on  the  development  of  the  QMS/HSE  instrument. 

Quadruoole  Mass  Spectrometer/High  Speed  Electrometer  - 

Accordingly,  the  primary  focus  of  our  activity  has  been  the  development  of  a 
Quadrupole  Mass  Spectrometer/High  Speed  Electrometer  that  is  capable  of  independently 
measuring  the  charge  to  mass  ratio  and  the  charge  of  individual  droplets  and  particulates 
in  the  size  range  of  interest  (order  of  ten  microns).  Figure  3  is  a  schematic  of  the 
instrument  we  have  built  and  which  is  fully  described  in  the  thesis  of  Appendix  B.  While 
the  schematic  depicts  the  droplet  source  as  a  diode  sprayer,  the  device  was  designed  to  be 
coupled  with  the  droplet  levitator  and  with  the  tethered  droplet  rig  of  Elghazaly.  As 
described  in  the  appendix,  which  should  be  consulted  for  technical  detail,  the  device  has 
been  designed  to  filter  droplets  having  charge  to  mass  ratios  in  the  0.01  to  lOOC/kg  (10  to 
10  C/m  )  range.  However,  limitations  imposed  by  the  availability  of  appropriate  high 
voltage  transformers  effectively  restricted  the  charge  to  mass  capability  to  ~0.1  to  ~3C/kg, 
corresponding  to  the  droplet  size  range  of  mos*  interest.  Commercially  available  Amptek 
charge  detectors  and  amplifiers  provide  a  maximum  charge  sensitivity  of  280  electrons. 

Droplets  or  particulates  having  charge  to  mass  ratios  within  the  bandwidth  of  the 
QMS,  which  typically  is  set  at  3  to  5%  of  the  preselected  value  in  the  range  noted,  will  pass 
to  a  collector  plate  where  their  charge  is  transferred  for  measurement.  As  will  be  noted, 
particular  attention  had  to  be  devoted  to  the  charge  collection  process  to  assure  that  all 
droplet  charge  is  drained  to  the  amplifier  circuit  within  the  response  time  of  the  detector. 
Once  collected,  the  amplifier  generates  an  output  signal  that  is  proportional  to  the  charge 
collected.  It  is  this  signal  that  produces  the  droplet  population  vs  charge  level  spectrum  on 
the  multichannel  analyzer  (MCA)  shown  in  Figure  4.  These  data  and  the  operating  state 
of  the  droplet  source  of  the  quadrupole  mass  spectrometer  form  the  raw  data  suite  that  is 
used  as  a  basis  for  analysis  of  droplet  behavior. 

Calihiatioi}  - 

The  most  crucial  and  time  consuming  aspect  of  the  experiment  involved  calibration. 
Low  mass  (high  charge  to  mass)  instruments,  which  can  be  calibrated  with  well  defined 
atomic  and  molecular  species,  have  exhibited  excellent  agreement  of  predicted  and  actual 
performance.  Although  comforting,  the  droplet  charge  to  mass  range  is  too  far  removed 
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F i gure  4 :  Channel  Number 


from  typical  experience  to  permit  confident  extrapolation.  In  the  absence  of  a  clearly 
defined  calibration  source,  in  the  range  of  interest,  it  is  necessary  to  resort  to  the  use  of 
droplet  data  obtained  from  an  independent  measurement  technique  for  comparison.  This 
other  technique,  a  time  of  flight  charge,  charge  to  mass  ratio  measurement  method,  [18] 
also  uncalibrated,  has  the  advantage  of  minimum  data  reduction  complexity.  Since  the 
only  such  data  involve  the  use  of  an  octoil  (dioctal  phthalate)  spray,  an  octoil  sprayer  was 
developed  for  use  with  the  QMS/HSE  experiment. 

Only  if  the  results  of  both  experiments  are  in  agreement  can  there  be  confidence 
that  the  QMS/HSE  is  providing  quantitative  data.  Otherwise  there  is  no  means  by  which 
a  judgment  of  the  validity  of  either  technique  can  be  made.  The  comparison  of  data  sets 

obtained  by  both  techniques,  shown  in  Figure  5,  displays  the  level  of  agreement  that  exists. 

3  3 

The  initial  comparison  made  with  high  charge  to  mass  ratio  (>  10  C/m  )  is  in  excellent 

agreement  with  the  time  of  flight  data  of  Hendricks.  This  congruence  of  the  two 
experiments  is  sufficiently  precise  for  us  to  conclude  that  the  QMS/HSE  experiment  and 
the  time  of  flight  measurements  are  both  quantitative.  It  is  worth  noting  that  whereas  the 
Hendricks  experiment  involved  manual  data  collection  and  a  very  limited  number  of  spray 
droplets  the  computerized  QMS/HSE  experiment  permitted  two  orders  of  magnitude  larger 
data  base  to  H  developed. 

O  O 

Subsequent  tests  conducted  at  charge  density  levels  below  10  C/m  ,  the  range  of 
most  interest,  produced  data  that  increasingly  diverged  with  decreasing  charge  density 
level  from  that  published  by  Hendricks.  This  behavior  was  ultimately  traced  to  charge 
relaxation  at  the  detector  surface,  as  discussed  in  Appendix  B.  Briefly,  what  we  discovered 
is  that  the  relatively  high  resistivity  of  the  octoil  used  in  the  test  produces  charge 
relaxation  times  for  the  transfer  of  droplet  charge  to  the  metallic  charge  detector  plate  that 
are  comparable  to,  or  larger  than,  the  gated  charge  acceptance  time  scale  of  the  charge 
detector  circuit.  The  reason  this  effect  only  occurs  for  the  larger,  less  highly  charged 
droplets  is  simply  related  to  the  relaxation  time  scale  being  a  function  of  conductivity, 
which,  in  turn,  is  dominated  by  the  mobility,  charge  density  term.  Once  this  factor  was 
taken  into  account,  excellent  agreement  between  the  two  data  sets  is  obtained  over  the 
entire  charge  density  range,  as  displayed  in  Figure  5. 

Octoil  Sprav  Testing  - 

Having  established  that  the  QMS/HSE  instrument  is  capable  of  measuring  both  the 
charge  and  the  mass  of  all  droplets  (and  deagglomeration  debris)  having  charge  to  mass 
ratios  between  3C/kg  and  O.lC/kg  and  with  total  surface  electron  counts  larger  than  280, 
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the  research  was  directed  toward  the  use  of  this  instrument  to  analyze  the  octoil  spray. 
This  effort  was  motivated  by  several  factors  First,  this  study  would  provide  a  convenient 
means  to  become  familiar  with  the  apparatus  and  to  assess  its  capabilities.  Second,  little 
quantitative  droplet  spray  data  exists  and  the  results  of  the  investigation  would  represent  a 
valuable  contribution  to  the  sparse  data  base.  Third,  data  would  be  obtained  in  the 
droplet  size  and  charging  realm  in  which  the  Statistical  Equilibrium  Model  of  electrostatic 
spraying  predicts  the  existence  of  an  enigmatic  phase  transition,  about  which  nothing  is 
known.  And,  fourth  and  finally,  it  was  anticipated  that  the  availability  of  quantitative 
information  concerning  droplet  behavior  in  this  charge  and  mass  interval  would  provide 
insight  into  the  physics  of  surface  charging  and  would  make  it  possible  to  understand  why 
it  is  necessary  to  invoke  the  disquieting  assumption  of  limited  surface  mobility  in  order  to 
model  the  Rayleigh  Bursting  phenomena  [8].  It  is  this  question  above  all  that  motivated 
the  decision  to  devote  the  bulk  of  our  resources  to  this  study. 

Background  — 

In  addition  to  the  octoil  data  previously  mentioned,  comparable  droplet  mean 
charge  to  mass  ratio  vs  radius  data  exist  only  for  Woods’  metal  [19]  and  glycerin  [20]. 
Despite  the  sparseness  of  this  data  base,  it  is  sufficient  to  set  as  a  first  test  of  the  validity 
of  the  Statistical  Equilibrium  Theory  of  electrostatic  atomization  [21-23].  This  is 
illustrated  by  the  comparison  of  the  octoil  and  Woods’  metal  data  and  theory  presented  in 
Figure  6.  The  Statistical  Equilibrium  Theory  is  a  thermodynamic  model  insofar  as  it  is 
based  on  the  premise  that  the  long-range  electrostatic  force  field  in  charged  fluids  will 
guarantee  that  all  physical  states  of  the  spray  system  are  "well-mixed",  are  in  intimate 
communication  and  result  in  the  most  probable  end  state  —  the  maximum  entropy 
(equilibrium)  end  state  [21-23].  These  states  are  shown  to  be  describable  in  terms  of  the 
product  of  two  coupled  distributions;  the  Maxwell— Boltzmann  distribution  for  droplet 
number  and  Fermi— Dirac  for  charging.  Use  of  the  Lagrangian  multiplier  method,  to  define 
the  most  probable  state,  results  in  the  following  two  descriptive  equations  for  all  Rayleigh 
Limit  restrained  electrostatic  atomization: 
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where  a',  /?',  and  are  the  Lagrangian  multipliers  representing  the  conservation  law 

restraints  of  charge,  energy  and  system  mass,  respectively.  The  system  is  defined  in  terms 

of  droplet  number  N  and  the  charge  parameter  Z,  defined  as  the  fraction  of  Rayleigh  limit 

charging.  Within  the  context  of  this  model,  the  surface  tension  (7),  fluid  density  ( p ),  and 

surface  electron  field  emission  limit,  E  ,  are  the  only  fluid  properties  of  fundamental 

s 

importance  to  the  description  of  electrostatic  spraying.  The  permittivity  t  and  the 

0 

electron  charge  q  are  the  only  other  parameters  in  the  model. 

Comparison  of  the  charging  equation  with  data,  as  illustrated  in  Figure  6,  yields  the 
startling  result  that  the  ratio  of  the  charging  and  energy  Lagrangian  multipliers  (a'//?)  is  a 
constant  and  is  numerically  equal  to  ~  -66eV.  Physically,  this  means  that,  on  average,  a 
droplet  in  either  the  Woods’  metal  and  the  octoil  sprays,  and  by  implication  all  sprays,  are 
charged  to  levels  that  require  an  expenditure  of  ~66eV  to  bring  another  charge  to  their 
surface.  This  unexpected  result  must  reflect  an  underlying  fundamental  mechanism  insofar 
as  all  of  the  physical  properties  of  these  materials  differ  by  at  least  an  order  of  magnitude 
and  the  individual  values  of  o'  and  0'  differ  by  approximately  four  orders  of  magnitude. 
Additional  work  has  verified  that  —a'  10'  is  a  constant  for  a  wide  variety  of  fluids 
including  No.  2  Heating  Oil  and  Jet-A,  to  name  but  two. 

Equally  startling  are  the  ramifications  of  the  fact  that  a'  /0'  is  a  "constant  of 
motion"  for  spraying.  This  can  be  seen  in  Figure  7,  which  illustrates  the  theoretical  mean 
droplet  charging  behavior,  derived  from  the  charging  equation  when  -a’ 10'  =  66eV. 
Droplet  charging  is  predicted  to  undergo  a  first  order  phase  transition  that,  in  the  case  of 
hydrocarbon  oils,  occurs  at  a  mean  diameter  of  about  one  micron.  Moreover,  the  behavior 
of  sprays  having  droplet  Bizes  larger  than  the  transition  diameter  is  described  by  the 
energetically  favored  middle  solution  branch.  In  fact,  all  of  the  available  data  have  been 
obtained  in  this  regime;  none  is  available  for  droplet  sizes  smaller  than  the  transition. 

In  addition,  the  theory  predicts  that  along  this  central  solution  branch  (branch  B) 
the  mean  droplet  size  will  be  independent  of  all  fluid  properties  and  will  be  determined 
uniquely  by  the  fluid  charge  density  level.  This  pertains  to  all  sprays  on  the  large  droplet 
side  of  the  transition  providing  the  physical  properties  of  the  spray  fluid  are  not  extreme 
(eg.  viscosity  <  200cp)  [23].  The  independence  of  mean  droplet  size  on  fluid  property 
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FIG. 


implies  that  the  spray  process  is  controlled  by  strictly  electronic  phenomena.  Most 
intriguing  is  the  indication  that  the  size  of  sibling  droplets  produced  by  Rayleigh  Bursting 
are  in  the  vicinity  of  the  transition  size.  Collectively,  these  factors  were  sufficiently 
compelling  to  make  the  study  of  this  transition  region  the  focus  of  our  research  effort.  It 
can  be  argued  that  since  the  transition  establishes  the  scale  of  spraying,  it  also  dominates 
the  Rayleigh  Bursting  process,  which  is  a  more  restricted  form  of  the  spray  process. 

QMS/HSE  and  the  Phase  Transition  - 

The  phase  transition  for  Woods’  metal  and  octoil  has  the  form  noted  in  Figure  8 
when  plotted  in  terms  of  the  charge  to  mass  ratio,  droplet  charging  parameter  space  of  the 
QMS/HSE  experiment.  Although  the  experiment  will  ultimately  be  capable  of  operation 
within  a  full  four  order  of  magnitude  charge  to  mass  ratio  range  (0.01  to  lOOC/kg),  the 
available  range  of  ~0.1  to  ~3C/kg  fortuitously  brackets  the  predicted  octoil  phase 
transition.  Note  that  only  the  mean  droplet  charge  level  is  displayed  in  this  figure. 
Consequently,  since  there  is  a  distribution  of  droplet  charging  levels  about  the  mean, 
operation  at  fixed  charge  to  mass  ratio  will  produce  a  "charging  spectrum"  (droplet 
number  vs  charge  level)  that  will  consist  of  three  peaks,  one  each  for  the  A,  B  and  C 
branches.  If  present,  these  peaks  will  be  t'*adened  by  the  small  but  finite  bandwidth  of 
the  QMS/HSE,  which  is  typically  maintained  at  3  to  5%  of  the  mean  charge  to  mass  ratio, 
as  explained  in  Appendix  B. 

Appendix  B  should  also  be  consulted  for  details  of  the  manner  in  which  the  charge 
spectrum  data  are  reduced.  Briefly,  the  charging  and  droplet  number  distributions  have 
been  recast  in  terms  of  the  measurables,  charge  to  mass  ratio  and  charge,  and  have  been 
reformulated  to  permit  the  fundamental  spray  (a',  0  ,  S')  and  droplet  property 
parameters  (c,  EJ  to  be  evaluated.  Attempts  to  fit  the  data  with  a  single  distribution 

proved  fruitless.  Independent  of  the  charge  to  mass  ratio  of  the  data  used,  the  SAS 
Institute  non-linear  regression  program  [24]  that  has  been  used  for  all  data  reduction  was 
always  able  to  obtain  a  more  accurate  fit  to  the  data  when  using  two  distributions  as 
opposed  to  when  one  distribution  was  fit.  AH  attempts,  in  which  the  program  was  asked  to 
fit  the  theoretically  expected  three  distributions  (one  each  for  branch  A,  B  and  C),  were 
indistinguishable  from  the  two  distribution  result.  This  outcome  is  to  be  anticipated  since, 
in  the  charge  to  mass  ratio  being  scanned,  the  presence  of  branch  C  would  be  in  the  noise 
level  of  the  experiment. 

Figures  9  through  13  show  the  results  of  the  automatic  fit  of  dual  distributions  to 
the  QMS/HSE  data.  Quantitatively,  the  results  conform  to  expectation.  At  the  highest 
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Figure  8:  SURFACE  ELECTRON  COUNT 
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Figure  9:  Radius 
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Figure  10:  Radius  (Microns) 
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Figure  II:  Radius 
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Figure  12:  Radius  (Microns) 
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charge  to  mass  ratio  accessible  with  the  QMS/HSE  (Figure  9,  2.93C/kg),  the  dominant 
feature  is  the  branch  A,  and  little  if  any  contribution  is  to  be  made  by  branches  B  and  C 
(cf.  Figure  8).  This  operating  state  is  close  to  the  theoretically  predicted  small  droplet, 
high  charge  density  side  of  the  transition.  Subsequent  figures,  obtained  at  successively 
lower  charge  to  mass  ratios,  show  the  emergence  of  the  central  "B"  branch.  In  addition, 
the  data  are  also  consistent  with  theoretical  expectations  with  respect  to  the  drift  of  the 
distribution  peaks  for  the  A  and  B  branches  to  successively  larger  droplet  sizes  as  a 
function  of  increasing  charge  density  level. 

It  is  obviously  desirable  to  probe  the  high  charge  density  realm  upward  from  the 
~3C/kg  level  of  the  present  experiment  and  away  from  the  transition  to  demonstrate 
unequivocally  that  spray  behavior  continues  to  be  monomodal  as  demanded  by  theory. 
Equally  desirable  would  be  the  collection  of  data  at  lower  charge  density  levels  than  the 
0.118C/kg  of  Figure  13  and  away  from  the  transition  in  the  large  droplet  direction  to 
establish  how  the  dominance  of  the  central  B  branch  is  attained  and  spray  monomodality  is 
restored.  However,  pending  the  modification  of  the  experiment  and  the  collection  of  these 
data,  sufficient  data  have  now  been  obtained  to  permit  a  number  of  definitive  statements 
concerning  the  phase  transition  and  the  behavior  of  charge  on  the  droplet  surface  to  be 
made. 

Results  — 

_ C 

Since  the  QMS/HSE  experiment  is  conducted  under  hard  vacuum  (~10  Torr),  and 

since  evaporation  of  the  octoil  is  negligible,  the  permittivity  of  background  medium  was 

taken  to  be  that  of  vacuo.  By  fixing  e  =  £o,  the  distribution  function  coefficients  obtained 

from  the  SAS  Institute  non-linear  regression  fit  program  can  be  used  to  extract  values  for 

the  limiting  electric  field  strength  E  .  The  results  of  this  analysis  are  displayed  in  Figure 

0 

14  where  the  limiting,  emission  electric  field  E  is  shown  plotted  verse  charge  to  droplet 

s 

mass  ratio  as  determined  by  the  quadrupole  mass  spectrometer. 

Several  features  of  these  data  are  noteworthy.  First,  the  behavior  of  Eg  for  each  of 

the  branches  (Ega  &  Eg^)  exhibit  the  same  trends  with  charge  to  mass  ratio.  Secondly,  the 

measured  values  are  related  by  a  factor  of  5/3.  That  is,  for  all  of  the  data  available, 
covering  approximately  a  twenty  fold  variation  in  Q/M  (0.12  to  2.2),  the  ratio  of  the  B 
branch  to  the  A  branch  limiting  electric  field  is  Eg^/Ega  =  5/3  ±  8%.  Although  the  data 

are  limited,  the  limiting  electric  field  can  be  reasonably  taken  to  be  invariant  with  charge 
to  mass  ratio  for  charge  to  mass  ratios  greater  than  about  lC/kg.  More  data  are  required 
in  the  range  between  1  and  ~3C/kg,  where  the  B  branch  disappears  and  the  quality  of  the 
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available  test  results  is  questionable.  However,  based  on  the  three  data  runs  we  have 
obtained,  the  B  branch  droplets  have  a  limiting  field  strength  of  ~9GV/m  in  contrast  to 
~5GV/m  on  the  high  charge  density  A  branch. 

In  the  charge  density  realm  below  ~lC/kg,  the  A  branch  is  truncated,  as 
anticipated,  at  a  charge  density  level  of  about  O.lC/kg.  On  the  basis  of  the  available  six 
data  points,  the  data  neatly  describe  a  linear  relationship  of  E.  and  Q/M.  It  is  interesting 

to  note  that  the  minimum  value  for  the  A  branch  (at  ~0.lC/kg)  is  2.8GV/m,  which  is 
comparable  to  4.2GV/m  value  projected  for  the  B  branch  for  low  Q/M  values.  These  low 
values  of  emission  field  are  compatible  with  published  values  for  negative  (electron) 
emission  [25—27]. 

As  yet  there  is  no  explanation  for  these  trends,  nor  do  we  have  any  basis  for 
understanding  the  distinct  regularity  of  the  E^/E^  ratio.  A  provocative  feature  of  these 

data  is  the  relationship  that  exists  between  dielectric  permittivity  e  and  the  limiting 

i 

electric  field  E.  These  two  parameters  appear  as  a  product  having  the  form  rE  .  When 

0 

interpreted  in  terms  of  a  fixed  immutable  limiting  electric  field  and  a  variable  permittivity, 
all  of  the  available  data  (for  0.1  <  Q/M  <  2.2)  can  be  summarized  in  terms  of  the  effective 
dielectric  constant  for  the  large  droplet  **  branch  being  some  2.1  ±  0.2  times  that  of  the 
small  droplet  A  branch. 

It  is  this  curious  result  that  has  strongly  suggested  that  a  fundamental  shift  in  the 
configuration  of  surface  charge  is  occurring  at  the  transition.  That  is,  the  data  imply  that 
the  surface  charge  is  in  a  plasma  state  (low  permittivity)  on  the  small  droplet  side  of  the 
transition  and  is  in  a  condensed  state  (crystallized)  on  the  large  droplet  side  —  the  side 
defined  by  the  B  branch.  At  this  point  our  lack  of  sufficient  data  makes  this  an  inferential 
conclusion.  However,  support  for  this  conclusion  can  be  found  in  what  has  been  established 
concerning  the  behavior  of  droplet  systems  on  the  large  droplet  B  branch.  The  evidence  is 
as  follows: 

1.  The  first  principles  Analysis  of  Rayleigh  Bursting  presented  in  the 
Roth/Kelly  paper  (8)  requires  that  the  surface  charge  mobility  be  small  or 
zero.  No  other  assumption  permits  a  model  to  be  developed  that  is 
consistent  with  the  available  data. 

2.  The  first  principles  electrostatic  spray  model  with  fixed  -a '  If}'  (verified  by 
this  work)  reveals  that  beyond  the  transition  mean  droplet  size  is  not  a 
function  of  fluid  properties.  It  is  solely  dependent  on  charge  density  level 
and,  is  thus,  a  purely  electronic  effect. 


26 


3.  Estimates,  discussed  below,  indicate  that  the  mean  electrostatic  energy  of  the 
droplet  surface  charge  along  the  B  branch  is  constant,  independent  of  droplet 
size  and  is  consistent  with  what  is  known  about  the  electron  crystallization 
process. 

The  Electron  Crystal 

One  of  the  more  perplexing  aspects  of  this  theoretical  work  is  the  predicted 
discontinuous  charging.  For  all  intents  and  purposes  the  trace  of  predicted  charging  level 
vs.  droplet  size  represents  a  first  order  phase  transition  (cf.  Figure  7).  This  phase 
transition  in  droplet  charging  level  is  explicable  within  the  context  of  Coulomb 
crystallization. 

Evidence  for  this  is  presented  in  Figure  15  where  the  ratio  of  electrostatic  to 
thermal  energy  for  individual  droplet  surface  electrons  is  shown  plotted  vs.  droplet  radius 
using  the  a',  0'  parameters  from  the  octoil  data  [21-23].  The  phase  transition  behavior  is 
clearly  evident.  As  droplet  radius  increases  the  energy  ratio  increases  reaching  a  peak  of 
136  prior  to  decreasing  along  the  single  solution  portion  of  the  plot.  This  value  is 
interestingly  close  to  the  critical  value  Tm  =  137,  noted  by  Brinkman,  et  al.  [28]  for  the 

formation  of  planar  electron  "crystals".  It  should  be  noted  that  a  similar  calculation  for 
Wood’s  metal  yields  a  peak  value  of  122  (cf.  the  table).  However,  this  value  is  based  on 
the  assumption  that  the  data  were  collected  at  the  melting  point.  The  actual  test 
temperature  is  uncertain. 

As  the  spray  droplet  size  is  increased,  the  energy  ratio  decreases  until  the  onset  of 
multiple  solutions  at  a  critical  droplet  radius  of  ~0.85/an.  It  is  the  central,  minimum 
energy  branch  to  which  actual  sprays  conform. 

Immediately  after  the  discontinuity  and  bifurcation  the  ratio  on  this  branch 
decreases  from  somewhat  above  the  asymptotic  level  of  437  to  slightly  below,  and  rises 
slowly  to  ultimately  attain  this  value  for  large  droplet  sizes.  The  virtual  constancy  of  the 
energy  ratio  is  apparent.  As  is  noted  on  the  plot,  the  mean  electrostatic  energy  of 
individual  surface  charges  is  {—a'/0')/6  or  lleV;  a  value  that  is  suggestively  close  to  the 
first  ionization  potential  of  hydrocarbons  [29,30].  By  comparison  it  requires  an  average 
energy  expenditure  of  -a’  /0'  =  66eV  to  transport  a  charge  to  the  surface  of  droplets  in 
the  plume. 

The  essential  feature  to  note  is  that  the  energy  ratio  is  predicted  to  change  abruptly 
from  a  value  that  is  below  the  putative  critical  value  of  137  to  a  level  that  is  pi  times  as 
high.  It  i:  this  behavior  that  prompts  the  suggestion  that  Coulomb  crystallization  is 
occurring  on  spray  droplet  surfaces  in  the  same  manner  that  it  is  known  to  occur  in  planar 
electron  clouds  that  can  be  formed  in  close  proximity  to  liquid  helium  surfaces. 


27 


ROOM  TEMPERATURE 
OCTOIL  SPRAY 


Figure  15:  DROPLET  RADIUS  •  Jim 


Wood’s  metal  sprays  are  sufficiently  similar  to  those  of  octoil  as  to  suggest  that  a 
universality  of  behavior  exists.  This  is  exhibited  in  Figure  16  where  the  ratio  of 
electrostatic  to  thermal  energy  per  surface  electron  is  normalized  by  the  asymptotic  value. 
This  normalization  eliminates  the  uncertainty  of  test  temperature  from  the  result  and 
reduces  the  ordinate  to  a  ratio  of  electrostatic  energies.  In  addition,  the  droplet  radius  has 
been  normalized  by  the  discontinuity  value.  The  detailed  similarity  of  the  relative 
behavioral  pattern  of  the  two  physically  disparate  liquids  is  apparent.  The  parameter 
listing  of  the  table  permits  an  absolute  comparison  to  be  made. 

Several  facts  are  certain:  1.  the  use  of  empirically  derived  values  for  a'  and  /?' 
(from  the  available  data  for  octoil  and  Wood’s  metal),  when  applied  to  the  statistical 
equilibrium  model  of  electrostatic  spraying,  predicts  discontinuous  spray  charge  behavior, 
2.  the  discontinuity  brackets  the  currently  accepted  energy  ratio  for  the  occurrence  of 
Coulombic  crystallization,  and  3.  the  asymptotic  spray  regime  corresponds  to  a  fixed 
energy  ratio  value  that  is  well  above  that  required  for  crystallization  to  occur.  We 
conclude  that  droplets  in  electrostatic  sprays  having  mean  sizes  larger  than  about  one 
micron  and  similarly  charged  isolated  droplets  are  encased  in  electron  crystals  and  that  the 
fluid  in  these  crystals  is  irrelevant  except  to  form  a  surface  upon  which  crystallization  can 
occur.  This  offers  an  opportunity  to  use  this  unique  form  of  matter  for  catalysts  or  in 
applications  such  as  oscillators  or  where  precise  control  of  the  state  of  electrons  is  of 
interest  and  importance. 


A.J.  Kelly 
Senior  Research  Engineer 


29 


NORMALIZED  TO  DISCONTINUITY  RADIUS 


DROPLET  SURFACE  CHARGE  PARAMETERS 


Pre-discontinuity  Peak 
Energy  ratio 
Droplet  radius  (pm) 
Charge  density  (m-2) 

Discontinuity 

Energy  ratio 
Droplet  radius  (pm) 
Charge  density  (m-2) 
Asymptotic  Energy  Ratio 

*  T  =  293°  K 
**  T  =  343®  K 


TABLE 


Fluid 


Octoi 1 


Wood's  Metal 


137* 
0.23 
6.2x10’ * 


122** 
0.020 
6.5x10’ * 


68  -  493 
0.87 

7 . 3x10’ 4  -  4.8x10’* 


65-415 

0.060 

1.1x10’*  -  6 


437 


373 


2x10’ • 
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A  REAL  TIME  TRACKER  CONTROLLED 
ELECTRODYNAMIC  BALANCE 

by 

K.Pradeep  Seneviratne 


ABSTRACT 


Implementation  of  a  'Electrodynamic  ‘Balance  to  levitate  electrostatically 
charged  large  liquid  droplets  is  presented.  ‘Ihe  levitator  consists  of  two 
horizontally  placed  electrodes  and  a  nozzle  to  eject  liquid  droplets.  ‘Mien  a  droplet 
is  ejected  to  the  space  Between  the  electrodes,  a  real  time  ‘Pyramid  ‘Tracker 
determines  the  position  of  the  center  of  the  droplet  ‘Picture  frames  of  the  droplet  is 
recorded  By  a  CCIV  camera,  and  trasfered  to  the  tracker  through  a  VCR,  A 
COMPAQ,  3 86  obtains  the  position  of  the  droplet  from  the  tracer  and  using 
position,  velosity,  and  acceleration  feedback,  determines  the  appropriate  voltage 
difference  required  to  stabalize  the  droplet.  ‘Ihe  voltage  diffemce  Between  the  two 
electrodes  is  constantly  changed  to  compensate  for  the  mass  loss  of  the  droplet  due 
to  evaporation.  ‘Ihe  ultimate  object  of  the  experiment  is  to  observe  sibling 
formation  characteristics  of  the  suspend  droplet  when  it  evaoporates  and  loose 
enogh  mass  to  collapse  in  to  many  droplets,  Results  of  the  experiment  verify  that 
the  tracker  is  fast  enough  to  calculate  the  center  of  the  droplet  up  to  four  times  as  it 
travels  across  the  view  of  the  camera.  Since  the  velosities  of  the  droplets  in  these 
studies  have  Been  comparable  to  that  of  ejection  velosities  of  droplets,  the  tracker  is 
fast  enough  for  die  COMPAQ  to  calculate  the  required  voltage  difference  Between 
the  electrodes.  ‘Ihe  results  of  these  tests  further  indicate  that,  when  a  fast  moving 
liquid  spray  is  ejected  to  the  space  Between  the  electrodes,  the  tracker  is  capable  of 
identifvng  the  number  of  pixels  contained  in  some  of  the  Bigger  droplets  of  the  spray. 
This  ability  of  the  tracker  is  very  usefull  in  studying  the  process  of  sibling 
formation. 
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INTRODUCTION* 

Levitation  of  electrostatically  charged  liquid  droplets  has  generated  a  great  deal 
of  interest  in  recent  years  because  of  its  possible  applications  in  material  processing 
and  its  potential  usefulness  in  the  understanding  of  properties  of  droplets  under  zero 
gravity  conditions.  Containerless  processing  of  materials  in  space  as  well  as  in 
ground-based  laboratories,  supercooling  and  superheating  of  droplets,  dynamics  of 
non-stationary  objects,  and  characteristics  of  evaporating  liquid  droplets  are  some  of 
the  applications  and  experiments  of  current  interest  that  will  require  an  efficient 
levitator. 

While  experiments  on  large  liquid  droplets  can  be  performed  free  from  physical 
contact  in  space  under  various  conditions,  the  same  experiments  are  difficult  to 
duplicate  on  ground  due  to  the  fact  that  the  size  of  the  drop  to  be  levitated  against 
one-gravity  is  limited  by  its  density  and  surface  tension.  Levitation  of  large  droplets, 
up  to  several  millimeters  in  diameter,  has  been  achieved  in  wind  tunnels  since  the  late 
1960's^,  but  such  levitated  drops  invariably  suffer  from  shape  distortion  due  to  air 
drag.  Presently,  there  are  two  common  methods  used  in  levitating  liquid  droplets 
without  a  wind  tunnel;  unfortunately,  both  methods  have  inherent  drawbacks  as  well. 
In  die  first  method,  a  droplet  is  suspended  by  forces  and  torques  generated  by  applying 


*  The  first  three  paragraphs  of  this  chapter  was  taken  from  reference  3:  a  paper 
written  by  the  author  and  another  student. 
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acoustic  waves  on  the  surrounding  medium.  It,  however,  has  the  disadvantage  of 
requiring  a  gas  medium  around  the  droplet  to  be  suspended.  Moreover,  a  droplet 
levitated  using  this  technique  suffers  shape  distortion  due  to  forces  generated  by  the 
acoustic  waves  ^  .  The  second  method  employs  a  RF  (radio  frequency)  field  to 
suspend  the  droplet.  This  method  can  only  be  used  on  electrically  conductive 
materials  and  thus  limits  the  number  of  liquids  which  can  be  used  in  conjunction  with 
such  a  levitator. 

The  method  of  levitation  this  paper  propose  is  to  suspend  the  liquid  drop  by 
first  inducing  an  electrostatic  charge  on  the  droplet,  and  then  launching  the  droplet  into 
an  electric  field  generated  by  a  pair  of  horizontal  parallel  electrodes.  The  position  of 
the  droplet  is  located  by  a  CCTV  (closed  circuit  television)  camera  and  controlled  by  a 
feedback  loop  to  a  microcomputer  connected  to  the  voltage  sources.  This  method  is 
free  from  the  disadvantages  of  the  above  mentioned  techniques;  therefore,  it  increases 
the  scope  of  applications  of  liquid  drop  levitators.  This  type  of  levitator  for  a  large 
droplet  against  1-g  conditions  has  only  been  achieved  before  in  the  Jet  Propulsion 
Laboratory  at  California  Institute  of  Technology  1®. 

The  method  of  levitation  of  large  liquid  droplets  this  paper  propose  was 
originally  tried  out  to  suspend  droplets  by  the  author  and  another  classmate^.  This 
attempt  was  unsuccessfull  due  to  two  reasons.  Firstly,  the  frame  grabber  we  used, 
though  was  able  to  digitize  frames  in  real  time,  the  IBM  PC-XT  we  used  was  not  able  to 
process  these  frames  in  real  time.  Thus,  our  PC-XT  was  not  able  to  respond  fast 
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enough  to  change  the  voltage  diffemce  between  the  two  electrodes  in  order  to  stabalize 
the  droplet  Secondly,  the  power  supplies  we  used  did  not  have  a  fast  enough  response 
time.  In  this  experiment  I  attempted  to  aviod  the  first  of  these  two  problems  by  using  a 
Pyramid  tracker  instead  of  the  frame  grabber  and  the  powerfull  COMPAQ  386 
computer  instead  of  a  IBM  PC-XT.  The  second  problem  will  be  avoided  simply  by 
building  a  very  fast  response  high  voltage  power  supply. 

Other  than  these  two  major  differences  and  few  other  modifications  in  some  of 
the  equipment,  I  will  be  using  the  same  levitator  system  for  this  experiment.  Thus, 
description  of  system  components  in  the  sections  3.0  through  3.6  of  chapter  3  was 
taken  directly  from  the  report  submitted  by  me  and  my  class  mate(  these  sections  were 
changed  slightly.  For  instance,  instead  of  the  IBM  PC-XT,  the  COMPAQ  386  and 
instead  of  the  frame  grabber  the  Pyramid  tracker  are  described  in  these  sections.)  The 
modifications  I  did  to  the  system  is  explained  in  the  section  3.7  and  the  fast  response 
power  supply  is  detailed  in  section  3.8. 

Chapter  4  details  charging  the  droplet,  ejecting  it  and  the  feedback  control 
sequence  to  stabalize  the  droplet.  Chapter  5  is  a  summary  of  the  experiment  detailing 
some  of  the  results  obtained  and  the  conclusions  derived  from  these  results.  Appendix 
A  describes  the  vision  system  of  the  system.  The  main  component  of  the  vision 
system-the  pyramid  tracker  is  detailed  in  this  section.  Capabilities  of  the  tracker  as 
well  as  the  advantages  of  using  the  pyramid  tracker  instead  of  a  conventional  frame 
grabber  are  described  in  section  6.1.  In  section  6.2  die  algorithm  used  to  calculate  the 
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center  of  die  droplet  is  explained.  Section  6.3  contains  the  software  programs  which 
provide  information  on  tracking  the  liquid  droplet. 

Much  time  was  spent  in  obtaining  the  necessary  equipment  and  understanding 
and  reprograming  the  pyramid  traker.  Now  that  the  tracker  is  completely 
reprogramed  to  suit  this  experiment  and  most  of  the  equipment  is  built,  the 
electrodynamic  balance  could  be  assembled  and  tested. 


DESCRIPTION  AND  DESIGN  OF  EQUIPMENT 


3.1  Background: 

An  evaporating  liquid  droplet  suspended  in  an  electrostatic  field  is  in 
equilibrium  under  two  forces;  namely,  the  gravitational  force  acting  on  the  droplet 
and  the  electrostatic  force  provided  by  the  two  horizontal  parallel  electrodes.  At 
equilibrium,  we  have  the  following  equation: 

V(t)Qld  =  m(t)g  (i) 

where:  V(t)  =  Variational  voltage  difference  between  the  electrodes. 

Q  =  Charge  carried  by  the  droplet. 
d  =  Distance  between  the  two  electrodes. 
m(t)  =  Variational  mass  of  the  droplet. 
g  =  Acceleration  of  gravity. 

Due  to  the  loss  of  mass  from  evaporation,  m(t)  will  be  a  function  of  time.  The 
rate  of  mass  loss  of  the  droplet  depends  on  such  external  factors  as  temperature  and 
humidity  of  the  surrounding  medium.  When  m(t)  changes,  the  voltage  difference 
V(t)  has  to  be  changed  accordingly  to  provide  a  corresponding  compensating 
electrostatic  force.  Initial  hand  control  of  V(t)  has  been  attempted  unsuccessfully  at 
the  Electric  Propulsion  Laboratory  of  Princeton  University.  In  fact,  due  to- the  high 
rate  of  mass  1ms  as  well  as  the  different  V(t)'s  required  for  the  launching  and 
levitating  phases  of  the  experiment,  it  is  necessary  that  the  voltage  difference  be 
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continuously  monitored  and  varied  by  means  of  a  computer  controlled  feedback 
system. 

32  Overall  System  Configuration: 

The  overall  experiment  set  up  is  shown  in  Figure  1  with  die  main  controller 
and  CPU  of  die  system  being  a  COMPAQ  386.  The  input  voltage,  V(t),  to  the 
levitator  is  determined  by  the  microprocessor  through  continuous  tracking  of  die 
position  of  the  droplet,  using  an  image  acquisition  system  consisting  of  a  CCTV 
camera,  amonitor,  a  VCR,  and  a  Pyramid  tracker  built  at  David  Saranoff  Research 
Center.  The  digitized  image  by  the  tracker  can  be  displayed  on  the  monitor  in  real 
time  for  focusing  and  positioning  purposes  during  the  set  up  phase  of  the 
experiment.  Once  V(t)  is  calculated  by  the  COMPAQ  386 ,  it  is  send  to  the  system's 
high  voltage  power  supplies  via  an  IBM  PC  parallel  interface  board  connected  to 
digital-to-analog  converters  and  inverting  operational  amplifiers.  The  high  voltage 
power  supplies  are  monitored  by  a  NICOLET  oscilloscope,  which  can  dump 
incoming  signal  to  a  Hewlett  Packard  microcomputer  connected  to  a  plotter  for 
voltage  vs.  time  plots.  Finally,  the  appropriate  voltage  V(t)  is  send  to  die  electrodes 
mounted  on  a  multi-purpose  room  temperature  chamber. 

33  Microprocessor  To  Power  Source  Interface: 

The  circuit  diagram  of  the  microprocessor  to  high  voltage  power  supply 
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interface  is  shown  in  Figure  2.  Since  two  different  power  supplies  are  used  for  die 
two  electrodes,  the  output  byte  of  the  parallel  interface  is  split  into  two  4-bit  nibbles 
feeding  into  two  MC1408  8-bit  D-to-A  converters.  The  4  data  bits  are  connected  to 
bits  3, 4, 5,  and  6  of  the  D-to-A,  while  bits  0, 1, 2,  and  7  of  the  D-to-A  are  jumpered 
to  ground.  This  allows  each  power  supply  to  have  16  possible  input  control  voltage 
levels.  Each  D-to-A  converter  has  a  3k  Q  resistor  connected  between  it  and  die  +5V 
power  supply  to  provide  the  necessary  drive  current  and  a  1000  pF  capacitor 
connected  between  it  and  the  -1SV  power  supply  for  stabilization.  Die  output  of 
D-to-A  converters  are  connected  to  TL081  inverting  differential  operational 
amplifiers,  since  each  D-to-A  acts  as  a  current  source.  The  timing  of  the  parallel 
port  is  controlled  via  software  delay  loops,  because  D-to-A  converters  are  dumb 
peripherals  that  do  not  take  advantage  of  the  status  and  control  ports  provided  on  the 
parallel  communication  board.  The  only  line,  other  than  data  lines,  that  is  used  for 
interfacing  is  the  BUSY  line,  which  is  jumpered  to  the  +5V  power  supply  so  it  will 
never  be  asserted  during  normal  operation. 

3.4  Voltage  Monitoring  System: 

To  observe  the  high  voltage  outputs  through  die  NICOLET,  special  voltage 
divider  circuits  had  to  be  made  because  the  maximum  allowable  input  voltage  to  die 
MCOLET  is  40V.  Since  the  outputs  of  the  high  voltage  devices  have  to  be  fed  into 
die  electrodes  as  well,  two  T-junctions  were  made  to  interface  each  of  the  voltage 
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devices  with  the  NICOLET  and  the  electrodes.  The  T-junctions  themselves  are  well 
covered  with  epoxy  to  avoid  high  voltage  discharges.  The  connection  between  die 
high  voltage  device  and  the  T-junction,  as  well  as  the  connection  between  the  voltage 
divider  circuit  and  the  T-junction,  are  both  temporary  connections,  so  alligator 
clamps  are  used. 

It  should  be  noted  that  even  after  taking  great  care  in  making  these  connections 
and  grounding  the  table  on  which  the  high  voltage  devices  are  placed,  it  is  not 
possible  to  avoid  minor  electric  shocks  when  we  touch  the  table  or  the 
room-temperature  chamber.  Only  after  suspending  all  high  voltage  cables  in  air 
using  clamps  and  keeping  the  cables  far  away  from  each  other  was  this  experiment 
performed  safely. 

The  voltage  dividers  used  in  our  experiment  consist  of  a  400MO  resistor  in 
conjunction  with  a  lOkQ  resistor  to  produce  an  approximate  1:40,000  voltage 
divider  for  the  power  supply  which  control  the  voltage  of  the  bottom  electrode,  and 
a  100M11  resistor  in  conjunction  with  a  lOOkft  resistor  to  produce  an  approximate 
1 : 1,000  voltage  divider  for  the  power  supply  which  control  die  voltage  of  die  top 
electrode.  The  400M&  and  the  lOOMft  resistors  are  high  quality  ceramic  resistors. 
These  resistors  are  necessary  since  ordinary  resistors  cannot  be  guaranteed  to  retain 
their  specified  ohmic  values  when  subjected  to  high  temperatures  caused  by  drastic 
drops  in  high  voltages.  The  values  lOOkfl  and  10kft  are  chosen  since  these  are 
negligible  when  compared  to  the  input  impedance  of  the  NICOLET.  The  voltage 
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outputs  of  the  high  voltage  devices  can  then  be  observed  by  connecting  the  voltage 
dividers  to  the  NICOLET  as  shown  in  fiqure  1,  and  hardcopies  of  the  voltage  vs. 
time  plot  can  be  obtain  by  dumping  data  from  the  NICOLET  to  a  Hewlett  Packard 
microcomputer  connected  to  a  HP  plotter. 

3J  Electrode  Design  And  Specifications: 

hi  this  experiment,  the  concern  is  stabilization  of  the  liquid  drop  along  the 
vertical  axis,  since  the  electric  field  created  by  the  parallel  electrodes  can  be  utilized 
to  confine  the  path  of  the  droplet  in  lateral  directions  to  a  reasonable  degree.  This  is 
easily  accomplished  by  choosing  properly  shaped  electrodes.  As  a  first  attempt,  a 
circular  convex  upper  electrode  and  a  matching  concave  lower  electrode  was  used. 
The  symmetry  of  the  electric  field  and  the  strong  upward  electric  force  provide  die 
necessary  centering  force,  as  the  levitated  object  follows  the  horizontal  potential 
well  along  the  axis  of  symmetry.  This  design  was  also  used  in  an  earlier  electrostatic 
positioner  at  JPL^.  As  a  direct  consequence  of  equation  (i)  given  earlier,  the 
centering  force  created  by  such  a  pair  of  electrodes  is  proportional  to  (V QldRjX, 
where  X  is  die  separation  between  die  object  and  the  axis  of  symmetry,  and  R  is  die 
curvature  radius  of  the  electrodes.  This  type  of  electrode  configuration,  however, 
requires  a  large  constant  gravitational  force;  therefore,  it  can  only  be  operated  in 
ground-based  laboratories. 

The  electrodes  actually  used  in  the  experiment  are  shown  in  Figure  3.  It 
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follows  basically  the  same  design,  except  for  dimensions,  as  that  used  in  JPL^,  and 
it  was  built  at  the  Electric  Propulsion  Laboratory  of  Princeton  University.  The 
electrodes  are  a  pair  of  cylindrical  aluminum  disks.  The  top  ring  is  tapered  to 
provide  the  necessary  centering  force,  and  both  rings  are  connected  to  the 
COMPAQ  386  controlled  high  voltage  power  supplies.  The  electrodes  are 
physically  mounted  on  plexiglas  boards  that  fit  into  the  multi-purpose  room 
temperature  chamber  used  in  this  experiment.  An  injection  nozzle  can  be  mounted 
into  the  chamber  through  a  0.1S  inch  opening  at  the  center  of  the  bottom  electrode. 

3.6  Multi-Purpose  Room  Temperature  Chamber: 

The  plexiglas  levitation  chamber  used  in  Hus  experiment  was  also  constructed 
here  at  the  Electric  Propulsion  Laboratory.  The  electrodes  described  above  are 
mounted  at  the  centers  of  both  die  top  and  the  bottom  plates  of  the  chamber,  and  the 
chamber  itself  has  an  inside  dimension  of  6"  x  6"  x  1-13/32".  At  the  present  time, 
the  multi-purpose  room  temperature  chamber  is  designed  only  for  the  use  of 
suspending  a  charged  droplet;  however,  it  can  be  easily  modified  to  have  the 
capability  of  inducing  levitated  droplet  rotation  and/or  oscillation,  like  the  final 
version  used  in  the  experiment  at  JPL.  The  basic  method  to  induce  specific  motion 
on  an  electrostatically  levitated  droplet  is  to  apply  either  an  acoustic  torque^  or  an 
oscillating  acoustic  pressure1 .  One  possible  approach  is  to  have  two  acoustic 
transducers  mounted  on  adjacent  side  walls  of  the  levitation  chamber  so  that 
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mutually  orthogonal  acoustic  standing  waves  can  be  created  within  die  chamber1®. 
These  acoustic  transducers  can  be  driven  by  amplified  outputs  from  two  phase  and 
frequency  locked  synthesizers.  Mathematically,  we  know  that  the  acoustic  torque  is 
given  by  the  following  equation,  according  to  Busse  and  Wang^: 

T  =  (312) lri(PxPy/2pc2)A  sin  <p •  (ii) 

where:  =  Viscous  length. 

Px  =  Pressure  amplitude  of  the  standing  wave  in  the  x-direction. 
Py  =  Pressure  amplitude  of  die  standing  wave  in  die  y-direction. 

p  =  Density  of  the  medium. 

c  =  Velocity  of  sound. 

A  =  Surface  area  of  the  sphere. 

<p  •  =  Relative  phase  angle  between  the  orthogonal  waves. 

Clearly,  a  varying  level  of  acoustic  torque  can  be  produced  by  either  holding 
the  relative  phase  fixed  while  varying  the  sound  pressure  levels  or  holding  the  level 
of  sound  pressure  constant  and  adjusting  the  relative  phase.  According  to 
Annamalai,  Trinh,  and  Wang1,  we  can  induce  oscillation  cm  a  non-rotating  liquid 
droplet  by  fixing  die  relative  phase  at  0*  and  sinusoidally  modulating  die  ampliude 
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of  die  sound  wave.  Similarly,  both  rotation  and  oscillation  can  be  induced  at  the 
same  time  by  applying  a  non-zero  relative  phase  with  amplitude  modulation.  Such  a 
modified  chamber  can  be  used  to  carry  out  a  range  of  experiments  on  charged  drop 
dynamics,  since  surface  tension  and  viscosity  of  a  liquid  droplet  can  be  calculated 
from  observations  of  the  oscillating  droplet  at  its  natural  frequency.  One 
experiment  of  interest  is  to  observe  the  behavior  of  a  large  evaporating  electrically 
charged  liquid  droplet  near  the  threshold  of  Rayleigh's  criteria  for  stability. 

3.7  Modifications  of  the  System: 

During  the  experiment  which  I  and  my  class  mate  carried  out  ( see 
introduction)  it  was  observed  that  after  operating  the  levitator  for  many  hours,  the 
electrodes  can  no  longer  be  predicted  to  operate  in  a  consistent  manner.  To  avoid 
this  problem  and  make  the  levitaor  more  reliable  and  efficient,  I  did  the  following 
modifications  to  the  levitaor. 
a)  The  reservoir  of  liquid: 

When  the  levitator  is  operated  for  a  long  time,  a  considerable  amount  of 
liquid  from  the  reservoir  is  drawn  out.  After  a  fair  amount  of  liquid  has  flown  out 
of  the  reservoir,  the  internal  pressure  acting  on  the  liquid  at  the  tip  of  the  needle  will 
have  a  different  value  from  before;  thus  resulting  in  the  need  for  a  new  ejection 
voltage.  The  solution  to  this  problem  is  to  have  a  reservoir  with  a  near  constant 
pressure  near  the  tip  of  the  needle  at  all  times.  Thus,  in  place  of  the  fixed  reservoir. 
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a  reservoir  of  liquid  which  could  slide  approximately  2"  in  the  vertical  direction  was 
installed.  Now,  by  sliding  the  reservoir  the  liquid  pressure  on  die  tip  of  the  ne  die 
can  be  kept  constant. 

b.)  The  position  of  the  needle  : 

One  of  the  drawbacks  that  was  affecting  the  levitator  was  the  relatively 
unstable  stand  used  to  hold  the  needle  in  position.  "After  carrying  on  the  experiment 
for  a  few  minutes,  some  liquid  would  accumulate  on  the  bottom  electrode.  When  an 
attempt  is  made  to  wipe  out  this  liquid  in  order  to  preserve  the  surface  shape 
characteristics  of  the  electrode,  often  die  needle  would  be  displaced  slightly;  hence 
varying  a  dimension  that  should  be  kept  constant  if  consistent  results  are  to  be 
expected."^ 

In  order  to  solve  this  problem,  die  needle  positioning  system  was  designed  in 
the  late  part  of  the  above  mentioned  experiment.  However  die  method  implemented 
proved  to  be  of  no  great  promise.  Thus,  to  obtain  a  stable  needle,  a  clamping 
mechanism  which  tightly  position  die  glass  tube  which  leads  to  the  needle  was 
employed.  Even  when  considerable  force  is  employed  on  die  glass  tube  .this  design 
seems  to  be  able  to  maintain  the  needle  tip  in  a  stable  position. 

c)  The  electrodes: 
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aluminum  rings  which  supported  two  central  metal  coated  glass  disks.  It  was 
observed  that  after  some  time  this  metal  coat  would  wear  away  leaving  behind 
electrodes  of  unsuitable  surface  characteristics.  To  avoid  this  problem,  alumunum 
disks  were  placed  in  place  of  the  glass  disks. 

d)  The  droplet  ejection  nozzle: 

I  plan  to  cover  the  tip  of  the  ejection  nozzle  with  a  non  wetting  cloth  so  that 
surface  tension  force  between  the  droplet  and  the  needle  could  be  kept  to  a  minimum. 
Such  an  arrangement  would  reduce  the  voltage  necessary  to  eject  a  droplet  form  the 
needle  and  consequently  the  impotance  of  the  settling  time  of  the  power  supply  would 
be  reduced.  In  addition  to  that,  to  further  reduce  die  surface  tension  between  water 
and  the  needle,  small  amount  of  a  liquid  which  would  reduce  the  surface  tension  of 
water  will  be  added  to  the  reservoir. 

e.)  The  voltage  difference  between  the  electrodes: 

In  the  present  configuration  of  the  electrodes,  the  top  electrode  is  grounded 
while  the  voltage  of  the  bottom  electrode  is  kept  at  an  appropriate  negative  value.  It 
can  be  shown  that  under  this  arrangement,  negatively  charged  droplets  ejected  from 
the  nozzle,  could  experience  a  non  converging  electric  filed  By  grounding  the  top 
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the  nozzle,  could  experience  a  non  converging  electric  filed.  By  grounding  the  top 
electrode  and  placing  an  appropriate  po&sitive  voltage  on  the  bottom  electrode  this 
probelm  could  be  avoided.  However  the  droplets  ejected  under  this  arrangement 
will  be  positively  charged. 

3.8  The  Fast  Response  High  Voltage  Power  Supply : 

The  circuit  configuration  of  the  high  voltage  power  supply  is  illustrated  in 
figure  4.  The  Hitachi  flyback  transformer,  produces  a  high  voltage  output 
corresponding  to  the  input  voltage  at  the  low  voltage  section  of  the  power  supply,  as 
illustrated  in  graph  1.  There  are  many  advantages  of  using  a  Hitachi  flyback 
trasformer  to  generate  the  necessary  high  voltage  for  the  electrodes  of  the  levitator. 
The  most  important  of  these  is  that  this  transformer  is  capable  of  producing  an 
output  of  0  to  25kV  to  a  corresponding  input  voltage  with  a  transient  time  of 
approximately  1ms.  As  seen  from  graph  1  this  output  is  linear  with  the  input 
voltage. 

The  low  voltage  section  of  this  unit  was  seperated  from  the  high  voltage 
section  by  enclosing  the  high  voltage  section  in  a  metal  box.  The  high  voltage 
transformer  is  provided  with  a  square  wave  input.  The  output  of  the  transformer 
depends  cm  the  magnitude  and  the  frequency  of  this  square  wave.  Since  we  need  to 
have  the  dynamic  range  of  the  high  voltage  output  as  large  as  possible  the  frequency 
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of  the  square  wave  provided  to  the  transformer  has  to  be  the  resonant  frequency  of 
the  transformer.  The  Hitachi  flyback  trasformer  resonates  around  55kHz,  and  thus, 
the  frequency  of  the  square  wave  generated  by  the  NEC  555  oscillator  is  adjusted  to 
55kHz. 

The  amplitude  of  the  square  wave  determines  the  output  voltage  of  the 
transformer.  Since  the  output  voltage  of  the  trasformer  has  to  be  a  function  of  die 
input  voltage  to  the  low  voltage  section  of  the  power  supply,  ( This  input  comes 
frame  the  output  of  the  circuit  of  figure  3.  The  output  of  the  circuit  of  figure  3  is 
controled  by  the  COMPAQ  386.)  the  amplitude  of  the  square  wave  has  to  be 
controled  by  this  input  voltage.  The  input  voltage  to  the  low  voltage  section  of  the 
power  supply  is  first  amplified  and  properly  adjusted  by  means  of  a  buffer  amplifier 
and  a  level  adjuster.  This  properly  scaled  input  voltage  when  applied  across  the 
clamping  diode  1N916,  controls  the  amplitude  of  the  square  wave.  This  controled 
square  wave  is  combined  with  a  voltage  stabalizing  feedback  signal  at  the 
comparator  TL074.  The  output  of  the  comparator  is  fed  to  the  primary  winding  of 
the  flyback  transformer  to  generate  the  corresponding  output  voltage.  This  output 
voltage  of  the  trasformer  when  properly  half-wave  rectified  and  filtered  should 
produce  a  positive  dc  voltage  with  a  1  %  55-kHz  ripple.1 1 

In  addition  to  the  elements  shown  in  the  circuit,  to  obtain  a  symmetric  square 
wave,  the  values  of  RA  and  Rfi  and  c  have  to  be  properly  chosen.  According  to  plot 


2  to  obtain  a  square  wave  of  55kHz  frequency,  die  values  of  RA,  Rfi,  and  c  have  to  be 
such  that^, 

RA  +  2RB  =  10kft. 

In  addition  since  th  and  tL  have  to  satisfy?; 

V*L  =  V<Ra  +  V 

Ra  ;>  lkfi. 

and; 

c  ~  .003  pF 

Solving  these  equations  RA  and  Rg  were  taken  to  be  Ak£i  and  2kQ  respectively. 


17 


STABILIZATION  OF  THE  DROPLET 


4.1  Charging  and  Ejecting  the  Droplet: 

Before  the  start  of  the  experiment,  we  have  a  column  of  liquid  located  at  the 
lower  section  of  the  injection  nozzle.  This  column  is  then  subject  to  an  electrostatic 
field  produced  by  the  two  electrodes,  with  the  bottom  electrode  being  negatively 
charged  with  respect  to  the  top  electrode.  Since  the  liquid  column  is  in  contact  with  die 
reservoir,  the  relatively  positive  top  electrode  will  cause  the  liquid  column  to  be 
induced  with  a  negative  charge.  The  mutual  attraction  between  the  top  electrode  and 
the  negatively  charged  liquid  column,  as  well  as  the  liquid  pressure  acting  on  the 
column  produced  by  a  reservoir  located  at  a  higher  level  than  the  tip  of  the  nozzle,  will 
cause  the  column  to  flow  toward  the  tip  of  the  needle. 

At  the  tip  of  the  needle,  the  first  speck  of  liquid  which  attempts  to  eject  itself  out 
from  the  nozzle  toward  the  top  electrode  will  be  stopped  from  doing  so  by  the  surface 
tension  between  the  needle  and  itself.  Thus,  a  small  amount  of  liquid  will  accumulate  at 
die  top  of  the  ejection  needle.  This  droplet  of  liquid  will  gather  more  and  more 
negative  charge  due  to  the  influence  of  the  positive  top  electrode.  Invariably,  such  an 
accumulation  of  charge  will  lead  to  more  liquid  being  pushed  to  the  top  of  the  needle. 
This  charge  will  cause  an  ever  increasing  electrostatic  force  to  push  the  now  bigger 
and  highly  charged  droplet  upwards.  When  this  force  is  large  enough  to  overcome  the 
surface  tension  between  the  droplet  and  the  needle,  a  spherical  droplet  will  be  observed 


18 


to  leave  the  tip  of  the  needle. 

As  mentioned  before,  I  and  my  class  mate  tryed  out  numerous  control 
algorithms  in  an  earlier  experimnt  to  stop  the  droplet  from  reaching  the  top  electrode 
after  its  ejection  from  the  needle.  During  that  experiment  we  found  a  particular 
algorithm  that  is  capable  of  stopping  droplets  from  reaching  the  top  electrode.  That 
algorithm  is  as  follows3: 

(1)  Provide  a  control  voltage  of  known  magnitude  to  the  bottomm  electrode. 
Provide  the  control  voltage  long  enough  so  that  the  output  of  the  power 
supply  will  reach  a  constant  value  and  will  remain  at  that  value  for  a  few 
seconds.  (The  method  we  used  to  observe  the  outputs  of  the  high  voltage 
devices  is  outlined  in  the  section  Voltage  Monitoring  System.) 

(2)  Increase  the  control  voltage  to  the  power  supply.  Let  the  computer  run 
through  a  loop  of  few  useless  instructions  to  allow  the  power  supply  to  reach 
an  ejection  voltage.  (Note  that  once  a  control  voltage  is  fed  into  the  power 
supply,  even  if  die  computer  goes  through  a  wait  loop,  the  control  voltage  to 
the  power  supply  will  remain  the  same,  since  data  is  latched  into  die  parallel 
port  of  die  PC.  To  change  the  control  voltage,  another  control  voltage  must 
be  explicidy  fed  into  the  power  supply.) 

(3)  Provide  a  zero  control  voltage  to  the  power  supply. 

(4)  While  steps  1 X  and  3  are  being  executed,  apply  a  zero  voltage  to  the  top 
electrode.  As  soon  as  step  3  is  executed,  apply  a  large  control  voltage  to  the 
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top  electrode. 

(5)  Wait  for  a  few  seconds  and  then  provide  a  zero  control  voltage  to  the  top 
electrode. 

(6)  Repeat  steps  1  through  6. 

The  function  of  the  individual  steps  in  this  algorithem  is  as  follows: 

(1)  Provide  voltage  large  enough  to  bring  up  a  droplet  to  the  needle  tip  and 
provide  enough  charge  to  that  droplet  so  that  if  a  higher  voltage  is  applied  (in 
the  next  step),  the  droplet  will  be  ejected  from  the  needle.  It  is  important  not 
to  provide  too  big  a  voltage  in  this  step  since  we  do  not  want  the  droplet  to  be 
ejected  while  this  step  is  being  carried  out  We  had  to  experiment  with  the 
levitator  for  a  while  before  selecting  such  a  control  voltage. 

(2)  Provide  a  higher  control  voltage  so  that  the  droplet  will  be  ejected.  It  is 
important  that  the  time  spent  on  the  wait  loops  is  not  too  long,  for  what  we 
want  is  to  provide  just  the  ejection  voltage  and  no  more. 

(3)  Switch  off  the  voltage  supply  to  the  bottom  electrode  so  that  the  voltage  on 
this  electrode  will  decay  to  zero. 

(4)  While  die  above  steps  are  being  performed,  the  voltage  on  die  top  electrode 
should  be  zero.  Step  3  will  eject  a  droplet  Thus,  by  supplying  a  large 
control  voltage  to  the  power  supply  immediately  after  step  3,  we  will  repel 
die  ejected  droplet  from  the  top  electrode. 


20 


(5)  Wait  to  make  sure  that  the  droplet  is  stopped  and  shut  off  die  voltage  to  the 
top  electrode. 

(6)  Start  again  to  eject  another  droplet. 

While  executing  this  program,  we  observed  that  after  ejection,  the  droplet  would 
travel  up,  reach  zero  velocity,  and  then  fall  back.  I  hope  to  use  this  same  algorithm 
to  charge,  launch  and  stop  the  droplet  from  reaching  the  top  electrode. 

42  Feedback  Control  Sequence: 

Since  there  are  4  bits  of  data,  which  can  easily  be  expanded  into  8  bits  by 
sending  out  separate  controlling  signals  to  the  top  and  the  bottom  electrodes,  for 
each  power  supply  at  disposal,  a  range  of  control  algorithms  are  available  for 
stabilizing  the  droplet.  The  obvious  extension  from  a  bang-bang  control  is  a  simple 
proportional  control  scheme,  where  variable  control  voltages  are  sent  to  the 
electrodes,  depending  on  the  actual  position  of  the  drop  within  the  frame  rather  than 
its  presence  in  a  particular  region.  Furthermore,  if  the  velocity  of  the  droplet  can  be 
measured  based  on  the  previous  position  of  the  droplet,  a  straight  forward  PID,  i.e. 
proportional-integral -derivative,  control  algorithm  can  be  implemented  to  model  a 
simple  second-order  linear  feedback  system.  The  relationship  in  discrete  form  that 
governs  such  a  control  algorithm  is  as  follows: 
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n 

S[n]  =  GpE[n]  +  G/  (I  E[j]  )  +  GJE[n] -E[n-1])  (iii) 

j=0 

where:  S[n]  =  Output  at  n^1  sample  time. 

E[n]  =  Ptf  -P[n]  =  Position  error. 

Pj  »  Desired  position  in  pixels. 

P[n]  =  Actual  position  in  pixels  at  n^1  sample  time. 


Gp  =  Proportional  gain. 

G/  =*  Integral  gain. 

Gj  as  Derivative  gain. 

One  advantage  that  the  PID  control  does  have  over  the  simple  bang-bang  control  is 
greater  flexibility,  where  various  design  tools,  such  as  root-locus,  Nyquist,  or  Bode 
plots  can  be  used  to  pick  the  desired  gains  for  different  situations. 
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SUMMARY  AND  CONCLUSION: 


Success  of  this  project  depends  on  two  major  capabilities  of  the  levitator  system. 
The  first  is  the  capability  of  the  high  voltage  power  source  to  have  a  fast  enough  decay 
time  such  that  the  droplet  will  not  hit  the  top  electrode.  This  means  that  the  power 
supply  has  to  have  a  very  fast  transient  time.  The  second  is  die  capability  of  the  image 
processor  of  the  system  to  process  frames  fast  enough  so  that  die  COMPAQ  386  could 
control  the  voltage  difference  between  the  two  plates  fast  enough  to  stabalize  the 
droplet. 

The  first  of  these  two  problems  will  be  avoided,  as  mentioned  before  by  building 
a  fast  response  high  voltage  unit.  If  the  transient  time  of  the  power  supply  is  within  its 
guranteed  range  of  1ms,  then  this  power  supply  would  certainly  be  fast  enough  for 
our  purpose.  The  low  voltage  section  and  the  control  section  of  the  power  supply  has 
been  already  built.  Building  the  high  voltage  section  of  the  power  supply  has  been 
halted  due  to  the  delay  in  obtaining  a  high  voltage  transformer  and  lack  of  vital 
information  regarding  the  outlets  of  this  high  voltage  tranformer.  Once  this 
information  is  available  the  power  supply  can  be  built  and  tested  for  its  transient  time. 

One  of  the  major  major  problems  of  the  above  mentioned  experiment  carried 
out  by  me  and  my  classmate,  was  that  once  die  droplet  is  stopped  by  the  top  electrode 
after  its  ejection,  the  droplet  showed  a  lateral  movement  in  its  falling  trajectory.  The 
reason  for  this  movement,  we  then  believed,  was  that  the  top  electrode  causing  too 
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great  a  damping  force  to  the  oncoming  droplet.  Since  we  were  able  to  provide  only  16 
different  control  voltages,  at  that  time,  to  the  power  supply,  we  could  not  achieve  the 
exact  damping  force  that  was  necessary  to  smoothly  turn  back  the  oncoming  droplet, 
however,  if  the  power  supply  which  I  am  in  the  process  of  building  prove  to  have  a 
fast  enough  transient  time,  there  is  no  need  to  control  the  voltage  of  the  top  electrode. 
In  this  case  not  only  can  we  avoid  the  problem  of  top  electrode  diverting  the  droplet, 
but  we  will  also  be  able  to  control  the  bottom  electrode  wirh  256  different  control 
voltages  levels. 

The  second  of  the  major  concerns  of  the  project  has  been  avoided  by  using  the 
Pyramid  tracker  as  the  image  processor  of  the  system.  The  capabilities  of  the  tracker 
and  how  these  capabilities  are  modified  so  that  the  tracker  functions  as  a  real  time 
frame  grabber  are  detailed  in  Appendix  A.  Results  of  some  of  the  tests  performed  to 
determine  the  speed  at  which  the  tracker  could  calculate  the  center  of  the  droplet  are 
outlined  below. 

In  these  experiments  a  reservoir  of  liquid  was  placed  at  a  considerable  height 
above  the  view  of  the  camera.  The  liquid  in  the  reservoir  is  water-colored  by  mixing 
with  black  ink.  A  small  hole  was  punched  in  the  bottom  of  die  reservoir  so  that  very 
small  liquid  droplets  would  emanate  from  this  hole  and  fall  down  passing  die  view  of 
die  camera.  It  should  be  noted  that  the  vertical  distance  between  die  camera  and  die 
bottom  of  the  reservoir  is  such  that  the  droplet  will  pass  the  view  of  the  camera  at  a 
considerable  velocity.  This  was  necessary  if  we  are  to  determine  whether  the  tracker 
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is  capable  of  detecting  a  fast  moving  droplet.  ( Note  that  the  tracker  has  to  detect  die 
droplet  when  the  droplet  is  ejected,  and  that  at  the  time  of  ejection  the  droplet  has  a 
considerable  velocity.) 

The  results  of  the  experiment  indicated  that  the  tracker  was  capable  of 
calculating  the  center  of  the  droplet  upto  four  times  while  the  droplet  was  moving 
across  the  view  of  the  camera.  Note  that  for  velocity  feedback  control  of  the  droplet 
as  explained  in  chapter  3,  the  tracker  has  to  detect  the  droplet  only  2  times.  Thus,  the 
tracker  is  certainly  fast  enough  for  the  purpose  of  this  experiment. 

One  of  the  other  experiments  carried  out  to  determine  die  capabilities  of  the 
tracker  is  as  follows.  The  reservoir  of  liquid  was  brought  closer  to  the  camera  and  a 
metal  needle  was  connected  to  the  reservoir  of  liquid  by  a  glass  tube.  The  glass  tube  is 
also  connected  to  a  high  voltage  power  supply.  Thus,  when  the  power  supply  is  turned 
on  it  will  inject  charge  to  the  liquid  resulting  in  a  spray  of  liquid  at  the  mettle  needle. 
Depending  on  the  output  voltage  of  the  power  supply  and  the  diameter  of  the  needle 
the  size  distribution  of  the  droplets  in  the  spray  could  be  varied.  The  tracker  was 
initialized  at  the  begining  of  the  experiment  such  that  the  frame  obtained  from  the 
camera  contains  no  pixel  above  a  prescribed  threshold  value.  (  The  process  of  doing 
this  is  explained  in  Appendix  A)  It  is  interesting  to  note  that  when  the  needle  was 
spraying  the  liquid  in  the  view  of  the  camera  the  tracker  was  able  to  identify  die  total 
number  of  pixels  which  are  above  the  threshold  value.  This  capability  of  the  tracker  is 
vital  in  observing  the  sibling  formation  characteristics  of  liquid  droplets.  For,  by 
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programming  the  tracker  to  group  all  the  pixels  which  are  adjacent  to  each  other  and 
which  are  above  the  threshold  value  we  can  determine  the  number  of  droplets  and  the 
size  distribution  of  these  droplets  when  an  evaporating  droplet  collapses  into  many 
droplets. 

Note  that  to  calculate  the  center  of  a  fast  moving  droplet  the  tracker  has  to 
calculate  the  center  of  the  droplet  at  least  twice  and  thus  be  programmed  to  operate  in 
Pyramid  levels  1  or  2.  ( See  Appendix)  The  need  in  this  application  is  then,  to  have  a 
veiy  fast  image  processor  at  the  expense  of  great  resolution.  To  calculate  the  number 
of  siblings  formed  when  a  bigger  droplet  collapses,  the  tracker  has  to  be  able  to 
identify  whether  a  group  of  pixels  is  separated  from  another  group  of  pixels  and  thus 
have  to  be  programmed  to  operate  in  Pyramid  levels  5  or  6.  The  need  here  is  to  have 
a  image  processing  system  of  high  resolution  at  the  expense  of  speed.  Note  that  I  was 
able  to  program  the  same  tracker  to  switch  between  speed  and  resolution.  In  this 
respect  perhaps  the  Pyramid  tracker  is  the  only  single  image  processor  which  is 
capable  of  successfully  stabalizing  a  droplet  and  obtain  information  on  sibling 
formation. 

The  effectiveness  of  the  modifications  done  to  the  electrodes,  reservoir,  and  the 
injection  needle  will  be  revealed  only  when  the  actual  experiment  is  carried  out. 
Much  time  was  spent  in  trying  to  obtain  the  necessary  equipment,  understand  the 
Pyramid  tracker  and  reprograming  it  to  suit  the  purpose  of  this  experiment  Now  that 
the  tracker  has  been  completely  modified  and  tested  for  its  capabilities,  with  positive 
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results,  once  the  high  voltage  power  supply  is  built  the  actual  experiment  of 
attempting  to  suspend  large  liquid  droplets  could  be  carried  out. 
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APPENDIX  A  -  THE  VISION  SYSTEM 


6.1  THE  PYRAMID  TRACKER 

The  Pyramid  Tracker  is  a  real  time  frame  grabber  which  can  be  programmed  to 
operate  in  6  "Pyramid  levels".  A  Pyramid  level  identifies  a  frame  transferred  by  a 
camera  as  a  picture  consisting  of  following  number  of  pixels. 

Pyramid  level  width  X  length  of  the  picture  in  pixels 


1 

2 

3 

4 

5 

6 


7x8 
15x16 
30x32 
60X64 
120  X  248 
240  X  256 


The  advantage  of  having  Pyramid  levels  is  the  ability  to  analyze  objects  in 
different  resolution  levels  and  at  diffent  speeds.  For  instance,  if  we  want  to  stop  a  fast 
moving  droplet  from  reaching  the  top  electrode ,  we  need  to  know  its  position, 
velocity,  and  acceleration  in  a  short  time.  In  this  case,  identifying  the  droplet  with 
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good  resolution  is  unnecessary.  Thus,  by  considering  the  whole  frame  transferred  by 

the  camera  as  a  7  X  8  pixel  picture  the  position  of  the  droplet  can  be  calculated  in  a 
maximum  of  56  instructions.  To  perform  the  same  calculation  in  Pyramid  level  6  a 
maximum  of  61 ,440  instructions  have  to  be  performed.  Now,  consider  observing  the 
droplet  at  the  moment  of  sibling  formation  with  the  objective  been  observing  the 
number  of  droplets  and  their  size  distribution.  It  is  obvious  that  in  this  case,  to  identify 
the  correct  number  of  siblings  and  to  identify  the  size  of  these  droplets  with  a  high 
degree  of  accuracy,  we  need  to  capture  only  one  high  resolution  frame.  Thus,  for  this 
application  the  number  of  frames  the  tracker  can  process  per  second  is  of  no  concern. 
In  this  case,  there  are  two  advantages  of  using  Pyramid  level  6  over  Pyramid  level  1. 
The  first  is,  by  dividing  the  frame  transfered  by  the  camera  into  a  grid  consisting  of 
240  x  256  pixels  the  size  of  the  individual  droplets  can  be  measured  with  high 
precesion.  The  second  is,  by  dividing  the  picture  into  many  pixels,  even  the  very  small 
droplets  can  be  identified. 

Note  that  it  is  posible  to  stabalize  the  droplet  using  Pyramid  level  1,  and  then 
once  the  droplet  is  stabalized,  the  droplet  and  its  immediate  surrounding  area  can  be 
scaned  using  Pyramid  level  6.  Thus,  by  using  the  Pyramid  tracker  the  need  to  have  two 
frame  grabbers  and  supporting  equipment-  A  fast  frame  grabber  to  suspend  die 
droplet,  and  a  high  resolution  frame  grabber  to  observe  sibling  formation  —  can  be 
avoided. 

There  is  another  advantage  of  using  the  Pyramid  Tracker  instead  of  a 
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Note  that  it  is  posible  to  stabalize  the  droplet  using  Pyramid  level  1 ,  and  then 
once  the  droplet  is  stabalized,  the  droplet  and  its  immediate  surrounding  area  can  be 
scaned  using  Pyramid  level  6.  Thus,  by  using  the  Pyramid  tracker  die  need  to  have  two 
frame  grabbers  and  supporting  equipment-  A  fast  frame  grabber  to  suspend  the 
droplet,  and  a  high  resolution  frame  grabber  to  observe  sibling  formation  -  can  be 
avoided. 

There  is  another  advantage  of  using  the  Pyramid  Tracker  instead  of  a 
conventional  frame  grabber  for  this  experiment.  When  I  was  attempting  to  suspend  a 
droplet  using  a  30  frames  per  second  frame  grabber3,  though  the  frame  grabber  was 
transfering  30  frames  per  second  to  the  IBM-XT  processor,  the  XT  was  able  to  process 
only  6  frames  per  second.  Thus,  not  counting  the  time  it  would  take  to  calculate  the 
necessary  voltage  difference  between  the  two  plates,  the  XT  would  take  l/5s  to  respond 
to  a  frame  transfered  by  the  camera.  As  will  be  explained  in  the  next  section,  by 
taking  advantage  of  the  special  features  of  the  tracker,  the  COMPAQ  386  could 
respond  to  a  frame  transferred  by  the  tracker  in  l/30s:  a  six  fold  increase  in  the 
response  time  of  the  feedback  loop. 

62  DETECTING  THE  CENTER  OF  THE  DROPLET 

The  following  algorithm  was  used  to  calculate  the  center  of  die  droplet  in  die 
two  experiments  explained  in  die  last  section.  This  same  algorithm  can  be  usd  to 
calculate  the  center  of  die  droplet  in  the  actual  experiment 
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The  basic  algorithm  used  in  calculating  the  center  of  the  droplet  is  as  follows. 

VI,  V2  and  V3  are  integer  variables  with  zero  initial  value. 

1)  Paste  bright  white  paper  to  the  interior  walls  of  the  room  temperature 
chamber.  Now  illuminate  the  inside  of  the  chamber.  By  constantly  checking  the 
highest  grey  leveled  pixel  in  the  frame,  ( the  tracker  should  be  operated  in 
Pyramid  level  6  to  maximize  the  resolution.  The  gray  level  of  a  pixel  can  be 
directly  obtained  from  the  tracker  by  giving  the  address  of  that  pixel.)  provide 
light  to  the  chamber  until  the  highest  grey  level  found  in  the  frame  is  a 
minimum.  Note  this  threshold  level  Tl.  At  this  time  fix  the  position  of  the 
illuminators.  The  light  intensity  level  inside  the  chamber  must  not  be  varied 
hereafter. 

2)  Now  select  a  threshold  level  for  the  droplet.  This  threshold  level  is  such  that  it 
is  slightly  higher  than  Tl .  Thus,  if  we  now  ask  the  tracker  to  find  pixels 
whose  grey  level  is  above  the  threshold  value  the  tracker  should  not  be  able  to 
find  any  such  pixels.  Program  the  tracker  to  await,  expecting  a  dropet. 

3)  Since  the  water  in  the  reservoir  has  been  colored  black,  if  a  droplet  of  this  liquid 
appears  in  the  view  of  the  camera  now,  the  tracker  should  detect  some  pixels 
above  the  threshold  value,  program  the  tracker  to  check  each  pixel  of  the  frame 
starting  from  the  upper  left  hand  and  all  the  way  upto  the  upper  right  hand 
comer.  At  the  end  of  this  row  the  tracker  will  start  checking  the  pixels  starting 
from  the  left  hand  corener  of  the  second  row.  This  process  continues  untill  the 
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tracker  has  checked  all  die  pixels.  All  the  pixels  in  the  first  row  have  x  address 
1  and  all  the  pixels  in  the  second  row  have  x  address  2  etc.  Thus,  in  Pyramid 
level  6,  the  x  address  of  any  pixel  in  the  last  row  is  240.  Similarly  the  y  address 
of  any  pixel  in  the  first  column  is  1 

3)  If  a  pixel  is  encountered  above  the  threshold  value  increment  VI.  Add 
the  x  address  of  this  droplet  to  V2  and  y  address  to  V3.  Repeat  this  procedure 
for  all  the  pixels  in  the  frame. 

4)  At  die  end  of  the  frame  dividide  the  value  in  V2  by  Viand  V3  by  VI.  This 
gives  the  approximate  x  and  y  cordinates  of  the  center  of  the  droplet. 

There  is  actually  no  need  to  scan  the  whole  frame  to  look  for  the  position  of  the 
droplet.  Once  we  determine  the  x  and  y  cordinates  of  the  center  of  the  droplet  from 
the  first  frame,  to  track  the  droplet  thereafter,  though  we  need  to  scan  all  the 
horizontal  lines,  only  few  of  the  vertical  lines  need  to  be  scaned.  As  the  droplet 
stabalizes  we  need  not  scan  all  the  horizontal  lines  either.  In  this  manner  the  COPMAQ 
has  more  time  to  spend  for  the  feedback  loop. 

Another  feature  was  incorporated  into  this  algorothm.  While  scanning  the 
image  if  we  encounter  high  grey  level  pixels  that  implyies  that  the  tracker  is  currendy 
scanning  the  droplet.  In  such  a  case,  if  we  started  scanning  from  the  top,  and  if  we  do 
not  encounter  any  high  value  pixel  for  an  entire  horizontal  line  that  means  that  we  have 
completely  scaneed  the  dropletr.  Tn  this  case  we  can  stop  scanning  the  droplet  and 
start  calculating  the  necesssary  control  voltages  to  the  power  supplies  etc.  Hence  by 
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using  this  method  the  COMPAQ  has  more  time  to  spend  on  calculating  the  feedback 
outputs. 


63  The  Soft  Ware  Listings: 

Listed  below  are  the  general  software  routines  which  die  tracker  requires,  and 
the  special  purpose  routines  which  are  required  to  detect  the  droplet.  The  programs  in 
pages  37-44  were  written  by  me.  The  rest  of  the  programs  were  originally  written  by 
Dr.  Van  der  Wal  of  David  Samoff  Research  Center  and  modified  by  me  to  best  suit  the 
purposes  of  this  experiment. 

Modules  cursor.c  and  reg.c  contains  propriety  iitformation  of  RCA 
Coporation  and  should  be  treated  as  confidential 


XIIULAND.  JUNCTION  OF  .IHEJROG&  AM  PA OE 

Initializes  pyramid  microcode  memory 

after  power  on  37 

makefile  -  makefile  for  main.c  and  gaus.h  38 


t 

i 


34 


gaus.h  -  Gaussian  pyramid  header  module  39 

main.c  -  Generate  gausian  pyramid.  Calculate 
the  center  of  the  droplet.  Determine  the  required 
voltage  difference  between  the  two  electrodes 

to  stabalize  the  droplet.  41 

makefile  -  makefile  for  ips  directory  45 

bus.c  -  Image  processing  system  drivers  46 

cursor.c  -  mouse  functions  48 

gfx.c  -  Graphical  utilities  54 

ips.c  -  image  processing  system  drivers  59 

ips.h  -  macros  for  standard  ips  routines  62 

key.c  -  For  input/output  directed  through  the  keyboard  65 

reg.c  -  Functions  for  register  reading  and  writing  68 
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screen.c  -  Input/output  functions  and  functions  to 
display  the  screen 


state.c  -  Functions  to  manipulate  programs  in  ips 


wr.c  -  writes  code  from  'file.q'  to  control  board 
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Script  Gaus  setup  microcode  to  run  tracker 
#  Initializes  pyramid  microcode  memory  after  power  on 

cpy  -ab  <  gp4.q  !  cpy  -dm  0x2000  0x1000  1 
gaus 


SRC=main.c 

}BJ=main.o 

CFLAGS=-i  -Mse2  -0 

NAME=gaus 

EXE=$(NAME) 

ARCH = .  ./ips/ips .  a 
LIB= 

$(EXE):  $(0BJ)  $(ARCH) 

cc  $(CFLAGS)  -o  $(EXE)  $(0BJ)  $(ARCH)  $(LIB) 


lint: 


lint  -a  -x  $(SRC)  >  1 


WC" 

wc  $(NAME)-h  $(SRC)  makefile 


purge: 


/bin/rm  $(0BJ)  $(EXE) 


-38' 


/*  file:  gaus.h 

desc:  gaussian  pyramid  header  module 
date:  11/12/87  rose 
mod:  11/12/87  rose 

*/ 

/*  include  file  */ 

*if  def  M_XENIX 
^include  ”../ips/ips.h" 

#else 

^include  ./ips/ips .  h" 

#endif 

/*  NULL  fix  */ 

Kifdef  NULL 
#undef  NULL 
#endif 

^define  NULL  0 


/*  macro  function  redefinitions  */ 


#ifdef 

MJffiNIX 

^define 

IP_ADDR(offset) 

^define 

IP_PGM(  s  taddr) 

*#endif 

#ifndef 

M_XENIX 

^define 

IP_ADDR(offset) 

#define 

IP_PGM(staddr) 

#endif 

(ips_base[(offset)  >> 
ioc  tl(ipd ,  'p' ,  s  taddr) 


(ips_base  +  (offset)) 
ips_pgm(s  taddr) 


8]  !  ((offset)  &  Oxffff)) 


/*  macro  constants  */ 


^define  TRK_X  16  /*  width  of  tracker  window  */ 

^define  TRK_Y  15  /*  height  of  tracker  window  */ 


/*  structure  for  pyramid  constants  */ 


struct  pyrlevel  { 
short  uix.uiy, 
unx.uny; 


>; 


<  struct  thlock  { 

short  xDean,ymean; 
short  numpix; 
short  pbesum; 
short  x,y; 

struct  pyrlevel  *pp; 


/*  (xl,xh),(yl,yh)  offset  registers  */ 

/*  image  size  in  (x,y)  directions  */ 


/*  sum  of  x,  sum  of  y  */ 

/*  number  of  pixels  >  thr  */ 

/*  sum  of  pixels  >  thr  */ 

/*  x  and  y  offsets  within  level  */ 
/*  pointer  to  pyramid  level  */ 


/*  static  global  variables  */ 


struct  pyrlevel  Pyr[6]; 


/*  pyramid  constant  values 


*/ 


I 
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/*  external  functions  */ 


/*  main.c  */ 
int 


roainO; 


/*  file:  main.c 

desc:  generate  gaussian  pyramid 
func:  main(argc,argv) 

date:  11/01/87  rase 
rood:  11/06/87  rose 

*/ 

^include  <stdio.h> 

.^include  "gaus.h" 


/* 

static  globed 

variables 

*/ 

struct  pyrlevel  Pyra[] 

=  { 

{ 

7, 

3, 

256, 

240 

>, 

/* 

:goi 

{ 

277, 

0, 

128, 

120 

>, 

/* 

:gi: 

{ 

279, 

121, 

64, 

60 

>, 

/* 

!G2! 

{ 

281, 

182, 

32, 

30 

}, 

/* 

!G3! 

{ 

283, 

213, 

16, 

15 

}. 

/*  ! 

G4! 

{ 

285, 

278, 

8, 

8 

} 

/* 

:g5 

}; 

/*  func:  main(argc,argv) 

desc:  controller  function 
args:  <none> 

ret:  (int)  0  ==  successfull 

.  ==  error 

*/ 

main(argc,argv) 
int  arge; 
char  *argv[]; 


register  short  pix, 
thr, 
xx; 

register  pointer  bp; 
short  numpix  =  0, 
pixsura  =  0, 
xoff  =  0, 
yoff  =  0, 
deep  =  0, 
pixoff  =  0 . 
xroean, 
ymean, 

yy. 

xcenter, 

ycenter, 

dxcenter, 

dycenter, 

i, 

j. 

k; 

byte  stl[64], 
st2[64]; 


struct  thlock  blockB, 
♦B; 

struct  pyrlevel  *pp; 


/*  initial  values 

*/ 

for  (i  =  0;  i  <  5; 

i++) 

/* 

pyramid  loc's 

Pyr[i]  =  Pyra[i]; 
if  (!ips_init()) 

/* 

image  processing  system 

exit(l); 

/*  arguments 

*/ 

for  (i  =  1;  i  <  arge;  i++)  { 

if  (argv[i][0]  ==  ' 

-')  ( 

for  (j  =  1;  j  < 

strlen(argv[i]);  j++) 

switch  (argv[i][j])  { 

case  T: 

break; 

/* 

camera  1  (default) 

case  '2': 

camera  =  0x20; 

/* 

camera  1  (default) 

for  (k  =  0;  k  <  4; 

k++) 

apcval[k]  ,‘=  camera; 

♦ape  =  apcval[0]; 
break; 

default: 

fprintf(stderr,“invalid  flag 

-%c\n",argv[i][j]); 

}; 

}  else 

fprintf(stderr, "invalid  argument  '%s'\n",argv[i]); 

>; 


/*  input  lookup  table 

input  LUT  1,2:  linear 

input  LUT  3:  inverse  linear  ♦/ 

for  (i  =  1;  i  <  4;  i++)  { 

♦ape  =  apcval[0];  *slut  =  slutval[i]; 
for  (j  =  0;  j  <  256;  j++)  { 

♦alut  =  (byte)  j; 

♦dlut  =  (i  !=  3)  ?  j  :  255  -  j; 

}; 


/*  output  lookup  tables  ♦/ 

for  (i  =  0;  i  <  3;  i++)  { 

♦ape  =  apcvalfi  +  1];  ♦slut  =  Oxaa; 
for  (j  =  0;  j  <  256;  j++)  { 

♦alut  =  (byte)  j; 

♦dlut  =  j; 

}; 


/♦  control  registers 

set:  camera->pyr->(v0a,v0b  &  v2a)  ♦/ 

♦(v0  +  Ctrl)  =  ♦(vl  +  Ctrl)  =  ♦(v2  +  Ctrl)  =  0x47;  /♦  256x240 

♦ab  =  0x44;  *cd  =  0x34;  ^e  =  2;  *out  =  0x07;  /♦  paths 


£alu  =  0x03; 

*(v0  +  ctrh)  =  0x0c; 

♦slut  =  0x55; 

/♦  initialize  state  variable  program  ♦/ 

/♦  To  show  GO  -  use  v2.ctrh=0x08  when  sst==0  ♦/ 

/♦  To  show  Gn  -  use  v2.ctrh=0x08  when  sst==l  */ 

if  (!st_encode(stl,"sst!=l  sst==l  v0.pan=0xl0Q  ")) 

exit(l); 

if  (!st_encode(st2,"sst!=0  sst=  =  0  v0.pan=0x000  v2.ctrh=0x08")) 
exit(l); 

/♦  initialize  blocks  ♦/ 

B  =  &(blockB); 

B->pp  =  &(Pyr[4)); 

B->x  =  B->y  =  0; 
for  (i=0;i  <  18;i++) 
printfC  \n"); 

printfC  GAUSSIAN  PYRAMID  TRACKER  11/12/8  7\n"); 

printfC  \n“); 

printfC  \n“); 

printfCPress  Ctrl/Backspace  to  stop  tracker\n“); 
printfC  \n"): 

printf ("Initializing  tracker.  Keep  off  from  the  view  of  the  camera. \n“); 
printfC  \n"); 

/♦  write  .q  microcode  program  to  buffer 

#ifdef  N0T_N0W 
if  (!wrips("gp4.q")) 
exit(l); 

#endif  ♦/ 

pp  =  B->pp;  /*  this  and  next  line  were  taken  from  loop  ♦/ 
dO  =  dO  +  DELY  *  (long)(pp->uiy  +  B->y)  +  (pp->uix  +  B->x); 

/*  loop  forever  and  determine  the  center  of  the  droplet  ♦/ 

fb  =  d2; 
while  (TRUE)  { 

IP_PGM(stl); 

/*  do  calculations  here  ...  ♦/ 

IP_PGM(st2); 
bp  =  dO; 

xroean  =  ymean  =  numpix  =  pixsum  =  0; 
for  (yy  =  TRK_Y;  yy;  yy--)  { 

for  (xx  =  16;  xx;  xx- -) 

if  ((pix  =  *(bp++))  >  20)  { 
numpix++; 
pixsum  +=  pix; 
xroean  +=  xx; 
ymean  +=  yy; 
printf ("Xd\n" ,  pix); 

}; 

bp  +=  DELY  -  TRK_X; 

}; 


/*  alu  A+B  ♦/ 

/♦  cont  grab  ♦/ 
/*  linear  lut  ♦/ 
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TRACKER  READY\n"); 


/* 


if  (deep  ==  700)  { 
deep  =  701; 
printfC 

} 

if  (  deep  ==  701)  { 

numpix  =  nurapix  -  pixoff; 
if  (numpix  >  0  )  { 

ycenter  =  (ymean  -  yoff)/  numpix; 
ycenter  =  16  -  ycenter; 
dycenter  =  240*ycenter/15; 
xcenter  =  (xroean  -  xoff)/  numpix; 
xcenter  =  17  -  xcenter; 
dxcenter  =  256*xcenter/16; 

} 

> 

if  (deep  <  700)  { 
xoff  =  xroean; 
yoff  =  ymean; 
pixoff  =  numpix; 
deep  =  deep  +1; 

} 

drawline(dxcenter  -3,  dycenter  .dxcenter  +3, dycenter, 255); 

drawline(dxcenter  .dycenter  +3,  dxcenter  .dycenter  -  3,255); 


ips_close(); 
return(0);  */ 


Sfc=bus.c  cursor.c  gfx.c  ips.c  key.c  reg.c  screen.c  state.c  wr.c 
OBJ=bus.o  cursor.o  gfx.o  ips.o  key.o  reg.o  screen.o  state. o  wr.o 
CFLAGS=-i  -Mse2  -0 
IAME=ips 
ARCH=ips.a 
LIB= 

$(ARCH):  $(0BJ)  $(0BJ) 

ar  ru  $(ARCH)  $(0BJ) 
ranlib  $(ARCH) 

lint: 

lint  $(SRC)  >  1 


wc  $(NAME).h  $(SRC)  makefile 


purge: 


/bin/rm  $(0BJ)  $(ARCH) 


/*  file:  bus.c 

desc:  image  processing  system  drivers 
func:  bus_init() 

bus_close() 
date:  07/02/87  msc 

mod:  10/29/87  rose 

*/ 

^include  <stdio.h> 

^include  <fcntl.h> 

^include  "ips.h" 

#ifndef  M.JCENIX 
^include  <sys/hwdbus.h> 

^include  <sys/ioctl.h> 

#endif 

/*  func:  bus__init() 

desc:  initialise  hardware  bus 
args:  <none> 

ret:  bool:  FALSE  =  error 

TRUE  =  successful! 

*/ 

bool 

bus_init() 

{ 

byte  b; 
pointer  bp; 

4ifndef  M_XENIX 
struct  hwdbase  pb; 

#endif 

/*  open  image  processing  system  as  file  descriptor  */ 

tfifndef  M_XENDC 

/*  open  bus  device  */ 

if  ((ipd  =  open("/dev/bus",0_RDWR))  ==  ERR)  { 
fprintf(stderr, "cannot  open  bus\n”); 
retum(FALSE); 

>; 

/*  get  base  address  */ 

if  ((ioctKipd, HWDBASE, &pb))  ==  ERR)  { 
fprlntf(stderr,"ioctl  error\n“); 
retum(FALSE); 

ips_base  =  (byte  *)  pb.hwdmero  +  OxOOaOOOOOL; 

/*  wake  up  hardware  */ 

b  =  *(pb.hwdio  +  0x0200); 

•else 

/*  open  bus  device  */ 

if  ((ipd  =  open("/dev/ip",0_RDWR))  ==  ERR)  { 

fprintf(stderr, "cannot  open  device  '/dev/ip'\n");[ 


l 


-  V-6- 


retum(FALSE); 


r 


}; 

/*  get  initial  segment  descriptor  -  'inquiry'  */ 

if  (ioctl(ipd,'i',&ips_base))  { 

fprintf(stderr, "cannot  access  '/dev/ip'Xn"); 
retum(FALSE); 

}; 

#endif 

/*  read  and  write  a  test  byte  */ 

bp  =  ips_base  +  0x10;  b  =  *bp; 

*bp  =  0x00; 
if  (*bp  !=  0x00)  { 

fprintf(stderr, "could  not  access  image  processing  system\n“); 
retum(FALSE); 

}; 

*bp  =  Oxff ; 
if  (*bp  !=  Oxff)  { 

fprintf(stderr, "could  not  access  image  processing  system\n"); 
retum(FALSE); 

>; 

*bp  =  b; 
retum(TRUE); 

} 

/*  func:  bus_close() 

desc:  terminate  bus  interface  parameters 
args:  <none> 

ret:  <void> 

*/ 

void 

bus_close() 

{ 

if  (close(ipd)  ==  ERR) 

fprintf(stderr,"error  closing  bus\n"); 

> 


I 
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/*  file:  cursor.c 

desc:  mouse  functions 
func:  cursor_jLnit() 

cursor_close() 
rawxy(buttons  ,mx,roy) 
isbutO 

getxy(mx,ray  ,f  unc) 
drawcursor(mx,my,func) 
date:  07/07/87  msc 
mod:  11/29/87  msc 

*/ 

♦♦include  <stdio.h> 

♦♦include  <fcntl.h> 

#ifdef  M_XENIX 
“♦include  <termio.h> 

♦♦include  <sys/ioctl.h> 

♦♦else 

♦♦include  <sys/hft.h> 

#endif 

♦♦include  "ips.h" 


/*  macro  constants 

*/ 

♦♦define 

MAX_CURSOR_X 

8 

/* 

mouse  cursor  size 

*/ 

♦♦define 

MAX_CURSOR_Y 

14 

“♦define 

CURSOR_ERASE 

FALSE 

/* 

flags  sent  to  cursor_draw() 

V 

♦♦define 

CURSOR_DRAW 

TRUE 

#ifdef 

MJCENIX 

♦♦define 

MOUSE_DEVICE 

"/dev/ttyla" 

/* 

mouse  device  i/o  name 

*/ 

♦♦define 

MAX_RETRIES 

500 

/* 

max  number  of  error  retries 

*/ 

“♦endif 

/*  static  variables  */ 

#ifndef  MJCENK 

static  char  mousectrlt]  =  {  ESC.TYx, NULL, NULL, NULL, OxOe.HFKSRPROH, 

HFKSRPROL, 0x04, NULL, NULL,NULL,NULL,  0x80, NOLL, 0x80, 
ESC, TYx', NULL, NULL, NULL, OxOc.HFLOTHCH, 

HFLOTHCL, 0x02, NULL, NULL, 0x02, NULL, 0x02  >; 

#endif 

/*  func=  cursor_init() 

desc:  initialize  mouse  cursor 
ergs:  <none> 

ret:  bool:  FALSE  ==  unsuccessful! 

TRUE  -=  successfull 

*/ 

,bool 

cursor_init() 

{ 

#ifdef  HJXENXX 
int  i; 
byte  b  =  ' 
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static  struct  terraio  ttmouse; 


/*  for  communication  line  */ 


if  ((mfd  =  open(MOUSE_DEVICE,0_RDWR  {  0_NDELAY»  ==  ERR)  { 
fprintf(stderr, “could  not  open  '%s'\n",MOUSE_DEVICE); 
retum(FALSE); 

}; 


ttmouse. c_iflag  =  0;  ttmouse. c_oflag  =  0;  ttroouse.c_lflag  =  0; 
ttmouse.c_cflag  =  B1200  !  CS7  !  CSTOPB  !  CREAD  J  PAREHB  !  CLOCAL; 
if  (ioctl(mfd,TCSETA,&ttmouse)  ==  ERR)  { 
fprintf(stderr, "could  not  ioctl  info\n“); 
retum(FALSE); 

>; 


for  (i  =  0;  i  <  MAX_RETRIES;  i++)  { 
if  (wrrte(mfd,&b,l)  =-  -1)  { 

fprintf(stderr, "could  not  write  to  mouse\n“); 
retuim(FALSE); 

>; 

if  (rdchk(mfd» 
break; 

>; 

if  (i  ==  MAX_RETRIES)  { 

fprintf(stderr, "mouse  is  not  attached\n"); 
retum(FALSE); 

>; 

Seise 

wrate(tfd,mousectrl,32); 

Sendif 


return(TRUE); 

} 


/* 


%/ 


func: 

desc: 

args: 

ret: 


cursor„close() 

close  cursor  routine 

<none> 

<void> 


void 

cursor_close() 

f 

Sifdef  M_XENIX 
close(mfd); 

Seise 

mousectrl[16]  =  NULL; 

write(tfd,raousectrl,17); 

close(tfd); 

Sendif 

} 


/*  func: 
desc: 
args: 

ret: 


*/ 


rawxy(buttons , rax, my) 
returns  raw  mouse  data 

short  *buttons,  *  button  status  -  0x01  and  0x02  * 

*mx,*my;  *  relative  mouse  values  * 

(bool):  FALSE  ==  error 

TRUE  ==  successful! 


I 
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bool 

rawxy(buttons  ,mx,my) 
short  ♦buttons, 

*mx, 

♦my; 

r 

short  r; 

static  short  blast; 

#ifdef  M_XENIX 
static  char  ch[3], 
t; 

static  short  n; 

/♦  clear  old  values  */ 

♦buttons  =  blast;  *mx  =  ♦my  =  0; 

/♦  wait  for  new  value(s)  ♦/ 

while  ((r  =  read(mfd,ch  +  n,l»  ==  1)  { 
switch  (n)  { 

case  0:  if  (ch[0]  < 
continue; 

t  =  (ch[0]  &  OxOf); 

if  ((t  !=  0x00)  &&  (t  !=  0x03)  &&  (t  !=  0x0c)  &&  (t  !=  OxOf)) 
continue; 
break; 

case  1:  if  (ch[l]  >=  '@')  { 
ch[0]  =  ch[n]; 
n  =  0; 
continue; 

}; 

if  ((!(t  &  0x03)  '  !(ch[l]  &  0x20))  =  =  1)  { 
n  =  0; 
continue; 

}; 

break; 

case  2:  if  (ch[2]  >=  '@')  { 
ch[0]  =  ch[n]; 
n  =  0; 
continue; 

}; 

if  ((!(t  &  0x0c)  ~  !(ch[2]  &  0x20))  ==  1)  { 
n  =  0; 
continue; 

}; 

}; 

if  (++n  <  3) 
continue; 
n  =  0; 

♦buttons  =  blast  =  (ch[0]  >>  4)  &  0x03; 

♦wx  =  ch[l]  -  ((ch[l]  &  0x20)  <<  1); 

♦my  =  ch[2]  -  ((ch[2]  &  0x20)  «  1); 
return(TRUE); 

}; 

if  (r  ==  -1)  { 


♦ 


-  so  - 


fprintf(stderr," trouble  reading  mouse"): 
retum(FALSE); 

#else 

/♦  check  input  queue  */ 
iskey(); 

/*  clear  old  values  */ 

if  (mqbegin  ==  roqend)  { 

♦buttons  =  blast;  *mx  =  ♦my  =  0; 

/*  wait  for  new  value(s)  */ 

}  else  { 

blast  =  *buttons  =  mqbegin- >buttons; 

♦mx  =  mqbegin- >dx; 

♦my  =  -mqbegin  ->dy; 
if  (++mqbegin  >=  mqueue  +  MAX_QUEUE) 
mqbegin  =  mqueue; 

}; 

#endif 

retum(TRUE); 

> 

/*  func:  isbutO 

desc:  returns  button  value 
args:  <none> 

ret:  (bool):  button  up  (FALSE)  or  down  (TRUE) 

*/ 

bool 

isbutO 

r 

short  buttons,mx,my; 

if  (rawxy(&buttons,&rax,&my)) 
retum(buttons  !=  NULL); 
return(FALSE); 

■v 

/*  func:  getxy(mx,my,func) 

desc:  gets  x  &  y  coordinates 
draws  cursor  on  screen 

returns  when  left  or  right  button  is  pressed 
args:  (int)  *mx,*my;  ♦  x  &  y  values  ♦ 

(int)  func;  *  type  of  function  ♦ 

ret:  short:  button  value 

♦/ 

short 

,  qetxy(xx,y  y  ,f  unc) 
short  *xx, 

*yy; 

short  func; 

{ 

short  buttons, 

,  but, 


-S*- 


height, 

width, 

dx, 

dy, 
oldx, 
oldy; 

static  short  mx, 
my; 


/*  flush  cursor  input  stream  */ 

if  (func  ==  C(JRSOR_FLUSH)  { 
if  (!rawxy(&buttons,&dx,&dy)) 
retum(NULL); 
do  { 

but  =  buttons;  height  =  dx;  width  =  dy; 
if  (!rawxy(&buttons  ,&dx,&dy)) 
return(NOLL); 

>  while  ((buttons  !=  but)  !!  (dx  !=  height)  |!  (dy  !=  width)); 
/*  cursor  */ 

}  else  { 

drawcursor(oldx  =  mx,oldy  =  my,CURSOR_DRAW); 
do  { 

/*  get  raw  data  -  update  */ 

if  (!rawxy(&buttons,&dx,&dy)) 
retum(NULL); 
oldx  =  mx;  oldy  =  my; 
mx  +=  dx;  my  +=  dy; 
if  (mx  <  0)  mx  =  0; 
if  (my  <  0)  my  =  0; 
if  (mx  >=  xraax)  mx  =  xmax  -  1; 
if  (my  >=  ymax)  my  =  ymax  -  1; 

/*  redraw  cursor  */ 

if  ((oldx  !=  mx)  !!  (oldy  !"  my))  { 

drawcursor(oldx,oldy,CURSOR_ERASE); 

drawcursor(mx,roy,CURSOR_DRAW); 

}; 

/*  return  condition  */ 

}  while  (((func  ==  CURSOR_PRESS)  &&  (Ibuttons))  ;• 

((func  ==  CURSOR_RELEASE)  &&  (buttons))); 
drawcuraor(mx,my,CURSOR_ERASE); 

}; 

/*  return  values  */ 


*xx  =  mx; 

*yy  =  my; 

retum(buttona ); 


/*  func:  drawcursor(mx, ray, func) 

-desc:  draw  or  erase  cursor 
args:  (short)  rax, 

my; 

(short)  func;  0=DRAW,  1=ERASE 
ret:  <void> 

*/ 

void 

drawcursor(mx,my,func) 
short  mx, 
my, 
func; 

{ 

short  i, 

height; 

pointer  addr  =  fb  +  ADDR(mx,my); 
byte  *cptr, 
bit; 

static  byte  bit_data(MAX_CURSOR_Y  *  MAX_CURSOR_X], 

cursor_data[]  =  {  0x01, 0x03, 0x07,0x0f  ,0xlf  ,0x3f  ,0x7f , 

Oxf f , Oxlf ,  Oxlb, 0x19 , 0x31, 0x3 0 , 0x2 0  }; 

cptr  =  cursor_data;  i  =  0; 

for  (height  =  MAX_CURSOR_Y ;  height--;  cptr++)  { 
for  (bit  =  1;  bit;  bit  <<=  1)  { 
if  (*cptr  &  bit)  { 

if  (func  ==  CURSOR_DRAW)  { 
bit_data[i++]  =  *addr; 

*(addr++)  =  0x50  !  color_cursor; 

}  else 

*(addr++)  =  bit_data[i++]; 

}  else 
addr++; 

>; 

addr  +=  DELY  -  MAX_CURSOR_X 

}; 
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/*  file:  gfx.c 

desc:  graphical  utilities 
func:  drawline(xO,xl,yO,yl,color) 

dra  wf  rame(x  0  ,xl ,  y  0 ,  y  1 ,  c  olo  r) 
drawbox(xO,xl,yO,yl,color) 
r() 

date:  03/09/87  rjk 
mod:  12/10/87  msc 

*/ 

^include  "ips.h" 

/*  func:  drawline(xO,yO,xl,yl,color) 

desc:  draw  a  line  on  an  image/frame 

args:  (short)  x0,y0;  *  starting  point  * 

(short)  xl,yl;  *  ending  point  * 

(byte)  color;  *  value  to  use,  or  -1  for  exclusive-oring  * 

ret:  <void> 


yoid 

irawline(xO  ,y  0  ,xl,yl, color) 
short  xO, 

yo, 

xl, 

yi; 

byte  color; 

[ 


short  1, 

expr, 

el, 

e2, 

inc, 

ddx, 

ddy, 
len; 

pointer  addr; 

/*  vertical  line  */ 


if  (xO  ==  xl)  { 

if  ((xO  <  0)  !!  (xO  >=  xroax)) 
return; 
if  (yl  <  y0)  { 

1  =  y0; 
y0  =  yl; 
yi  =  l; 

}; 

if  ((y0  >=  ymax)  ,'!  (yl  <  0)) 

return; 

if  (y0  <  0)  yO  =  0; 

if  (yl  >-  ymax)  yl  =  ymax  -  1; 

addr  =  fb  +  ADDR(x0,y0); 

for  (1  =  yl  -  yO;  1--  >=  0;  addr  +=  DELY) 
if  (color  ==  Oxff) 

♦addr  '=  Oxff; 
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else 

*addr  =  (*addr  &  color_mask)  !  color; 

/*  horizontal  line  */ 

>  else  if  (yO  ==  yl)  { 

if  ((yO  <  0)  i!  (y 0  >=  xmax)) 
return; 
if  (xl  <  xO)  { 

1  =  xO; 
xO  =  xl; 
xl  =  1; 

>; 

if  ((xO  >-  xmax)  !!  (xl  <  0)) 
return; 

if  (xO  <  0)  xO  =  0; 
if  (xl  >=  xmax)  xl  =  xmax  -  1; 
addr  =  fb  +  ADDR(xO,yO); 
for  (1  =  xl  -  xO;  1--;  addr++) 
if  (color  ==  Oxff) 

♦addr  ~=  Oxff; 
else 

♦addr  =  (♦addr  &  color_mask)  |  color; 

/♦  bresenham  line  algorithm  */ 

)  else  { 

if  ((xO  <  0)  iJ  (xO  >=  xmax)  \\  (xl  <  0)  'j  (xl  >  = 
(yO  <  0)  !J  (yO  >=  ymax)  |!  (yl  <  0)  !i  (yl  > 
return; 

ddx  =  ABS(xl  -  xO); 
ddy  =  ABS(yl  -  yO); 
if  (ddy  <=  ddx)  { 
if  (xl  <  xO)  { 

1  -  xO; 
xO  =  xl; 
xl  =  1; 

1  =  yO; 
yO  =  yl; 
yl  =  1; 

>; 

inc  =  (yl  >  yO)  ?  DELY  :  -DELY; 
el  =  (expr  =  (e2  =  (ddy  <<  1))  -  ddx)  -  ddx; 
len  =  ddx  +  1; 
addr  =  fb  +  ADDR(x0,y0); 
while  (len — )  { 
if  (color  ==  Oxff) 

♦addr  ~=  Oxff; 
else 

♦addr  =  (♦addr  &  color_mask)  i  color; 
if  (expr  >  0)  { 

addr  +=  inc  +  1; 
expr  +=  el; 

}  else  { 

addr++;  expr  +=  e2; 

}; 

}; 

}  else  { 

if  (yl  <  y  0)  { 
l  =  yO; 


xmax)  !! 

=  ymax)) 
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yO  =  yl; 
yl  =  1; 

1  =  xO; 
xO  =  xl; 
xl  =  1; 

}; 

Inc  =  (xO  >  xl)  ?  -1  :  1; 
len  =  ddy  +  1; 

el  =  (expr  =  (e2  =  (ddx  <<  1))  -  ddy)  -  ddy; 
addr  =  fb  +  ADDR(xO,yO); 
while  (len--)  { 

if  (color  ==  Oxff) 

*addr  ~=  Oxff; 
else 

♦addr  =  (♦addr  &  color_mask)  !  color; 
if  (expr  >  0)  { 

addr  +=  DELY  +  inc; 
expr  +=  el; 

}  else  { 

addr  +=  DELY; 
expr  +=  e2; 

}; 

>; 

}; 

}; 

> 

'*  func:  drawfrarae(xO,yO, xl.yl, color) 

desc:  draw  a  frame  on  an  image/frame 
args:  (int)  xO.yO;  ♦  one  corner  ♦ 

(int)  xl.yl;  ♦  another  corner  ♦ 

(byte)  color;  ♦  value  to  use,  or  -1  for  exclusive-oring  ♦ 

ret:  <void> 

♦/ 

void 

dra  wf  rame(x0 ,  y  0  ,xl ,  y  1 ,  colo  r) 
int  xO, 

y  0, 
xl, 
yl; 

byte  color; 

{ 

drawline(xO  ,y0,xl,y0,color); 
drawline(xl,y  0  ,xl,yl,  color); 
drawline(xl,yl,xO  ,yl, color); 
drawline(xO  ,yl,x0  ,y  0 , color); 

> 

/♦  func:  drawbox(xO,yO,xl,yl,color) 

desc:  draw  a  box  on  an  image/frame 
args:  (Int)  xO.yO;  ♦  one  comer  ♦ 

(int)  xl.yl;  ♦  another  comer  ♦ 

(byte)  color;  ♦  value  to  use,  or  -1  for  exclusive-oring  ♦ 

ret:  <void> 

♦/ 
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void 

drawbox(xO,yO,xl,yl, color) 
int  xO, 
yo, 
xl, 
yl; 

byte  color; 

{ 

short  1, 

len; 

pointer  addr; 

/*  left  to  right  */ 

if  (xl  <  xO)  { 

1  =  xO; 
xO  =  xl; 
xl  =  1; 

>; 

if  (yl  <  yO)  { 

1  =  yO; 
yO  =  yl; 
yl  =  1; 

>; 

/*  loop  by  rows  */ 

xl  -=  (xO  -  1); 

for  (addr  =  fb  +  ADDR(xO.yO);  y0++  <=  yl;  addr  +=  DELY  -  xl)  { 
len  =  xl; 

if  (color  ==  Oxff) 
while  (len--) 

*(addr+"0  Oxff; 

else 

while  (len — )  { 

*addr  =  (*addr  &  color_raask)  J  color; 
addr++; 

}; 

>; 

> 

/*  func-'  clear(v) 

desc:  clears  frame  buffer,  redisplays  menu 

args:  pointer  v;  *  frame  to  grab  * 

ret:  <void> 

*/ 

void 
clear(v) 
pointer  v; 

{ 

short  i, 

len, 

start; 
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/*  clear  lower  three  bits  of  vO  */ 


*alu  =  0x00; 

*(v  +  mask)  =  0xf8; 

*(v  +  ctrh)  =  0x08; 
len  =  (*(v  +  Ctrl)  ==  0x40)  ?  11  :  7; 
start  =  *(v  +  pan)  +  512  *  *(v  +  pan  +  1); 
*(v  +  pan)  =  (512  +  start  -  len)  &  Oxff; 

*(v  +  pan  +  1)  =  (512  +  start  -  len)  >>  8; 
for  (i  =  0;  i  <  2;  i++)  { 

while  (!(*(v  +  ctrh)  &  0x40)); 
while  (*(v  +  ctrh)  &  0x40); 

}; 

*(v  +  mask)  =  0x07; 

*(v  +  pan)  =  start  &  Oxff; 

*(v  +  pan  +  1)  =  start  »  8; 


/*  file:  ips.c 

desc:  image  processing  system  drivers 
func:  ips_init() 

ips_close() 
date:  07/02/87  rose 

mod:  12/10/87  msc 


# include  "ips.h" 


/*  func: 
desc: 
args: 
ret: 


*/ 


ips_init() 

initialize  ips  hardware  parameters 
<none> 

bool:  FALSE  =  error 

TRUE  =  successful! 


bool 

ips_init() 

{ 

register  short  i, 

j; 

short  k; 
pointer  v[3]; 

static  byte  slutidx[]  =  {  0x00,0x55,0xaa.Gxff  }; 


/*  initialize  bus  */ 


if  (!bus_init()) 
return(FALSE); 

/*  initialize  registers  */ 

if  (!reg_init()) 
retura(FAI.SE); 


/*  default  values  */ 

xroax  =  MAXJL.OFFSET;  yraax  =  MAX_Y_OFFSET; 
color_jmask  =  0x00;  color_cursor  =  Oxff; 
fb  =  dO; 
camera  =  0x00; 


/*  screen  size  */ 
/*  color  scheme  */ 

/*  frame  address  */ 
/*  camera  input  */ 


/*  initialize  'apeval'  and  'slutval'  arrays  */ 


for  (i  =  0;  i  <  4;  i++)  { 

apcvalfi]  =  0x40  *  i; 
slutval[i]  =  slutidxfij; 

}; 


/* 

initial  values  of  registers  */ 

/* 

input  lookup  table 

*/ 

/* 

input  LUT  0:  zero 

*/ 

/* 

input  LUT  1,2:  linear 

*/ 

/* 

input  LUT  3:  inverse  linear 

*/ 

S9 


*apc  =  0x00; 

♦slut  =  slutval[i); 
for  (j  =  0;  3  <  256; 

♦alut  =  (byte)  j; 
*dlut  =  (i  ==  0)  ? 

}; 

>; 


j++)  ( 

0  :  ((i  !=  3)  ?  j  :  255 


/*  output  lookup  tables  */ 
/*  output  LOT  0:  zero  */ 
/*  output  LUT  1,2,3:  linear  */ 


for  (i  =  1;  i  <  4;  i++)  { 

♦ape  =  apcval[i]; 
for  (j  =  0;  j  <  256;  j++)  { 

♦alut  =  (byte)  j; 
for  (k  =  0;  k  <  3;  k++)  { 

♦slut  =  (byte)  slutvalfk]; 
♦dlut  :  (k  ==  0)  ?  0  :  j; 

>; 

}; 


/♦  alu  &  mux  registers  */ 

♦ape  =  0x00; 

♦slut  =  0x55; 

♦kl  =  *k 2  =  »k3  =  0; 

♦alu  =  ♦shft  =  ♦mult  =  0; 

♦offs  =  0x00; 

♦gain  =  0;v80; 

♦ab  =  0x44; 

♦cd  =  0x44; 

♦e  =  0x00; 

♦out  =  0; 

/*  reset  x,y, pan, control, mask  registers 


v[0]  =  v0; 
v[l]  =  vl; 
v[23  =  v2; 

for  (i  =  0;  i  <  3;  i++)  { 

/♦  set  for  full  screen  (512x480)  */ 

*(v[i]  +  Ctrl)  =  0x40; 

♦(v[i]  ♦  mask)  =  0x00; 

/♦  initialise  frame  registers  */ 

for  0=0;  d  <  2;  J++)  { 

*(v[i]  +  x  +  j)  =  0x00; 

♦(v(i]  ♦  y  +  j)  =  0x00; 

*(v[i]  +  pan  +  j)  =  0x00; 

♦(v[i3  +  scroll  +  5)  -  0x00; 


}; 

•v  /*  wait  one  frame  time  */ 

while(!(*(vO  +  ctrh)  &  0x40)); 
while(*(v0  +  ctrh)  &  0x40); 

/*  clear  screen  */ 

*(v[i]  +  ctrh)  =  0x08; 

}; 

/*  wait  two  frame  times  */ 

for  (i  =  0;  i  <  2;  i++)  { 

while(!(*(v0  +  ctrh)  &  0x40)): 
whlle(*(v0  +  ctrh)  &  0x40); 

>; 

return(TROE); 

} 

/*  func:  ips_close() 

desc:  close  ip  system 
args:  <none> 

ret:  <void> 

*/ 

void 

\ps_close() 

{ 

reg_close(); 

bus_close(); 

} 


I 
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/*  file:  ips.h 

desc:  macros  for  standard  ips  routines 
date:  07/01/87  msc 
mod:  11A0/87  msc 

*/ 

/*  type  definitions  */ 

typedef  unsigned  char  byte; 
typedef  char  bool; 

#ifdef  M_XENIX 

typedef  byte  huge  *  pointer; 

((else 

typedef  byte  *  pointer; 

#endif 


/*  macro  boolean  macros  */ 


#undef 

NULL 

((define 

NULL 

0 

((define 

FALSE 

0 

((define 

TRUE 

1 

((define 

ERR 

-1 

/*  macro  register  constants 

*/ 

((define 

REG 

1 

/* 

register  (kl,k2,ab,etc..) 

*/ 

((define 

FRAME 

2 

/* 

frame  (v0,vl,v2) 

((define 

OFFSET 

3 

/* 

offset  (x,y,pan,scroll,etc..) 

*/ 

/*  macro  constant  maximums  */ 


((define 

MAX_X_OFFSET 

512 

/* 

max  x  pixel  offset 

*/ 

((define 

MAX_Y_OFFSET 

480 

/* 

max  y  pixel  offset 

*/ 

((define 

MAX-LABEL 

16 

/* 

maximum  length  of  register  name 

*/ 

((define 

MAX-REGS 

32 

/* 

maximum  number  of  system  registers  */ 

((define 

MAX-QUEUE 

512 

/* 

max  queue  size 

*/ 

/*  macro  keyboard  constants 

V 

((define 

BACK 

' 

/* 

backspace 

*/ 

((define 

BRK 

0x7f 

/* 

break  key 

*/ 

{(define 

DEL 

/* 

delete  key 

*/ 

{(define 

ESC 

Oxlb 

/* 

escape  key 

*/ 

{(define 

RET 

An' 

/* 

return  key 

*/ 

/*  macro  cursor  constants 

*/ 

{(define 

CURSOR_FLUSH 

0 

/* 

flush  cursor  input 

*/ 

{(define 

CURSOR_GET 

1 

/* 

return  cursor  parameters  stat 

*/ 

{(define 

CURSOR_PRESS 

2 

/* 

wait  until  button  is  pressed 

*/ 

((define 

CURSOR_RELEASE  3 

/* 

wait  until  button  is  released 

*/ 

/*  macro  image  constants 

V 

{(define 

DELY  ! 

512 

/* 

default  value  for  image.dely 

*/ 

I 
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/*  macro  functions  */ 


^define 

MIN(x.y) 

((x  <=  y)  ?  x  : 

y) 

^define 

MAX(x,y) 

((x  >=  y)  ?  x  : 

y) 

#define 

ABS(x) 

(((x)  >=  0)  ?  (x)  :  -(x» 

^define 

ABS2(x,y) 

(((x)  >=  y)  ?  (x 

1 

ts 

W 

/-s 

>» 

1 

#define 

SIGN(x) 

((x  <  0)  ?  -1  :  1) 

^define 

ADDR(x,y) 

((DELY  *  (long) 

(y))  +  (x)) 

/*  image  processing 

data  structures  */ 

3truct  regtype  { 

char  label[MAX_LABEL];  /* 

label  of  register 

pointer  addr; 

/* 

register  address 

short 

type, 

/* 

register  type  (REG 

offset, 

/* 

register  offset 

•  f 

len; 

/* 

register  length  (1  ■ 

#ifndef 

M_XENIX 

struct  qtype  { 
short  dx.dy, 

buttons ; 

}; 

#endif 


/*  mouse  delta  x  & 
/*  mouse  buttons 


2) 


y  values 


*/ 

*/ 

*/ 

*/ 

*/ 


*/ 

*/ 


/*  global  variables  */ 

pointer  ips_base;  /* 

int  ipd;  /* 

#ifdef  M_XENIX 

int  mfd;  /* 

#else 

int  tfd;  /* 


struct  qtype  mqueue[MAX_QUEUE],  /* 


*roq begin, 

♦roqend; 

#endif 

char  cqueue[MAX_QOEUE],  /* 

♦cqbegin, 

♦cqend; 

pointer  d0,dl,d2,  /* 

v0,vl,v2,  /* 

pyrmero.sst,  /* 

dlut.alut, slut, ape,  /* 

kl,k2,k3,ab,cd,e,  /* 

alu.shft, out, mult,  /* 

offs, gain;  /* 

short  x,y, pan, scroll,  /* 

Ctrl,  ctrh,  mask;  /* 

pointer  fb;  /* 

byte  color_mask,  /* 

color_cursor;  /* 

i  short  xmax.ymax;  /* 

short  char_width,char_height;  /* 

byte  apcval[4],  /* 

slutval{4],  /* 

camera;  /* 

short  num_reg;  /* 

I 


base  address  for  image  proc.  system  */ 


image  proc.  system  file  descriptor  */ 

mouse  file  descriptor  */ 

terminal  file  descriptor  */ 

queue  for  mouse  cursor  stats  */ 


queue  for  character  input  */ 


frame  buffer  data  addresses  */ 

frame  buffer  control  addresses  */ 

pyramid  program, status  register  */ 

lookup  tables  */ 

multipliers ,  channels  */ 

alu  parms  */ 

camera  offset  &  gain  controls  */ 

frame  buffer  x  &  y  parms  */ 

frame  buffer  control  parms  */ 

current  frame  buffer  */ 

number  of  bits  for  color  ~  Oxff  */ 

color  of  scren  cursor  */ 

maximum  screen  size  for  gfx, cursor  */ 

size  of  character  font  */ 

index  for  ape:  0x00, 0x40, 0x80, OxcO  */ 

index  for  slut:  0x0 0,0x5 5, Oxaa, Oxff  */ 
camera  0x00  or  0x20  -  'or'ed  to  ape  */ 
number  of  registers  */ 
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struct  regtype  regs[MAX_REGS];  /* 
/*  system  functions  */ 


•shar  *strcpy(),*stmcpy().  /* 

*strcat(),*stmcat(),  /* 


*getenv(); 

A 

long  lseekO; 

/* 

int  atoi(); 

/* 

long  atol(); 

/* 

iouble  atof(); 

/* 

char  *malloc(); 

/* 

void  free(); 

/* 

/*  external  functions  */ 

y*  bus.c  */ 

bool 

bus_Jnit(); 

void 

bus_close(); 

/*  cursor.c 

*/ 

bool 

cursorJnitO; 

void 

cursor_close(); 

bool 

rawxy(); 

bool 

isbut(); 

short 

getxy(); 

void 

drawcursorO; 

/*  gfx.c  */ 

void 

drawlineO; 

void 

drawframeO; 

void 

drawboxO; 

void 

clear(); 

ips.c  */ 

bool 

ips_jLnit(); 

void 

ips_close(); 

/*  keyboard.c 

*/ 

bool 

key_init(); 

void 

key_close(); 

char 

getch(); 

bool 

iskey(); 

/*  main.c  */ 

int 

main(); 

/*  reg.c  */ 

bool 

reg_lnit(); 

void 

reg_dose(); 

/*  screen.c 

*/ 

bool 

screen_init()-, 

void 

screen_close(); 

void 

drawstr(); 

void 

drawch(); 

void 

®rr(); 

/*  wr.c  */ 

bool 

wripsQ; 

data  structure  for  register  info  */ 


string  library  functions  */ 

string  library  functions  */ 

return  home  environment  */ 

i/o  functions  */ 

string  to  int  conversion  function  */ 

string  to  long  conversion  function  */ 
string  to  float  conversion  function  */ 
dynamic  memory  allocation  routine  */ 

return  allocated  memory  to  system  */ 


- 
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/*'  file:  key.c 

desc:  keyboard  i/o 
func:  key_AnitO 

key_close() 

getchO 

iskeyO 

date:  07/14/87  rase 
mod:  11/11/87  rase 

*/ 

/*  include  files  */ 


^include  <fcntl.h> 
#ifndef  M_XENIX 
^include  <stdio.h> 
^include  <sys/hft.h> 
#endif 

^include  "ips.h" 


/*  func:  key_init() 

desc:  initialize  keyboard 
args:  <none> 

ret:  bool:  FALSE  ==  unsuccessful! 

TRUE  ==  successfull 

*/ 

bool 

’:ey_init() 

{ 

/*  initialize  */ 


#ifdef  MJIENIX  ...  *  ..v 

systeraCstty  -echo  -brkint  -isig  -istnp  -icanon  eof  ), 

Mpl  c  p 

char  kstr[80]; 


sprintf  (ks  tr ,  ”/de  v/hf  t/%d "  ,ioc  tl(  0 ,  HFG  CH  A  N ,  0 )) ; 

if  ((tfd  =  open(kstr,0_RDWR  !  0_NDELAY))  ==  ERR) 

if  ((tfd  =  open(‘7dev/tty",0_RDWR  !  0_NDELAY))  ==  ERR)  { 
fprintf(stderr ."cannot  open  vterm\n"); 
exit(l); 

}; 


systemC’stty  -enhedit  ignpar 
roq begin  =  mqend  =  requeue; 
ttendif 

cq  begin  =  cqend  =  cqueue; 


-echo  -brkint  -isig  -istrip  -icanon  eof  "); 
/*  reset  cursor  queue  */ 

/*  reset  character  queue  */ 


retum(TRUE); 

} 


£ 
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/*  func: 
desc: 
args: 
ret: 


*/ 


key_close() 
close  keyboard 
<none> 

<void> 
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void 

key_close() 

{ 

''ifndef  MJCENIX 
close(tfd); 

#endif 

system("stty  sane  echoe  ixany  tab3“); 

> 


'*  func:  getchO 

desc:  receive  character  through  keyboard 
args:  <none> 

ret:  char:  character  typed  in,  or  NULL  if  none  in  keyboard  buffer 

*/ 

'har 

getchO 

{ 


char  ch  =  NULL; 

if  (iskeyO)  { 
ch  =  *cqbegin; 

if  (++cqbegin  >=  cqueue  +  MAX_QUEUE) 
cqbegin  =  cqueue; 

>; 


retum(ch); 

} 

/*  func:  iskey() 

desc:  checks  keyboard  for  character  typed 
args:  <none> 

ret:  bool:  FALSE  ==  no  character  entered 

TRUE  ==  character  pressed 

*/ 

bool 

iskeyO 

f 


char  ch  =  NULL; 

kifdef  M_XENIX 

while  (rdchk(0)  >  0)  { 
read(0,&ch,l); 

♦cqend  =  ch; 

if  (++oqend  >=  cqueue  +  MAX_QUEUE) 
cqend  =  cqueue; 

}; 

Jtelse 
short  i; 
struct  nice  { 
short  deltax, 
deltay; 
char  sec[3]; 
char  sixty ths , 
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buttons, 

stype; 

}  mouse; 

while  (read(tfd,&ch,l)  ==  1) 
if  (ch  ==  ESC)  { 

if  ((readftfd.&ch.l)  ==  1)  &&  (ch  ==  '[')) 
if  ((read(tfd,&ch,l)  ==  1)  &&  (ch  ==  y'))  { 
for  (i  =  0;  i  <  sizeof(struct  mice);  i++) 
if  (read(tfd,((char  *)  &raouse)  +  i,l)  <  1) 
break; 

mqend->dx  =  mouse.deltax; 
mqend->dy  =  mouse.deltay; 
reqend->buttons  =  mouse. buttons; 
if  (++roqend  >=  requeue  +  MAX_QUEUE) 
reqend  =  requeue; 

}; 

ch  =  NOLL; 

}  else  { 

♦cqend  =  ch; 

if  (++cqend  ==  cqueue  +  MAX_QUEUE) 
cqend  =  cqueue; 

}; 

#endif 

retum(cqbegin  !=  cqend); 

} 


* 
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/*  file:  reg.c 

desc:  functions  for  register  reading/writing 
func:  reg_init() 

reg_close() 
date:  05/14/87  rose 
mod:  11/14/87  rose 

V 

♦♦include  <stdio.h> 

-♦include  "ips.h" 

/*  static  variables  */ 

static  struct  fldsettype  { 
char  *label; 
pointer  *addr; 
short  ^offset; 

}  fldsetf]  =  { 

{"v0",&v0, NULL}, f'd0",&d0, NULL},  {"vl",&vl, NULL}, {“dl",&dl, NULL}, 

{"v2",&v2, NULL},  {"d2",&d2, NULL},  {"pyrroem",&pyrroera, NULL},  {"sst",&sst, NULL}, 
{"dlut",&dlut, NULL}, {“alut",&alut, NULL}, {"slut", &slut, NULL}, {“ape”  ,&apc, NULL}, 
{"kl",&kl, NULL}, {“k2",&k2, NULL}, {"k3",&k3, NULL}, {"ab",&ab, NULL}, 

{“cd",&cd, NULL}, {"e",&e,NULL},{"alu",&alu, NULL}, {”shft",&shft, NULL}, 

{"out"  ,&out,  NULL},  {"mult"  ,&mult,  NULL},  {"of  fs",&offs,  NULL},  {"gain",  &gain,  NULL}, 
{"x”,NULL,&x},{"y",  NULL, &y},  {"pan",  NULL,  &pan},  {"scroll",  NULL, &scroll}, 

{"Ctrl"  .NULL  ,&ctrl} ,  {"ctrh"  .NULL  ,&ctrh} ,  {"mask" ,  NULL  ,&mask} , 

{"".NULL, NULL} 


/*  func:  reg_jnit() 

desc:  initialize  registers 
args:  <none> 

ret:  (bool)  FALSE  ==  error 

TRUE  ==  successful! 

*/ 

bool 

reg_init() 


char  *hp, 

str[80], 
type[40], 
saddr[40], 
attr[80], 
env[40]; 
short  tp; 
long  raddr, 

Ian; 

FILE  *fd; 

struct  fldsettype  *fs; 
struct  regtype  *rp; 

/*  clear  variables  */ 

num_reg  =  0;  /*  number  of  registers  */ 
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/*  get  home  directory  */ 

if  0(hp  =  getenv("HOME")))  { 

fprintf(stderr, "could  not  get  home  directory\n"); 
retum(FALSE); 

}: 

/*  initialize  registers  */ 

/*  get  hardware  stats  from  file  '$HOME/.cliphw'  */ 

sprint£(str,"fcs/.cliphw”,hp); 

if  (!(fd  =  fopen(str,"r“))) 

fd  =  fopen("/usr/clipAoroe/cliphw","r‘'); 
if  (fd)  { 

while  Ofeof(fd)) 

if  (fgets(s tr ,7  9  ,f d))  { 
if  (!strncmp(str,"/*",2)) 
continue; 

if  (sscanf(str,"%sS6s%s%s",type,env,saddr,attr)  ==  4)  { 
if  (!strncrop(saddr>"0x",2)) 

sscanf(saddr  +  2,"%lx",&raddr); 
else 

raddr  =  atol(saddr); 
if  (!stmcmp(attrI"0x",2)) 

sscanf(attr  +  2,“%lx"1&len); 
else 

len  =  atoKattr); 

/*  find  specific  variable  */ 

for  (fs  =  fldset;  (fs->addr)  j!  (fs->offset);  fs++) 
if  (!strcmp(fs->label,env)) 
brcsk* 

if  (.(fs->addr)  &&  !(fs->offset))  { 

fprintf(stderr, "unable  to  assign  register  '%s'\n",env); 
retum(FALSE); 

>; 

/*  assign  value  */ 

if  (!strcmp(type)"frame"))  { 
tp  =  FRAME; 

*((fs  +  l)->addr)  =  ips_base  +  len; 

*(fs->addr)  =  ips_base  +  raddr; 

}  else  if  (!strcmp(type, "offset"))  { 
tp  =  OFFSET; 

*(fs->offset)  =  raddr; 

}  else  if  (!strcmp(type,"reg")  ii  !strcmp(typel"addr"))  { 
tp  s  REG; 

*(£s->addr)  =  ips_base  +  raddr; 

}  else  { 

fprintf(stderr, "invalid  type  '%s'  for  label  'Xs'\n",type,env); 
return(FALSE); 

}; 

/*  add  to  register  table  */ 

rp  =  regs  +  num_reg++; 
strcpy(rp~>label,env); 
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rp->type  =  tp; 
rp->offset  =  raddr; 
rp->len  =  len; 

}  else 

fprintf(stderr,“sscanf  failed  on:  '%s*\n",str); 

}; 

>  else  { 

fprintf(stderr, “could  not  open  $HOME/.cliphw  £ile\n“); 
return(FALSE); 

>; 


retum(TROE); 

} 

/*  func:  reg  closeO 

desc:  close  registers 
args:  <none> 
ret:  <void> 

*/ 

void 

reg_close() 

> 


i 
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/*'  file: 

screen.c 

desc: 

functions  for 

screen  display 

func: 

screen_init() 

screen_clo  s  e() 

drawstr(sp,row,col,color) 

dra  wch(ch,  row  ,col,  color) 

date: 

06/26/87  rose 

mod: 

*/ 

11/16/87  msc 

# include 

<stdio.h> 

((include 

<fcntl.h> 

((include 

"ips.h" 

/*  macro  constants 

*/ 

((define 

F0NT_HDR  LEN 

0x2c 

((define 

DDDFSIZE 

0x00 

((define 

FNTCOLMN 

Oxlc 

((define 

FNTROWS 

Oxle 

((define 

FNTLKD 

0x28 

‘define 

FNTBUG 

0x60 

/*  static  variables  */ 

static  struct  fonts  {  /*  information  concerning  font 

byte  *hptr, 

*dptr, 

*tptr; 

}  font; 


/*  func: 
desc: 
args: 
ret: 


*/ 


screen_init(fontname) 

Initialize  screen  and  load  font 

char  *fontname;  *  name  of  font  to  use  * 

(bool)  FALSE  ==  error 

TRUE  ==  successful! 


bool 

>creen_Anit(fontname) 
char  *fontname; 

{ 


int  fd; 

char  str[80]; 

unsigned  int  len; 

/*  load  font  */ 

if  ((fd  =  open(fontname,0_RD0NLY))  ==  ERR)  { 
sprintf(str,"/usr/cl±p/fontAs",fontname); 
fd  =  open(str,0_RD0HLY); 

>; 

if  (f d  ==  ERR)  { 

fprintf(stderr, "cannot  open  font  file\n"); 
return(NULL); 

>; 

l  if  (lseektfd.dong)  FHTBUG.O)  ==  ERR)  { 


fprintf(stderr, "cannot  lseek  font  file\n“); 

close(fd); 

return(NULL); 

}; 

If  (font.hptr)  { 

free(font.hptr); 
font.hptr  =  NULL; 
free(font.dptr); 
font.dptr  =  NULL; 

}; 

font.hptr  =  (byte  *)  roalloc(FONT_HDR_LEN); 
if  (read(fd,font.hptr,(int)  FONT_HDR_LEN)  <  FONT_HDR_LEN)  { 
fprintf(stderr, "cannot  read  font  header\n"); 
f  ree(f  ont  .hpt  r) ; 
font.hptr  =  NULL; 
dose(fd); 
retum(NULL); 

>; 

/*  len'  note  -  assumes  <64K  files  */ 

len  =  256  *  *(font.hptr  +  DDDFSIZE  +  2)  +  *(font.hptr  +  DDDFSIZE  +  3); 
char_width  =  256  *  *(font.hptr  +  FNTCOLMN)  +  *(font.hptr  +  FNTCOLMN  +  1); 
char_height  =  256  *  *(font.hptr  +  FNTROWS)  +  *(font.hptr  +  FNTROWS  +  1); 
font.dptr  =  (byte  *)  roalloc(len); 
if  (read(fd,font.dptr,len)  ==  ERR)  { 

fprintf(stderr, "cannot  read  font  data\n"); 

free(font.hptr); 

font.hptr  =  NULL; 

close(fd); 

return(NULL); 

}; 

if  (close(fd)  ==  ERR)  { 

fprintf(stderr, "cannot  close  font  file\n"); 
retum(NULL); 

}; 

/*  lku'  note  -  assumes  <64K  files  */ 

len  =  256  *  *(font.hptr  +  FNTLKU  +  2)  +  *(font.hptr  +  FNTLKU  +  3); 
font.tptr  =  font.dptr  +  len  -  FONT_HDR_LEN; 

return(TRUE); 

> 

/*  func:  screen_close() 

desc:  close  screen,  free  memory 
args:  <none> 
ret:  <void> 

*/ 

/old 

screen_cloae() 

{ 

free(font.hptr); 

free(font.dptr); 


/*  func:  draws tr(sp,row,col,color) 
desc:  draws  string  on  screen 
args:  (char  *)  sp;  *  string  to  display  * 


» 
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(int)  x,y;  *  location  * 

(byte)  color;  *  color  * 

-  ret:  <void> 

*/ 

void 

urawstr(sp,xO,yO,  color) 
char  ♦sp; 
int  xO, 
y  0; 

byte  color; 

{ 


for  (;  *sp;  xO  +=  char_width) 
drawch(*(sp+ +),x0 ,  y  0  .color); 

> 

/*  func:  drawch(sp, row, col, color) 

desc:  draws  character  on  screen 
args:  (char)  ch;  *  string  to  display  * 

(int)  x.y;  *  location  * 

(byte)  color;  *  color  * 
ret:  <void> 

*/ 

void 

drawch(ch,xO  ,y  0  .color) 
char  ch; 
int  xO, 
y  0; 

jyte  color: 

{ 


pointer  addr; 
short  i, 

yi, 

width, 
height; 
byte  *tptr, 

♦cptr, 

bit; 

/*  calculate  table  location,  width,  height  */ 

tptr  =  font.tptr  +  4  *  ch; 
width  =  *(tptr  +  1)  &  0x3f; 
height  =  char_height  -  (*tptr  »  3)  - 

((♦tptr  &  0x07)  <<  2)  -  (*(tptr  +  1)  >>  6); 
cptr  =  font.dptr  +  256  ♦  *(tptr  +  2)  +  *(tptr  +  3) 

/♦  center  character  ♦/ 

xO  +=  (char_width  -  width)  /  2; 
yl  =  yO  +  (♦tptr  »  3); 

/♦  loop  ♦/ 

addr  =  fb  +  ADDR(xO,yl); 

for  (bit  =  128;  height — ;  addr  +=  DELY  -  width) 
for  (1=0;  i  <  width;  i++,  addr++)  { 

•  if  (*cptr  &  bit)  { 
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if  (color  ==  Oxff) 

♦addr  “=  Oxff; 
else 

♦addr  =  (*addr  &  color_mask)  J  color; 

}; 

if  ((bit  »=  1)  ==  0)  { 
cptr++; 
bit  =  128; 


/*  file: 
desc: 
func: 


date-* 

mod: 


*/ 


state.c 

functions  for  ips  program  manipulation 
sklpch(sp) 
st_enccde(st,sp) 
ips_pgm(st)  *  A IX  only  * 

05/08/87  rase 
11/17/87  rase 


•include  <stdio.h> 
^include  <ctype.h> 
^include  "ips.h" 


/*  global  variables  */ 
itatic  char  *valid[]  =  {  M=+-*/X&!“  <>•#=!=".  _ 

•=+-*/%&r<>@#*'!:'‘ 


/*  func:  skipch(sp) 

desc:  skip  over  space,  tabs,  and  carriage  returns 

args:  (char  *)  sp;  *  string  *  L 

ret:  (char  *):  pointer  to  NULL  or  non-white  space  character 

*/ 

char  * 
skipch(sp) 

:har  *sp; 

{ 


while  (*sP  &&  ((*sp  ==  '  1  !!  (*sP  ==  At')  !!  (*sp  ==  '\n'))) 
sp++; 

return(sp); 


/*  func: 
desc: 
args: 

ret: 


*/ 


s  t_encode(s  t ,  s  trptr) 
encode  state  information  from  string 
(byte  *)  st;  *  state  to  encode  to  * 

(str  *)  str;  *  string  with  info  * 
variable:  FALSE  ==  error 

TRUE  ==  successful! 


bool 

st_encode(bp,sp) 
byte  *bp; 
char  *sp; 

{ 


short  offset, 
len, 
i, 

yul  * 

char  label[403, 
opl, 
op2, 
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*lp; 

byte  rg; 

/*  decode  */ 
while  (*sp)  { 

/*  clear  */ 

*bp  =  Oxff ; 

/*  optional  semicolon  */ 

if  (*sp  &&  (*sp  == 
sp  =  skipch(sp  +  1); 

/*  register  */ 

sp  =  skipchfsp); 

for  (lp  =  label;  *sp  &&  (isalnura(*sp));  *(lp++)  =  *(sp++)); 

*lp  =  NULL; 
sp  =  skipch(sp); 

for  (i  =  0;  i  <  num_reg;  i++) 
if  (!strcmp(label,regs[i].label)) 
break; 

if  (i  ==  num_reg)  { 

fprintf(stderr, "register  '%s'  not  defined\n“, label); 
return(FALSE); 

>; 

offset  =  regs[i]. offset; 
len  =  regs(i].len; 
if  (regs[i].t ype  ==  FRAME)  { 
if  (*sp  ==  '.')  { 

sp  =  skipch(sp  +  1); 

for  (lp  =  label;  *sp  &&  (isalnum(*sp));  *(lp++)  =  *(sp++)); 
*lp  =  NULL; 
sp  =  skipch(sp); 

for  (j  =  0;  j  <  num_reg;  j++) 
if  (!strcmp(label,regs[j].label)) 
brGdk* 

if  ((j  ==  num_reg)  !!  (regs[j].type  !=  OFFSET))  { 

fprlntf(stderr,"'%s'  is  not  an  offset  variableNn”  .label); 
return(FALSE); 

>; 

offset  +=  regs[j]. offset; 
len  =  regst  j].len; 

}  else  { 

fprlntf(stderrl"'.'  required  after  frame  variable\n"); 
return(FALSE); 

>; 

>  else  if  (regs[i].type  !=  PEG)  { 

f printf ( a  tderr , s '  is  not  a  register\n"  .label); 
return(FALSE); 

>; 

rg  =  ((offset  &  Oxlff)  ==  0x100)  ?  0x80  :  offset  &  Oxff; 

*(bp  +  1)  =  len; 

/*  operation  */ 


I 
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opl  =  *(sp++); 

op2  =  (*sp  ==  '=')  ?  *(sp++)  :  ' 
for  (i  =  0;  valid[0][i];  i++) 

if  ((opl  ==  valid[0][i])  &&  (op2  ==  valid[l][i])) 
break; 

if  (!valid[0][i])  { 

fprintf(stderr, "illegal  operation  *%c'\n",*sp); 
return(FALSE); 

*(bp  +  2)  =  valid[2][i]; 
sp  =  skipch(sp); 

/*  value  */ 

for  (lp  =  label;  *sp  &&  (isalnum(*sp));  *(lp++)  =  *(sp++)); 

*lp  =  NULL; 
sp  =  skipch(sp); 

/*  value  or  mask  */ 
if  (!strncmp(label,"0x",2)) 

sscanf (label  +  2,"%hx“,&val); 
else  if  (isdigit(*label)) 
val  =  atoi(label); 
else  if  (*(bp  +  2)  ==  '=')  { 

/*  register  (from  copy  register)  */ 

*(bp  +  2)  = 
val  =  -1; 

for  (i  =  0;  i  <  num_reg;  i++) 
if  (!strcmp(label,regs[i].label))  { 
val  =  regs[i]. offset; 
if  (regs[i].type  ==  FRAME)  { 
if  (*sP  ==  V)  { 

sp  =  skipch(sp  +  1); 

for  (lp  =  label;  *sp  &&  (isalnum(*sp));  *(lp++)  =  *(sp++)); 
*lp  =  NULL; 
sp  =  skipch(sp); 

for  (j  =  0;  j  <  num_reg;  j++) 
if  (!  s trcmp(label , reg s [  j] .label))  { 
vail  .+=  regs[j], offset; 
break; 

}; 

if  (j  ==  num_reg) 
val  =  -1; 

}  else  { 

fprintf(stderr/'Y  required  after  frame  variable\n”); 
retum(FALSE); 

}; 

>; 

}; 

if  (regs[i].type  !=  REG) 
val  =  -1; 

else 

val  =  ((val  &  Oxlff)  ==  0x100)  ?  0x80  :  val  &  Oxff; 

}; 

if  (val  ==  -1)  { 

f prlntf ( s tderr , "illegal  entry  *Xs'\n" .label); 
retum(FALSE); 

}; 

*(bp  +  3)  =  val  &  Oxff; 
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*(bp  +  4)  =  val  »  8; 

/*  next  */ 

*bp  =  rg; 
bp  +=  5; 

♦bp  =  Oxff ; 

>; 

retum(TRUE); 

} 

wifndef  M_XENIX 
/♦  fun c:  ips_jpgro( table) 

desc:  program  image  processing  system  hardware  registers 
args:  byte  ♦table;  *  state  variable  buffer  address  * 

ret:  <void> 

*/ 

void 

ips_pgro(  table) 
byte  ♦ table; 

{ 


short  num; 

register  pointer  addr; 
while  (♦ table  !=  Oxff)  { 

addr  =  ips_base  +  ((*table  ==  0x80)  ?  0x100  :  *tahle); 

switch  (♦(table  +  2))  { 

case  /♦  set  value  */ 

♦addr  =  *(table  +  3); 
if  (*(table  +  1)  ==  2) 

♦(addr  +  1)  =  ♦(table  +  4); 
break; 

case  /♦  wait  while  equal  to  constant  ♦/ 

num  =  ♦(table  +  3); 
while  (♦addr  ==  num); 
break; 

case  T:  /♦  wait  while  NOT  equal  to  constant  */ 
num  =  ♦(table  +  3); 
while  (♦addr  !=  num); 

}; 

table  +=  5; 

}; 

} 

tfendif 
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/*  file:  wr.c 

desc:  write  code  from  'file.q'  to  control  board 
func:  wrips(wname) 
date:  05/08/87  rose 
rood:  11/01/87  rose 

*/ 

♦♦include  <stdio.h> 

♦♦include  <fcntl.h> 

♦include  <ctype.h> 

♦♦include  "ips.h" 

/*  macro  constants  */ 


♦♦define 

MAX_STRING 

80 

/* 

maximum 

string  length 

*/ 

♦define 

MAX_PYR_MEMORY 

4096 

/* 

maximum 

size  of  ips  program  ram 

*/ 

♦♦define 

MAX_RETRIES 

25 

/* 

maximum 

number  of  write  retries 

*/ 

/*  local  variables  */ 
byte  buffer[2  *  MAX_PYR_MEMORY]; 

/*  func:  wrips(wnaroe) 

desc:  download  file.q  to  pyramid  board  code  space 

-  this  program  adds  code=0  to  the  end  of  the  file 
args:  (char  *)  wnarae:  *  file.q  name,  eg.  ‘Ip3.q"  * 

ref-  (bool)  FALSE  ==  error 

TRUE  ==  successfull 

*/ 

bool 

wrips(wnaroe) 
char  *wnaroe; 


short  1, 

n  =  0, 
pos  -  0; 
byte  *bp; 
char  bO, 
bl, 

str[MAX_STRING]; 
register  byte  *np; 
int  fd; 

/*  open  file  */ 

if  ((fd  =  open(wnaroe,0_RD0NLY))  ==  -1)  { 
sprln  tf ( s  t  r ,  "/u  s  r/clip/qA  s  " .  wn  a  roe ) ; 
if  ((fd  =  open(str,0_RD0NLY))  ==  -1)  { 

fprintf(stderr, “unable  to  open  'Xs*\n",wname); 
retum(FALSE); 

}; 

}; 

/*  read  file  */ 


-  79 


if  ((1  =  read(f d ,buf f e  r ,  2  *  MAX_PYR_MEMORY))  <=  0)  { 
fprtntf(stderr, "error  %d  while  reading  '%s'\n",l,wname); 
retum(FALSE); 

}; 


/*  convert  */ 

for  (bp  =  np  =  buffer;  bp  -  buffer  <  1;  *(np++)  =  16  *  (short)  bO  +  bl)  { 

bO  =  *(bp++);  bl  =  *(bp++); 
if  (!isxdigit(bO)  !i  iisxdigit(bl))  { 

fprintf(stderr,"invalid  char  at  location  %d  in  file  %s\n",l,wnaroe); 
return(FALSE); 

>; 

bO  =  (bO  >  '9')  ?  toupper(b0)  -  'A'  +  10  :  bO  -  'O'; 
bl  =  (bl  >  '9')  ?  toupper(bl)  -  'A'  +  10  :  bl  -  'O'; 


/*  file  with  zeroes  */ 

for  (1  /=  2;  np  <  buffer  +  MAX_PYR_MEMORY;  *(np++)  =  0); 
do  { 

/*  check  for  too  many  times  */ 
if  (n++  >=  MA3CRETRIES)  { 

fprintf(stderr, "could  not  write  to  ips  program  meroory\n“); 
return(NOLL); 

}; 


/*  make  sure  the  pointer  starts  on  a  WORD  boundary  */ 

pos  &=  Oxfffe; 

/*  write  to  pyramid  memory  */ 

#ifdef  MJCENIX 

lseek(ipd,(long)  (pyrmem  -  ips_base)  +  pos,0); 

write(ipd,bp  =  buffer  +  pos,MA)CPYR_MEMORY  -  pos); 

((else 

memcpy(pyrmem  +  pos, bp  =  buffer  +  pos,MAX_PYR_MEMORY  -  pos); 
((endif 

/*  read  from  pyramid  memory  */ 

(dfdef  tUCBMX 

lseek(ipd,(long)  (pyrmem  -  ips_base)  +  pos.O), 

read(ipd,np  =  buffer  +  MAX_PYR_MEMORY  ,1  -  pos); 

((else 

meacpy(np  =  buffer  +  MAX_PYR_MEMORY .pyrmem  >  pos,l  -  pos); 
((endif 

/*  compare  */ 

for  (;  pos  <  1;  pos++) 
if  (*(bp++)  !=  *(np++» 

break; 

>  while  (pos  <  1); 


/*  close  file  */ 


w 


if  (close(fd)  ==  ERR)  { 

fprintf(stderr, "could  not  close  file  '%s'‘,wnaine); 
return(FALSE); 

}; 

/*  else,  OK!  */ 
retum(TRUE); 


I 
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Figure  1:  System  Block  Diagram 


Microprocessor  (PC-XT)  Parallel  Port 


Figure  2:  Microprocessor  to  High  Voltage  Pover 
Supply  Interface  via  D/A  Converters  and  Inverting 
Op.  Amps. 


Figure  3:  A  Schematic  Draving  of  Electrodes 


NOTE: 

(1)  is  the  tope  ted  aluminum  line,  (3)  and  (6)  we  plexiglas  top  and  bottom  vails 
of  da  chamber,  (4)  and  (5)  aie  connections  to  the  aluminum  zincs,  (7)  is  an 
aluminum  zinc ,  (8)  is  the  aluminum  disk ,  (9)  is  a  plexiclas  support,  (10) 
is  da  tubing  vhich  leads  to  vetor  supply,  and  (11)  is  da  injection  nozzle 
made  Dy stracning  out  ua  glass  tube. 
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FIGURE  4  :  Circuit  Diagram  of  the  High  Voltage  Unit 


( taken  from  reference  11) 
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Abstract 


Theoretical  modelling  using  the  maximum  entropy  formalism  indicates  that 
droplet  charging  in  electrostatic  sprays  will  exhibit  phase  transition-like 
behavior.  This  research  is  directed  to  the  collection  of  detailed  electrostatic 
spray  data  and  their  interpretation  in  terms  of  theoretical  models  of  charged 
sprays.  It  involves  the  design  and  construction  of  a  low  charge-to-mass  ratio 
quadrupole  mass  spectrometer  and  an  accompanying  vacuum-sprayer  facility. 
Via  a  high  speed  electrometer,  the  electrostatic  spray  is  analyzed  at  various 
charge-to-mass  ratios  to  determine  its  charge  and  size  distributions.  A 
comparison  of  the  initial  Octoil  data  with  the  maximum  entropy  model  provides 
the  first  direct  evidence  for  the  predicted  phase  transition-like  behavior  in  spray 
charging.  In  addition,  the  previously  published  value  for  the  characteristic 
electrostatic  spray  constant  (-a'/P')  is  confirmed.  The  ability  of  this  facility  to 
provide  detailed  charge  and  size  distributions  can  facilitate  close  examination  of 
other  factors  that  may  affect  the  characteristics  of  the  electrostatic  sprays. 


Chapter  I 

Electrostatic  Atomization 


1.1  Introduction 

Fluid  atomization  is  a  key  process  in  engineering.  Its  applications  include  combustion,  chemical 
deposition,  colloid  production,  space  propulsion,  and  spray  painting  to  mention  a  few.  Various 
devices  are  available  for  spray  production  including  air  atomizers,  swirl  atomizers,  and  impact 
atomizers.  These  devices  use  aerodynamic  shear  force,  centrifugal  force,  and  impact  disruption 
respectively  to  achieve  atomization.  Electrostatic  atomization,  in  which  stresses  are  electrically 
applied  to  cause  atomization,  is  yet  another  means  of  producing  sprays. 

In  these  sprays,  unique  properties  are  observed  as  a  result  of  the  free  electrical  charges 
present  in  the  spray  droplets.  The  spray  cloud  becomes  self -dispersive  due  to  the  mutual 
repulsive  forces  among  charged  droplets  of  similar  polarity.  The  motion  of  the  spray  cloud  is 
limited  by  space  charge  effects  and  the  droplets  have  an  affinity  for  electrically  grounded  objects. 
The  electrostatic  spray  must  be  characterized  by  both  a  charge  distribution  and  a  size  distribution. 

The  primary  goal  of  this  research  is  to  examine  the  two  fundamental  properties  of  the 
electrostatic  sprays,  the  charge  and  size  distributions.  For  this  purpose  a  new  experiment  has 
been  constructed  involving  a  quadrupole  mass  spectrometer  and  a  high  speed  electrometer  that 
permit  both  the  spray  droplet  charge-to-mass  ratio  and  the  charge  to  be  independently  measured. 

1.2  Electrostatic  Atomization  Phenomena 

Electrostatic  atomization  can  be  easily  achieved  with  simple  apparatus.  Figure  1.1  illustrates  the 
operation  of  a  diode-type  sprayer4.  Spray  fluid  is  contained  in  the  reservoir  (or  a  capillary) 
having  an  orifice  at  the  apex  into  which  a  conducting  wire  is  dipped.  This  arrangement  allows 
the  potential  of  the  fluid  to  be  raised  with  respect  to  the  electrical  ground.  At  high  potentials,  the 
fluid  meniscus  at  the  orifice  is  drawn  into  fluid  filaments.  Instabilities  then  develop  along  these 
filaments  causing  them  to  break  up  into  charged  droplets. 
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Figure  1.1  Simple  Diode  Sprayer 


This  phenomena  was  first  described  by  Bose  in  1745  and  more  recently  by  Zeleny5-6, 
who  was  interested  in  the  breakdown  phenomena  of  various  liquid  and  metal  point  conductors 
when  they  were  raised  to  high  potentials.  Zeleny  found  that  the  discharge  current  from  the  liquid 
points  is  carried  mainly  by  charged  liquid  droplets  (the  remaining  cunrent  being  carried  by  ions 
and  electrons).  It  was  an  early  belief  that  the  applied  potential  on  these  liquid  points  governed  die 
atomization  process.  Further  experiments  performed  by  Drozin7  concluded  that  it  was  the 
electric  field  rather  than  the  applied  potential  that  was  causing  the  spraying. 

Various  theoretical  studies  have  been  published  on  electrostatic  atomization.  Most  of 
these  invoke  stability  analysis  of  the  charged  fluid  surfaces,  from  which  stability  criterion  for  the 
surface  and  the  most  probable  charge-to-mass  ratio  of  the  droplet  is  derived  as  a  result  of 
unstable  surface  perturbations**. 10.  The  stability  criterion  derived  for  spherical  liquid  droplets  is 
due  to  Rayleigh  11.  Rayleigh  demonstrated  that  a  limit,  known  to  us  now  as  the  Rayleigh  Limit, 
exists  for  the  charge-to-mass  ratio  of  a  droplet  of  given  size.  At  this  limit  the  stabilizing  surface 
tension  force  is  just  balanced  by  the  destabilizing  electrostatic  force.  This  is  an  important  limit 
for  electrostatic  atomization  becuase  it  restricts  the  size  of  a  droplet  if  its  charge-to-mass  ratio  is 
fixed.  Since  Drozin’s  experiments  showed  that  the  charge-to-mass  ratio  of  the  droplet  formed  is 
proportional  to  the  charging  electric  field,  it  follows  that  the  Rayleigh  Limit  is  a  physical  limit 
which  governs  the  maximum  size  of  a  liquid  droplet  to  be  formed  under  a  specified  set  of  sprayer 
operating  conditions. 

Another  charged  fluid  surface  study  was  performed  by  Taylor12.1^.  He  found  that, 
under  the  influence  of  an  electric  field,  a  charged  liquid  meniscus  deforms  from  a  hemispherical 
geometry  to  that  of  a  cone  (now  known  as  the  Taylor  Cone).  This  cone  sustains  a 
characteristic  semi-vertical  angle  of  49.3  degrees.  The  Taylor  Cone  has  been  experimentally 
verified  for  both  hydrocarbon  fluids  and  liquid  metals.  At  the  tip  of  these  cones,  liquid  is  drawn 
into  filaments  which  subsequently  break  up  into  droplets.  This  mode  of  charged  droplet 
formation  is  to  be  contrasted  to  the  droplet  break  up  when  a  charged  droplet  reaches  die  Rayleigh 
Limit.  In  an  electrostatic  spray,  both  mechanisms  are  believed  to  be  operating  at  the  same  time. 
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The  sprayer  used  in  this  thesis  is  one  that  utilizes  the  Taylor  Cone  for  droplet  formation  and 
dispersion. 

1.3  Models  and  Experimental  Data 

Three  classes  of  models  of  the  electrostatic  sprays  have  been  published  —  the 
electrohydrodynamic,  the  minimum  energy,  and  the  maximum  entropy  models.  The 
electrohydrodynamic  model  is  accurate  in  describing  the  break  up  of  a  single  charged  liquid 
filament  under  the  influence  of  external  perturbation.  However,  this  model  is  not  very  good  for 
describing  electrostatic  sprays  because  in  a  sprayer  like  the  one  shown  in  Figure  1.1,  many  fluid 
filaments  are  usually  formed.  Since  the  actual  number  of  filaments  cannot  be  predicted 
beforehand,  it  is  not  practical  to  use  this  model  to  predict,  for  example,  the  charge-to-mass  ratios 
of  the  droplets  formed. 

The  other  models  were  formulated  to  circumvent  the  difficulties  involved  in  the  modelling 
of  the  detailed  breakup  process  during  atomization.  The  characteristic  of  these  models  is  that 
they  ignore  the  complicated  macroscopic  details  during  breakup,  using  instead  thermodynamic 
conditions  about  the  system  to  derive  the  equilibrium  state.  Therefore,  these  models  are  end-state 
models  in  which  only  the  initial  and  the  final  state  of  the  spray  system  is  of  importance. 

The  minimum  energy  model  was  developed  by  Hendricks17  using  a  concept  proposed  by 
Vonnegut16.  The  model  assumes  that  all  observed  droplets  are  formed  from  larger  droplets 
which  are  at  the  Rayleigh  Limit.  As  a  result,  the  surface  of  these  droplets  becomes  unstable  to 
external  perturbations.  They  eventually  break  into  droplets  6f  smaller  diameter  and  having  a 
different  charge-to-mass  ratio  than  the  parent  droplet.  Assuming  the  final  droplets  are 
monodisperse  (i.e.  equal  in  size)  and  each  having  the  same  charge-to-mass  ratio,  an  expression 
for  the  combined  electrostatic  and  surface  energy  can  be  obtained.  The  combined  energy  may  be 
minimi  wd  with  respect  to  the  radius  by  taking  its  derivative.  This  thermodynamic  condition  then 
yields  the  most  probable  droplet  size  and  charge-to-mass  ratio. 

This  derivation  is  valid  when  the  spray  is  monodisperse,  but  the  production  of  such 
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sprays  is  extremely  difficult  in  practice.  Usually,  some  spread  is  observed  around  the  mean 
droplet  size.  Therefore,  this  model  can  only  give  rough  estimates  of  the  size  and  charge-to-mass 
ratio  due  to  the  simplifying  assumptions  made  in  the  derivation.  (The  minimum  energy  and  the 
electrohydrodynamic  model  are  presented  in  Appendix  A  and  B) 

The  maximum  entropy  model,  developed  by  Kelly2,  considers  the  electrostatic  spray  to 
be  an  isolated  thermodynamic  system.  At  equilibrium,  the  entropy  of  the  system  is  maximum. 
Since  the  entropy  of  a  system  is  proportional  to  its  degrees  of  freedom,  an  analytical  expression 
for  the  entropy  can  be  obtained  by  enumerating  all  the  possible  states  of  the  system.  The 
equilibrium  end  state  is  then  obtained  by  maximizing  the  entropy  with  the  system's  conservation 
equations  as  constraints  (i.e.  conservation  of  mass,  energy,  and  charge).  Additional  details  of 
this  theory  will  be  discussed  in  the  following  chapter.  A  further  point  to  be  made  for  the  three 
models  is  that  the  minimum  energy  and  the  electrohydrodynamic  models  predict  only  the 
charging  characteristic  of  the  sprays  while  the  maximum  entropy  model  is  capable  of  predicting 
both  the  charging  and  size  charactertistics. 

Despite  the  growing  interest  in  and  increasing  applications  of  electrostatic  sprays,  little 
fundamental  data  are  available,  and  very  few  fluids  have  been  studied.  The  major  difficulties  in 
obtaining  data  are  due  to  the  large  number  of  droplets  produced  in  an  electrostatic  spray  device 
and  their  generally  small  size  -  typically  in  the  micron  range.  To  produce  a  set  of  meaningful 
data,  it  is  necessary  to  make  charge  and  size  measurements  simultaneously  at  a  reasonable  rate. 
As  a  result,  most  studies  are  limited  in  scope.  Precise  measurements  of  charge-to-mass  ratio  and 
mass  of  a  single  droplet  are  sparse.  Experimentalists  often  have  had  to  resort  to  crude  estimates 
of  these  important  parameters  to  infer  the  characteristics  of  the  spray. 

Only  three  sets  of  data  for,  Octoil1,  Wood's  metal18,  and  glycerine19  have  been 
published.  Of  the  three  sets,  the  Octoil  data,  while  the  most  detailed,  is  still  far  from  complete. 
These  data,  published  by  Hendricks,  are  shown  in  Figure  1.2.  Each  point  on  the  figure 
represents  a  single  charged  droplet  measured.  This,  unfortunately,  represents  the  sole  data 
available.  To  obtain  a  better  understanding  of  the  electrostatically  driven  atomization  process,  a 
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CHARGE  —TO  -MASS  RATIO  (COULOMBS/ KILOGRAMS) 


Figure  1 .2  Hendricks'  Data 


more  comprehensive  data  base  is  required 


1.4  Summary  of  Results 

All  of  the  parameters  introduced  in  the  maximum  entropy  model  can  now  be  numerically 
evaluated.  Examination  of  these  results  shows  that  the  model  is  self-consistent.  Comparison  of 
the  present  data  with  existing  Hendrick's  data  reveals  that  the  finite  electrical  conductivity  of 
Octoil  can  cause  underestimations  of  droplet  radius.  A  correction  factor  based  on  conductivity 
variation  with  volumetric  charge  density  is  introduced  to  correct  the  error.  The  size  and  charge 
distribution  at  various  charge-to-mass  ratios  for  an  Octoil  spray  are  measured  in  this  research. 
Comparison  of  the  initial  Octoil  data  with  Hendrick's  shows  good  agreement. 

Regression  analysis  shows  that  the  experimentally  determined  size  distributions  correlate  . 
well  with  the  maximum  entropy  model.  The  close  agreement  confirms  the  multimodal  nature  of 
the  size  distribution  claimed  by  the  model.  This  indirectly  implies  the  existence  of  a  phase 
transition-like  behavior  in  the  charging  characteristic.  However,  more  conclusive  evidence  will 
require  the  examination  of  the  size  distribution  at  very  high  charge-to-mass  ratio  (>  5  C/kg) 
which  are  inaccessible  at  the  moment  due  to  instrumental  limitations.  Consistently  higher 
limiting  surface  electron  Held  emission  field  strengths  are  found  in  certain  distributions  of 
droplets.  This  result  suggests  possibly  different  surface  properties  of  the  droplets  in  these  spray 
distributions  and  should  be  the  focus  of  additional  research. 
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Chapter  II 

The  Maximum  Entropy  Model 


In  electrostatic  sprays,  the  quantities  that  are  of  interest  to  engineering  applications  are  the 
distributions  of  the  charge-to-mass  ratio,  the  charge,  and  the  size  of  the  droplets.  In  addition,  the 
variance  of  these  distributions,  the  spray's  penetration  depth,  and  the  size  range  are  important 
parameters  to  observe.  The  three  classes  of  models  mentioned  in  the  Chapter  I,  however,  do  not 
attempt  to  model  all  of  these  properties  of  the  spray.  Ultimately,  what  is  desired  in  a  spray  model 
is  its  ability  to  predict  these  properties  from  specification  of  just  a  few  parameters  such  as  sprayer 
voltage,  fluid  flow  rate,  and  fluid  properties. 

The  minimum  energy  model  describes  the  variation  of  charge-to-mass  ratio  as  a  function 
of  average  radius  in  a  charged  spray.  In  the  electrohydrodynamic  model,  the  charge-to-mass 
ratio  is  modeled  as  a  function  of  average  droplet  radius,  mass  flow  rate,  and  the  geometry  of  the 
sprayer.  In  the  maximum  entropy  model,  the  charging  characteristic  of  an  equilibrium  spray 
cloud  is  related  to  the  droplet  radius  and  other  fluid  properties.  Applying  the  same  formalism, 
Kelly  is  able  to  model  the  atomization  process  thereby  deriving  the  size  distribution  model. 
However,  none  of  these  models  has  yet  been  able  to  describe  the  distribution  of  the 
charge-to-mass  ratio  in  a  spray.  In  areas  such  as  space  propulsion  and  non-impact  printing, 
accurate  control  of  the  particle  beam's  charge-to-mass  ratio  is  crucial  for  the  purpose  of  particle 
focusing.  Thus,  attaining  a  model  capable  of  describing  the  charge-to-mass  ratio  distribution  in  a 
spray  is  the  next  step  in  the  theoretical  modeling  work. 

2.1  Maximum  Entropy  Model 

In  the  maximum  entropy  model,  no  effort  is  made  to  model  the  macroscopic  breakup  process 
during  spray  formation.  The  model  seeks  to  obtain  the  most  probable  size  distribution  of  a 
charged  spray  when  it  is  in  thermodynamic  equilibrium.  The  problem  may  be  posed  in  the 
following  way:  Given  a  fixed  mass  of  liquid  M,  total  energy  E  and  total  free  charge  Q,  what  is 
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the  most  likely  equilibrium  size  and  charge  distribution  when  the  liquid  is  atomized?  The  model 
takes  advantage  of  the  long  range  coulombic  interaction  between  charged  droplets  which 
produces  a  well-mixed  state  of  fluid  and  charge.  The  electrostatic  interaction  is  assumed  to 
promote  rapid  equilibration  within  die  spray. 

The  model,  developed  by  Kelly2,  assumes  that  two  independent  processes  are  occuring 
simultaneously;  the  atomization  and  charging  process.  From  this  assumption,  the  droplet  size 
distribution  can  be  modeled  independently  of  the  charge  distribution.  The  derivation  is  very 
similar  to  the  statistical  mechanical  derivation  of  the  most  probable  state  of  an  isolated  system  of 
ideal  gas  molecules. 

From  statistical  mechanics,  an  isolated  system  attains  thermodynamic  equilibrium  when 
its  entropy  is  maximized.  The  system  entropy  is  related  to  the  degrees  of  freedom  of  the  system. 
Thus,  to  obtain  an  expression  for  the  system  entropy,  one  must  be  able  to  write  an  analytical 
expression  for  the  degrees  of  freedom.  The  fundamental  expression  for  the  entropy  s,  is  given 
by  s=k  ln(w),  where  k  is  the  Boltzmann's  constant  and  w  is  the  degrees  of  freedom  of  a 
particular  system  state.  In  this  case,  a  system  state  is  identified  as  a  particular  distribution  of 
droplet  charge  and  size.  The  degrees  of  freedom  of  a  system  state  is  also  its  degeneracy.  For  the 
atomization  process  the  degeneracy  is  prescribed  by  the  Maxwell-Boltzmann22  statistics, 
while  the  degeneracy  of  the  charging  process  is  prescribed  by  Fermi-Dirac22  statistics.  These 
statistics  enable  us  to  express  w  analytically. 

From  s=k  ln(w),  we  can  find  the  most  probable  state  by  looking  for  a  state  of 
maximum  w  since  s  is  a  monotonic  increasing  function  of  w.  Mathematically,  the  most 
probable  state  is  determined  by  the  standard  Lagrangian  Multiplier  Method.  Such  a  method 
requires  a  constant  to  be  introduced  with  each  conservation  law  applied  to  the  system.  In  this 
case  three  constants  (or  Lagrange  Multipliers)  a',  (5',  and  S'  are  assigned  to  the  conservation  of 
charge,  energy,  and  mass  respectively.  The  resulting  set  of  equations  from  the  Lagrangian 
method  are  as  follows: 
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where  equation  (2.1)  is  an  implicit  equation  describing  the  charging  characteristic  and  equation 
(2.2)  describes  the  size  distribution.  These  equations  define  a  statistically  favorable  state  of  the 
observed  spray  cloud  at  equilibrium. 

In  the  maximum  entropy  model,  the  Rayleigh  Limit1 1  defines  a  non-dimensional  charging 
parameter,  Z,  which  is  the  ratio  of  the  number  of  droplet  surface  electrons  to  the  maximum 
number  allowed.  The  surface  electron  field  emission  strength,  Es,  identifies  another  limit  in  the 
electrostatic  spray.  The  charge  on  a  droplet  is  such  that  it  cannot  cause  a  surface  field  in  excess 
of  Es.  At  small  droplet  radius,  the  amount  of  charge  allowed  on  the  droplet  is  limited  by  the 
critical  emission  field  strength.  However,  at  larger  radius,  the  charge  on  the  droplet  is  limited  by 
the  Rayleigh  Limit.  This  is  because  the  Rayleigh  Limiting  charge  is  proportional  to  the  3/2 
power  of  the  radius  while  the  field  limiting  charge  is  proportional  the  square  of  the  radius.  Using 
Eg  and  the  surface  tension  y,  a  natural  scaling  length  can  be  introduced  into  the  model,  of  the 
form: 


r 
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where  rx  is  known  as  the  cross-over  radius.  Physically,  the  cross-over  radius  is  the 
characteristic  radius  of  a  charged  droplet  at  which  it  is  equally  likely  to  undergo  a  Rayleigh  type 
breakup  or  surface  electron  field  emission  to  relieve  its  surface  electrostatic  stress  (for  Octoil,  the 
cross-over  radius  is  approximately  3.3  x  10*4  microns).  The  droplet  radius  in  the  maximum 
entropy  model  is  normalized  with  the  cross-over  radius, 
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2.2  Preliminary  Results  of  Maximum  Entropy  Model 

The  maximum  entropy  model  was  compared  to  Hendrick’s  Octoil4  and  Wood's  metal18  data  by 
Kelly.  Because  experimental  size  distributions  were  not  available,  only  equation  (2.1)  was 
tested.  In  this  research,  the  theoretical  size  distribution  equation  (2.2)  is  checked  with  the 
experimental  size  distributions. 

Figure  2.1  depicts  the  data  for  Octoil  and  Wood's  metal,  each  data  point  on  the  graph 
represents  an  averaged  radius  from  the  group  of  droplets  detected  at  the  same  charge-to-mass 
ratio.  Using  these  data  and  equation  (2.1),  the  two  Lagrangian  multipliers  a'  and  p’  can  be 
determined  for  the  fluids.  Their  values  are  listed  below: 

Wood's  metal  a’  ■  -  0.18  p'  *  1.7  x  10l  6 

Octoil  a’  -  - 0.82  x  10- 4  p'  •  7.5  x  101 2 

Using  these  values  for  a'and  p',  the  charging  equation  (2.1)  can  be  solved  numerically.  It  turns 
out  that  for  both  Wood's  metal  and  Octoil,  a  critical  charge-to-mass  ratio  exists  below  which 
multiple  solutions  of  average  droplet  radius  are  found  for  a  given  charge-to-mass  ratio.  Above 
this  critical  ratio,  a  single  solution  for  average  radius  is  found  for  each  charge-to-mass  ratio  of  the 
charging  equation.  Each  of  these  solutions  for  the  radius  represents  the  averaged  radius  of  an 
individual  distribution  of  droplets.  Thus,  for  the  multiple  solution  case,  one  would  expect  the 
observed  distribution  to  be  composed  of  three  smaller  distributions. 

Figure  2.1  shows  the  solution  for  the  charging  equation  that  best  Fits  the  published 
experimental  data.  The  volumetric  charge  density  is  plotted  against  the  radius  because  it  is  a 
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more  appropriate  unit  of  engineering  interest  (the  volumetric  charge  density  is  defined  as  the 
product  of  the  charge-to-mass  ratio  and  the  fluid  density).  This  solution  branch,  for  the  moment, 
is  designated  as  the  B  Branch.  The  other  solutions  are  each  designated  as  the  A  and  C  Branch 
respectively. 

The  B  Branch  is  of  special  interest  because  on  it,  the  charging  equation  (2.1)  may  be 
simplified  for  large  droplets.  For  the  observed  data  of  Octoil  and  Wood's  metal,  the  assumption 
of  large  droplets  holds  for  the  range  of  sizes  observed30.  This  simplification  allowed  the 
charging  equation  to  be  written  into  the  form: 

,  .  12  e 

<3  -  -  1&-) - -  (2.5) 

where  d  is  the  diameter  of  the  droplet,  q  is  the  elemental  charge,  eQ  is  the  permittivity  of  free 
space  and  pe  is  the  volumentric  charge  density. 

Figure  2. 1  reveals  that  the  simplified  charging  equation  (2.5)  correlates  well  with  the 
published  data.  Therefore,  this  equation  can  predict  the  average  droplet  size  when  the  volumetric 
charge  density  is  specified.  For  hydrocarbon  fluids  such  as  diesel,  the  "large"  radius  regime  that 
facilitate  the  simplification  corresponds  to  droplet  sizes  of  10  microns  and  above.  These  are  sizes 
of  special  interest  in  engineering,  particularly  in  the  area  of  combustion  research. 

Equation  (2.5)  relates  the  diameter  of  the  droplet  to  its  volumetric  charge  density.  For  a 
given  fluid,  the  ratio  (-  a'/p')  is  a  constant  Two  points  were  stressed  in  Kelly’s  work  on  this 
simplified  charging  equation.  Firstly,  the  value  -  a'/p'  evaluated  for  Wood's  metal  and  Octoil 
were  very  close.  This  is  counter  intuitive  because  the  two  fluids  are  very  different  in  terms  of 
their  densities,  viscosities,  and  surface  tensions,  furthermore,  Wood's  metal  is  electrically 
conducting  while  Octoil  is  a  dielectric  fluid.  The  physical  significance  of  -  a'/p'  should  therefore 
deserve  further  examination.  Secondly,  equation  (2.5)  is  not  a  function  of  fluid  properties.  This 
is  in  contrast  to  the  case  of  uncharged  sprays  where  the  spray  properties  are  dependent  on  the 
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Figure  2.1  Octoil  and  Wood's  Metal  Data 
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properties  of  the  fluid  being  atomized. 

Therefore,  at  large  droplet  radius,  the  size  of  the  droplet  can  be  correlated  to  the 
volumetric  charge  density  alone.  For  conducting  fluids  like  Wood’s  metal  more  charge  can  be 
injected  into  the  fluid.  This  causes  a  larger  electrostatic  repulsive  force  on  the  surface  and  the 
Rayleigh  Limit  is  attained  at  a  smaller  radius;  as  a  result,  smaller  droplets  are  observed.  For 
Octoil,  charge  injection  is  not  as  efficient  due  to  its  dielectric  properties.  This  causes  larger 
droplets  and  therefore  lower  charge  densities  to  be  observed  in  the  spray. 


2*3  Further-AkrificattQn  of  the  Maximum  Entropy  Model 

We  have  seen  that  usii.g  the  existing  data,  only  two  of  the  three  Lagrangian  Multipliers  can  be 
evaluated.  To  obtain  further  information  on  the  model,  the  published  data  are  insufficient.  From 
this  facility,  full  characterization  of  the  spray  can  be  obtained  in  terms  of  detailed  size  and  charge 
distributions.  This  not  only  allows  us  to  obtain  a  more  accurate  average  of  the  previous  data,  it 
also  facilitate  the  testing  of  the  size  distribution  model.  By  appealing  to  the  size  distribution 
model  and  the  size  distribution  data,  we  can  proceed  to  evaluate  the  5'  multiplier.  In  addition,  a’ 
and  P'  can  also  be  found  from  this  exercise  thereby  allowing  comparison  of  their  values  with 
those  obtained  from  Kelly's  work. 
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Chapter  III 

The  Electrostatic  Atomization  Experiment 


This  chapter  describes  the  experimental  facility  constructed  to  collect  electrostatic  spray  data.  The 
results  and  predictions  of  the  maximum  entropy  model  is  further  examined  in  the  light  of  the  new 
data.  Data  for  Octoil  from  this  experiment  is  tested  against  Hendrick's;  this  enables  a  calibration 
of  the  facility.  Other  fluids  can,  in  turn,  be  sprayed  so  as  to  build  up  a  detailed  set  of 
fundamental  spray  data. 

The  facility  described  here  resembles  closely  to  the  Hogan  and  Hendrick's  apparatus20. 
A  schematic  diagram  of  the  experiment  and  a  photograph  of  the  experimental  set-up  are  shown  in 
Figure  3. 1  and  3.2.  The  electrostatic  spray  is  examined  at  various  charge-to-mass  ratios  with  the 
quadrupole  mass  spectrometer.  Droplets  with  the  pre-selected  charge-to-mass  ratio  are  extracted 
from  the  main  droplet  beam.  The  free  electron  charges  on  these  droplets  are  then  determined 
with  a  high  speed  electron  sensitive  pre-amplifier24.  These  charge  data  are  cumulatively  collected 
on  a  multichannel  analyzer  (MCA)  which  yields  the  charge  distribution  of  the  spray  at  the 
selected  charge-to-mass  ratio.  By  assuming  that  the  droplets  are  spherical  the  size  distribution  can 
in  turn  be  calculated. 

A  vacuum  facility  is  required  for  the  operation  of  the  quadrupole  mass  spectrometer.  The 
complete  sprayer-mass  spectrometer-charge  analysis  system  is  maintained  under  a  vacuum  of  1  x 
10‘6  torr.  Data  collected  on  the  MCA  are  transferred  via  fiber  optic  data-link  to  a  HP9816 
personal  computer  for  preliminary  reduction.  This  processed  data  is  then  sent  to  the  Princeton 
University  mainframe  IBM  3081  computer  for  regression  analysis.  The  purpose  of  the 
regression  analysis  is  to  examine  the  goodness-of-fit  of  the  model  with  the  measured  data. 

3.1  Charged  Droplet  Production/Vacuum  Facility 

Charged  droplets  are  produced  using  the  simple  diode-type  sprayer  shown  in  Figure  3.3.  The 
sprayer  consists  of  a  metallic  capillary  needle  of  inner  diameter  0.01  "(1/4  mm)  and  outer  diameter 
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Figure?  3.1  Schematic  of 
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Figure  3.3  Diagram  of  Reservoir  and  Sprayer 


0.02"(l/2  mm),  inserted  into  a  plexiglass  reservoir.  The  flow  rate  of  the  sprayer  can  be 
controlled  by  regulating  the  pressure  over  the  fluid  in  the  reservoir,  pressures  up  to  760  tore  can 
be  applied.  The  capillary  needle  is  connected  to  a  high  voltage  source  capable  of  voltages  up  to 
20  kV. 

To  achieve  atomization  the  potential  of  the  needle  is  raised  to  14  kV  with  a  reservoir 
pressure  of  45  tore.  Since  it  is  the  electric  field  that  is  responsible  for  the  atomization,  a 
electrically  grounded  ring-electrode  (also  known  as  an  extractor  electrode)  is  positioned  coaxially 
around  the  needle  to  provide  the  required  field  for  spraying.  To  prevent  corona  discharge  at  the 
needle  tip,  an  insulating  jacket  made  from  heat  shrink  tubing  is  placed  over  die  entire  length  of 
the  capillary  needle.  It  is  found,  through  extensive  tests,  that  the  sprayer  works  most 
satisfactorily  if  this  jacket  is  extended  over  the  tip  of  the  needle  by  a  millimeter  or  two.  This 
enables  a  column  of  fluid  to  insulate  the  needle  tip  preventing  breakdowns.  Smooth  operation  of 
the  sprayer  requires  careful  adjustments  of  the  flow  rate  (i.e.  reservoir  pressure)  and  the  charging 
voltage.  If  the  charging  voltage  is  too  high  or  if  the  flow  rate  is  too  low,  breakdown  will  occur 
due  to  the  depletion  of  spray  fluid  which  exposes  the  sharp  needle  tip.  Breakdowns  may  also 
occur  in  the  dielectric  fluid  if  the  electric  field  at  the  tip  is  too  high.  Considerable  effort  has  been 
devoted  to  avoiding  arcing  since  it  damages  the  surface  of  the  capillary  needle  and  causes  false 
recordings  in  the  multichannel  analyzer. 

The  charged  droplets  are  accelerated  by  the  extractor  electrode  toward  the  quadrupole 
mass  spectrometer.  Droplets  are  typically  accelerated  to  a  velocity  of  150  -  200  m/s.  The  8mm 
diameter  entrance  aperture  to  the  quadrupole  limits  the  radial  velocity  component  of  the  charged 
droplets  thereby  ensuring  proper  filtering  in  the  mass  spectrometer34*35. 

3.2  Quadrupole  Mass  Spectrometer 

Two  pieces  of  information  are  required  to  characterize  a  single  charged  droplet,  namely  its  charge 
and  size.  From  these  data,  other  quantities  such  as  mass,  surface  electron  density,  electrostatic 
stress,  and  surface  energy  can  be  derived.  Early  versions  of  electrostatic  spray  experiment  utilize 
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time  of  flight  method  to  deduce  these  quantities9.  This  method  is  admissible  when  the  charge  on 
the  droplet  is  relatively  large  because  at  low  charge  levels  the  measurements  have  very  small 
signal-to-noise  ratios.  Since  the  droplets  of  engineering  interest  have  low  charge  and 
charge-to-mass  ratios,  other  means  of  measurement  must  be  used  to  obtain  accurate  data. 

An  alternative  is  to  use  a  quadrupole  mass  spectrometer.  By  using  the  quadrupole  as  a 
mass  filter  only  droplets  with  the  selected  charge-to-mass  ratio  are  filtered  from  the  droplet  beam. 
If  the  charge  of  these  droplets  are  subsequently  measured,  other  information  can  be  computed 
(e.g.  droplet  size  and  surface  energy).  Since  the  quadrupole  can  be  operated  continuously,  less 
time  is  required  to  obtain  a  meaningful  charge  distribution  than  the  time  of  flight  instrument. 
Detailed  operating  principles  and  design  criterion  of  the  quadrupole  mass  spectrometer/jpower 
supply  used  in  this  research  is  described  in  Appendix  C. 

3.3  Charge  Detection  Assembly 

The  charge  analysis  system  consists  of  an  electrometer-amplifier  combination  manufactured  by 
Amptek  Inc.  (model  number  A-225  and  A-206  respectively).  The  electrometer  has  a  charge 
sensitivity  threshold  of  280  electrons.  For  each  charged  droplet  analyzed,  a  voltage  pulse  is 
generated  with  a  characteristic  pulse  width  of  10  jisec  and  a  rise  time  of  2.5  p.sec.  The  pulse 
height  is  proportional  to  the  amount  of  charge  detected  in  the  droplet  The  output  voltage  pulses 
are  recorded  on  a  multichannel  analyzer  (MCA)  which  then  displays  a  population  versus  droplet 
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charge  distribution.  Since  a  minimum  threshold  voltage  is  required  on  the  MCA;  therefore,  small 
signals  from  the  electrometer  must  be  amplified.  This  unfortunately  also  amplifies  die  noise  in 
the  system. 

While  the  droplet  detection  plate  is  being  partially  shielded  from  the  quadrupole  as  shown 
in  Figure  3.1,  it  is  not  sufficient  to  eliminate  completely  its  noise  pickup  from  the  quadruple's 
r.f.  field.  Because  the  applied  frequency  to  the  quadrupole  rods  is  known,  a  filter  is  built  to 
remove  this  "noise"  contribution.  Other  random  noise  in  the  system  cannot  be  completely 
eliminated.  This  problem  is  particularly  serious  at  high  gain  settings  on  the  amplifier.  This 
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difficulty  prohibits  the  collection  of  data  at  charge  levels  below  8000  electrons.  Under  this 
restriction,  the  smallest  droplet  size  detectable  is  about  one  micron.  Fortunately,  this  is  still 
within  the  range  of  interest 
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In  a  typical  experiment,  the  quadrupole  voltage  signals  are  calibrated  first.  This  involves 
balancing  die  two  a.c.  voltage  amplitudes  and  verifying  their  phase  difference.  A  variable 
resistor  is  used  to  balance  their  amplitudes.  The  required  d.c.  bias  voltages  are  computed  from 
the  program  "MCAHP_SM"  and  fed  in  with  the  a.c.  signals.  These  voltages  are  verified  from  a 
digital  oscilloscope.  Due  to  the  bit  resolution  on  die  oscilloscope,  an  error  between  1  to  3  %  is 
usually  associated  with  these  signals. 

Upon  completing  the  voltage  verification,  the  sprayer's  fluid  reservoir  is  pressurized  to 
approximately  40  torr.  This  provides  a  positive  flow  rate  through  the  10  mil  (i.d.)  capillary 
needle.  The  needle  charging  voltage  is  increased  slowly  in  small  increments  to  permit 
equilibration.  The  quadrupole  is  de-activated  during  this  initial  stabilization  process  to  evaluate 
the  charge  detector  count  rate,  which  provides  a  qualitative  measure  of  the  sprayer  output 

Past  experience  indicates  spraying  in  vacuum  starts  at  approximately  13  kV.  For  this 
sprayer-extractor  electrode  configuration,  voltages  beyond  14  kV  produces  arcing,  which  if 
sustained  can  demage  the  capillary  needle  and  the  electronics  of  the  MCA. 

In  a  typical  run,  the  sprayer  is  allowed  to  stabilize  for  a  few  minutes  before  the  mass 
spectrometer  is  turned  on.  The  droplet  detection  rate  usually  drops  by  a  factor  of  ten  when  the 
mass  spectrometer  is  activated  because  unwanted  droplets  are  being  filtered  off.  Approximately 
1000  droplets  are  required  to  produce  a  useful  distribution,  this  usually  requires  3  to  5  hours  of 
experiment  time. 


Chapter  IV 

Data  Analysis  and  Discussion 


4.1  Introduction 

The  size  distribution  of  the  spray  is  experimentally  determined  at  various  charge-to-mass  ratios. 
Figure  4.1  shows  the  capability  of  the  experiment  in  obtaining  such  data.  The  scanning 
charge-to-mass  ratio  of  the  quadrupole  mass  spectrometer  is  set  at  a  value  and  bandpass  ratio. 
The  resultant  MCA  spectrum  would  then  correspond  to  a  narrow  horizontal  band  extending 
rightward  from  the  detector  sensitivity  limit  in  Figure  4.1.  Plotted  in  Figure  4.1  also  are  the 
theoretical  results  of  equation  (2.1).  The  first  order  phase-transition  behavior  can  be  seen  when 
the  equation  transits  from  single  to  multiple  solutions. 

For  Octoil,  the  transition  occurs  approximately  at  3.4  C/kg.  Since  we  cannot ,  at  Jie 
moment,  examine  charge-to-mass  ratios  greater  than  2.9  C/kg  we  should  therefore  anticipate 
seeing  multimodal-type  size  distributions.  It  is  important  to  realize  that  the  curves  of  Figure  4. 1 
represents  the  behavior  of  mean  droplet  size,  thus,  the  observed  MCA  spectra  should  consist  of 
more  than  one  distribution  centered  on  each  these  means. 

4.2  Experimental  Data 

A  typical  spectrum  taken  from  the  MCA  is  shown  in  Figure  4.2  where  the  y-axis  is  the 
population  number  and  the  x-axis  is  proportional  to  the  amount  of  charge  detected  from  the 
incoming  droplet.  Droplets  of  low  charge  levels  usually  have  a  large  population.  It  then  drops 
off  rather  rapidly  at  higher  charge  levels.  Hie  number  of  data  channels  defining  the  peak  of  the 
distribution  is  small  in  comparison  with  the  total  number  of  channels  in  the  distribution.  This 
implies  that  the  distributions  have  rather  sharp  peaks  and  is  in  constrast  with  the  rather  broad  and 
smooth  size  distributions  seen  in  uncharged  sprays.  These  spectra  also  indicate  required 
improvments  to  be  made  on  low  charge  level  data  acquisition. 

At  higher  charge-to-mass  ratios,  typically  near  the  present  limiting  charge-to-mass  ratio  of 
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Figure  4. 1  Capability  of  Experiment 


Raw 
(  Oc to i  1 


s}unc>3  j.o  •  Ofvj 


Figure  4.2  Raw  Data  from  Multichannel  Analyzer  (MCA) 


2.9  C/kg,  some  detected  droplets  were  found  to  possess  charging  levels  well  in  excess  of  the 
Rayleigh  limit  determined  by  the  quadrupole  charge-to-mass  setting.  This  suggests  that  a  portion 
of  the  droplets  have  not  been  effectively  filtered  by  the  quadrupole  field.  However,  this 
phenomenon  is  only  seen  when  the  quadrupole  is  operating  at  its  upper  limit  Through  more 
stable  power  supply,  this  deficiency  can  be  eliminated.  The  small  number  of  droplets  having 
anomalous  charging  levels  has  no  impact  on  the  proper  interpretation  of  the  data  as  they  are 
usually  small  in  population  and  are  not  statistically  influential. 

4.3  Comparison  with  Published  Work 

Collected  from  the  experiment  is  a  set  of  size  distribution  at  various  charge-to-mass  ratios.  As  a 
first  test  of  the  experiment's  consistency,  the  MCA  spectrum  data  were  simply  averaged  in 
accordance  with  the  procedure  used  by  Hendricks  to  represent  his  manually  recorded  data. 
Figure  4.3  shows  a  comforting  level  of  agreement  of  the  initial  results  with  Hendrick's 
pioneering  work.  This  provides  confidence  that  the  experiment  is  operating  in  an  appropriately 
consistent  manner. 

A  slight  deviation  between  the  two  data  set  is  apparent  at  low  charge-to-mass  ratios. 
Analysis  of  these  data  indicate  that  this  deviation  may  be  explained  by  the  finite  conductivity  of 
Octoil  which  prevents  the  proper  measurement  of  droplet  charge  within  the  electrometer  response 
time  of  .002  msec.  It  is  computed  that  only  a  fraction  of  the  total  charge  can  be  transfered  if  the 
charge  density  of  the  droplet  is  low  resulting  in  low  electrical  conductivity.  As  a  result  of  the 
false  measurement  of  charge,  the  computed  size  would  be  in  error  since  the  size  of  the  droplet  is 
computed  from  the  knowledge  of  the  charge-to-mass  ratio  and  measured  total  charge  of  the 
droplet.  In  section  4.8,  a  method  of  correction  will  be  d'scussed. 

New  data  (shown  in  Figure  4.3)  is  checked  with  the  charging  equation  (2.1)  derived 
from  the  maximum  entropy  formalism.  Again,  of  the  three  solutions  predicted,  the  B  branch  was 
found  to  have  the  closest  correlation  with  the  data.  The  simplified  form  of  die  charging  equation 
on  this  branch  is  given  by  equation  (2.5).  Using  regression  analysis  the  constant  (-  a'/f$')  is 
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Figure  4.3  Comparison  with  Published  Data 


Droplet  Radius  (Microns) 


re-evaluated.  The  value  deduced  using  the  new  data  is  -53.7  eV  which  a  approximately  15% 
from  Kelly's  earlier  value.  The  relatively  large  difference  from  Kelly's  value  arises  because 
equation  (2.5)  is  not  very  sensitive  to  values  of  the  ratio  (-  a'/jj').  Furthermore,  there  are 
probable  errors  inherent  in  the  original  Hendricks  data  which  could  account  for  the  difference. 

4.4  ^ize  Distribution 

The  size  distribution  model  was  not  validated  due  to  the  lack  of  size  distribution  data  in  the 
literature.  With  the  experimental  charge  distribution  data,  it  is  possible  to  deduce  the  size 
distribution  since  the  charge-to-mass  ratio  at  which  the  chaige  distribution  is  obtained  is  known. 

The  size  distribution  model  is  specified  through  equation  (2.2).  However,  we  recall  that 
from  the  solution  of  the  charging  equation,  the  size  distribution  should  be  multimodal.  This  is  to 
say  that  the  observed  size  distribution  should  be  a  convolution  of  three  smaller  distributions.  The 
actual  size  distribution  model  is  then  a  summation  of  three  equations  each  with  the  form  of 
equation  (2.2): 


the  summation  represents  the  sum  of  the  three  branches  A,  B,  and  C.  Using  equation  (4.1)  and  a 
standard  non-linear  regression  routine,  the  size  distribution  model  is  tested  with  the  data.  The 
regression  routine  is  supplied  by  SAS  Institute  Inc.  and  implemented  on  the  Princeton  University 
IBM  370  mainframe  computer.  The  routine  uses  a  non-linear  Gauss-Newton  Least  Square 
method  in  its  analysis. 

Equation  (4.1)  is  a  function  of  the  charging  parameter  Z  and  the  normalized  radius.  The 


variable  Z  can  be  eliminated  by  relating  it  to  the  charge-to-mass  ratio  of  the  droplets.  By  writting 
the  charge-to-mass  ratio  as: 
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The  charging  parameter  may  then  be  introduced  by  dividing  the  numerator  and  denominator  of 
the  right-hand-side  by  Nray,  which  is  the  number  of  electrons  a  droplet  contains  at  the  Rayleigh 
Limit  at  a  given  radius.  Equation  (4.2)  thus  becomes: 
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The  advantage  of  using  the  above  expression  is  therefore  apparent  because  the 
experimental  data  have  fixed  charge-to-mass  ratios.  Using  equation  (4.3)  it  is  then  possible  to 
eliminate  one  of  the  independent  variables  in  the  size  distribution  model.  In  this  case,  the  size 
distribution  becomes  a  function  of  the  normalized  radius  only.  The  new  model  equation 
therefore  becomes: 
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4.5 .  Regression  Analysis,  Results 

The  regression  analysis  program  uses  equation  (4.4)  as  an  input  model  to  examine  its  degree  of 
correlation  with  the  observed  size  data.  All  the  MCA  spectra  examined  indicate  the  absence  of 
the  C  Branch.  In  other  words,  the  experimental  data  could  only  be  resolved  into  two  individual 
distributions.  These  distributions,  when  compared  to  the  theoretical  results  were  identified  to  be 
the  A  and  B  Branch  respectively.  While  it  is  unclear  why  the  C  Branch  is  absent,  it  is  possible 
that  the  C  Branch  population  is  too  low  to  be  resolved  by  the  analysis  routine. 

Figure  4.4  depicts  the  regression  result  of  one  of  the  observed  size  distributions.  Three 
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curves  are  generated  in  each  of  these  analyse.  The  curves  labeled  A  and  B  are  the  respective 
branches,  while  the  third  curve  is  the  sum  of  the  previous  two.  The  correlation  in  genual  is  very 
good  as  indicated  by  four  statistical  indicators  for  measuring  the  goodness*of-fit  of  the  theoretical 
model. 

The  four  indicators  are  the  correlation  coefficient,  the  sum  of  errors,  the  residual  plot  and 
the  standard  errors  of  the  estimated  model  coefficients.  Aside  from  two  cases,  all  the  size  spectra 
examined  yield  conelatior  coefficients  well  over  0.9.  This  shows  that  over  90%  of  the  data  may 
be  accounted  for  by  the  proposed  model.  The  effectiveness  of  the  model  is  further  ensured  by 
good  behavior  in  other  statistical  indicators. 

These  results  concluded  that  the  observed  spray  size  distribution  is  indeed  multimodal. 
The  possibility  of  a  monomodal  size  distribution  is  being  ruled  out  through  a  direct  comparison 
of  such  a  model  with  experimental  data.  Figure  4.5  shows  that  this  model  is  incapable  of 
reproducing  the  data.  In  fact,  the  monomodal  model  has  a  correlation  coefficient  of  about  0.7 
with  the  data.  This  result  reaffirms  the  multimodal  property  of  the  size  distribution. 

4.6  Multiple  Charging  Characteristics 

From  the  statistically  resolved  A  and  B  branch,  it  is  now  possible  to  compare  the  experimental 
charging  characteristics  to  the  theoretically  predicted  ones.  This  is  done  in  Figure  4.6  and  4.7. 
Discrepancies  were  observed  in  both  figures.  The  data  cannot  be  predicted  completely  from  the 
model  equations.  Gose  scrutiny  of  the  discrepancies  reviewed  that  the  deviation  between  theory 
and  data  were  not  constant  but  dependent  on  the  charge-to-mass  ratio.  Such  behavior  was  also 
observed  in  the  comparison  with  Hendrick's  data  in  the  earlier  section.  In  the  branch  A 
characteristic,  the  theory  tends  to  underpredict  the  data  while  the  opposite  occurs  for  the  branch  B 
characteristic. 

The  differences  between  theory  and  data  may  be  a  result  of  the  varying  electrical 
conductivity  of  Octoil  at  different  charge-to-mass  ratios.  This  idea  is  examined  through  the 
assumption  that  the  electrical  conductivity  of  Octoil  consists  of  two  components.  One  of  which 
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is  the  residual  conductivity  of  pure  Octoil  and  the  other  contributed  by  the  presence  of  excess 
charge  carriers  (i.e.  free  electrons).  As  noted  in  section  4.3,  this  idea  could  explain  the  deviation 
of  the  data  which  is  charge-to-mass  ratio  dependent.  A  conductivity  effect  correction  is  then 
introduced  to  the  experimental  charging  characteristics  in  Figure  4.6  and  4.7  to  check  the 
assumption. 

The  results  of  this  correction  are  plotted  in  Figure  4.8  and  4.9.  For  both  experimental 
characteristics,  the  correction  has  improved  the  correlation  with  the  model.  Aside  from  a  biased 
shift,  both  the  theoretical  and  experimental  charging  characteristics  have  the  same  dependence  on 
the  droplet  radius.  The  question  which  remains  is  the  cause  of  the  biased  shift  that  resulted  in  the 
non-overlapping  of  the  experiment  and  theoretical  charging  characteristics.  This  could  be 
resolved  by  using  different  values  for  Es  in  the  solution  of  the  charging  equation.  It  may  be 
recalled  that  to  solve  the  charging  equation,  several  constants  must  be  known  and  they  art  the 
two  Lagrange  multipliers  and  Es,  the  limiting  surface  emission  field  strength.  The  Lagrange 
multipliers  were  evaluated  by  making  comparison  to  the  Hendrick's  data  while  the  value  for  Es  is 
an  estimated  quantity. 

The  value  of  Es  is  difficult  to  obtain  because  surface  emission  depends  on  a  variety  of 
parameters  such  as  surface  conditions,  ambient  gas  properties,  and  geometry  of  the  surface 
concerned.  Experiments  conducted  to  determine  us  value  is  usually  difficult  to  duplicate  and  the 
values  fluctuate  from  test  to  test  because  it  is  not  possible  to  control  experimental  conditions 
precisely.  A  study  was  performed  to  analyze  the  effect  of  Es  on  the  charging  characteristics. 
This  is  accomplished  by  changing  the  value  of  Es  and  solving  for  the  charging  characteristic  for 
each  of  the  values32.  These  calculation  indicated  that  if  a  smaller  value  of  Es  is  assumed,  the 
resulting  charging  characteristic  will  be  shifted  to  higher  radius;  the  opposite  will  occur  if  a  larger 
value  of  Es  is  used.  Therefore,  the  discrepancies  observed  in  Figure  4.8  and  4.9  could  be  a 
result  of  incorrect  values  of  Es  used  in  solving  the  charging  equation. 

Since  the  value  of  Es  should  be  a  constant  for  a  given  surface  characertistic,  difference 
in  Es  could  mean  that  droplets  on  the  branches  are  exhibiting  different  surface  properties.  In 
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particular,  a  difference  in  dielectric  constant  would  account  for  the  observed  change  in  the  value 

of  Es. 

Since  the  size  distributions  have  been  resolved  into  two  separate  distributions,  it  is  now 
possible  to  estimate  the  Lagrange  multipliers  for  each  of  these  individual  distributions.  This 
provides  a  check  to  the  values  evaluated  from  the  averaged  data  in  the  earlier  section.  Using  the 
two  charging  characteristics  and  equation  (2.2),  a  regression  analysis  routine  allows  us  to  deduce 
a'  and  p'  for  branch  A  and  branch  B. 

Branch  £L_1  ILL  -ft.’/B' — (e,Y) 

A  -9.75  x  10-5  1.75  x  1013  34.8 

B  -4.29  x  10-5  4.92  x  1012  54.4 

The  ratio  a'/P'  for  the  B  Branch  is  very  close  to  the  value  53.7  eV  determined  from  the  averaged 
data  point  in  Section  4.3.  The  reason  for  this  is  a  result  of  the  dominant  population  of  the  B 
branch  droplets  compared  to  the  A  branch. 


4.7  Regression  Coefficients 

The  size  distribution  equation  is  given  by  equation  (4.4).  However,  for  the  purpose  of 
regression  analysis,  the  coefficients  are  lumped  into  a  single  number  for  ease  of  analysis.  The 
mathematical  form  of  the  model  used  in  the  computer  analysis  is  then  given  by, 

c  c 

N'ZrfNdx2rfOCp{Z2r  +  Z3r  +  Z5r)  <4-5) 

A  A 

where  the  Z's  are  the  lumped  coefficients.  The  regression  routine  then  determines  the  values  of 
these  coefficients  that  provide  the  best  fit  of  the  model  to  the  data.  Knowing  the  values  of  these 
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coefficients,  it  is  therefor  possible  to  determine  some  parameters  still  unknown  in  the  model,  for 
example  the  5'  Lagrange  multiplier.  Comparing  equation  (4.5)  and  equation  (4.4),  we  can  show 
that  the  Z  coefficients  have  the  following  forms. 


(4.6) 

£oE, 

Zj  -  4)  (-£I4)[  (-)  (£)  +  8'  ]  (4.7) 

“  3  t2J  q  M 

eo  $ 

Zj  *  (~)  P‘  1  }  (“■)  (4.8) 

J  9  5  c7  M 

E. 


Knowing  the  values  of  the  Z  coefficients  and  equation  (4.6),  (4.7)  and  (4.8),  all  three 
Lagrange  multipliers  and  the  value  of  Es  can  be  calculated.  Since  we  are  able  to  determine  two  of 
the  Lagrange  multipliers  from  the  charging  equation,  this  evaluation  should  provide  a  check  on 
the  accuracy  of  the  multipliers  and  the  consistency  of  the  model  equations.  The  valu.  of  the  third 
Lagrange  multiplier  8'  in  the  Kelly  model  is  evaluated  for  the  first  time  using  the  present  data. 
Through  manipulation  of  the  equations,  values  for  Es  can  be  deduced  independently.  At  the 
present  stage  of  the  theory,  there  are  no  analytical  forms  for  these  Lagrange  multipliers. 
Therefore,  their  numerical  values  on  the  different  branches  are  valuable  to  the  further 
understanding  of  the  model. 

Taking  the  ratio  2^  to  Z5,  the  following  expression  for  Es  may  be  obtained. 


This  quantity  should  be  a  constant.  However,  calculation  from  equation  (4.9)  indicated  that  its 
value  is  distinctively  larger  for  the  8  branch.  This  is  consistent  with  the  result  of  the  previous 
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section  where  we  could  only  find  close  correlation  with  the  model  if  the  Es  is  larger  on  the  B 
branch.  As  pointed  out  in  the  previous  section  a  different  surface  condition  such  as  a  difference 
in  dielectric  constant  could  possibly  change  the  value  of  Es.  Thus,  the  calculated  value  for  Es 
might  be  viewed  as  an  "effective"  limiting  surface  electron  field  emission  field  strength. 

Knowing  the  value  of  Es,  it  is  then  possible  to  deduce  the  value  of  P'  from  either 
equation  (4.6)  or  (4.8).  Equation  (4.8)  is  chosen  in  this  study.  For  each  size  distribution,  two 
values  of  P'  can  be  evaluated  corresponding  to  the  two  branches  of  droplets.  Figure  4.10  shows 
that  there  is  considerable  scatter  in  the  A  branch  values  while  the  B  branch  increases  steadily. 
Since  P'  has  dimensions  of  inverse  energy,  increasing  values  of  P'  on  the  B  branch  indicates  a 
decrease  in  "characteristic  energy"  of  the  spray.  The  precise  implication  of  this  energy  decrease 
can  only  be  studied  when  an  analytic  form  for  P'  is  available. 

Equation  (4.7)  is  a  coupled  equation  for  a'  and  S'.  Therefore  it  is  not  possible  to  deduce 
individual  values  of  a'  and  S'  for  each  charge-to-mass  ratio.  However,  a  single  value  for  a'  and 
S'  could  be  deduced  statistically  from  the  behavior  of  Z3.  a'  is  found  to  have  a  value  of  2.97  x 
10-4  for  branch  A  and  -4.46  x  1(H  for  branch  B.  These  values  are  three  to  five  times  larger  than 
the  previously  established  value  of  -0.82  x  10*4.  The  value  for  8’  is  however  quite  constant  for 
both  branches  having  the  values  5.02  x  1014  for  branch  A  and  3.47  x  1014  for  branch  B. 

The  results  of  this  section  draw  us  to  the  following  conclusion:  In  the  Kelly  model  it  is 
assumed  that  an  electrostatic  spray  can  be  considered  an  isolated  system  in  thermodynamic 
equilibrium.  From  which  a  state  of  maximum  entropy  must  be  attained.  The  soundness  of  this 
assumption  has  been  ascertained  from  the  various  analysis  performed  on  the  model  against 
measured  data.  The  model  equations  have  showed  consistency  and  good  correlation  with  the 
Octoil  data.  The  maximum  entropy  formalism  is  therefore  a  valid  approach  to  electrostatic  spray 
modelling. 

4.8  Correction  Term  For  Finite  Electrical  Conductivity 

Examination  of  Hendrick's  and  the  present  data  set  showed  discrepancies  that  are  dependent  on 
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Figure  4.10  Variation  of  Lagrangian  Multiplier  Beta' 
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the  charge-to-mass  ratio.  The  deviation  between  the  two  is  particularly  prominant  at  low  values 
of  the  charge-to-mass  ratios.  Experimentally,  the  two  experiments  differ  primarily  in  their  charge 
detection  devices.  In  Hendrick's  case,  an  open  Faraday  cage  is  used.  Charge  on  the  detected 
droplet  are  inferred  from  a  cathode  follower's  induced  voltage  while  in  the  present  experiment, 
an  electrometer  is  used.  In  this  case,  charge  is  measured  directly  when  they  are  transferred  to  the 
electrometer.  This  poses  no  problem  when  the  electrical  conductivity  of  the  Octoil  droplet  is  high 
however  it  is  not  the  case  for  low  charge-to-mass  ratio  droplets  because  the  electrical  conductivity 
of  the  Octoil  is  largely  responsible  by  the  presence  of  free  electrons.  When  conductivity 
becomes  very  small  at  small  values  of  charge-to-mass  ratios,  the  charge  transfer  process  is 
slowed  down.  Because  the  electrometer  allows  only  a  specific  time  for  charge  transfer,  errors  is 
introduced  if  the  charge  relaxation  time  is  comparable  to  the  charge  transfer  time. 

A  correction  factor  is  introduced  in  this  section  to  account  for  the  finite  charge  relaxation 
time  in  Octoil.  The  amount  of  correction  is  charge-to-mass  ratio  dependent  and  is  negligible  at 
high  charge-to-mass  ratios.  This  is  consistent  with  the  fact  that  the  observed  deviation  from 
Hendrick's  data  is  only  seen  at  small  values  of  charge-to-mass  ratios.  In  this  respect,  it  is 
assumed  that  Hendrick's  setup  is  not  affected  by  the  Finite  electrical  conductivity  because  the 
measurement  is  inductive  rather  than  direct. 

Consider  a  charged  droplet  having  a  total  charge  Q0  and  electrical  conductivity  a.  Upon 
contact  with  the  stopping  copper  plate,  charges  will  begin  transferring  to  the  electrometer.  The 
amount  of  free  charge  remains  after  time,  t,  is  given  by  the  expression, 

Q(t)  =  Q0e-°t/e  (4.10) 

where  e  is  the  permittivity  of  the  media.  Since  most  dielectric  fluids  have  a  dielectric  constant  of 
about  two,  we  will  approximate  e  by  2%  where  Eq  is  the  permittivity  of  free  space. 

The  conductivity  consists  of  two  components;  one  due  to  the  residual  conductivity  aQ, 
the  other  due  to  the  presence  of  excess  charge  carriers  given  by  PoPe  which  is  the  product  of  the 
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mobility  and  chaise  density  of  the  electrons  in  the  media.  The  conductivity  is  then  expressed  as. 


o  =  c0  +  PoPe  (4-n) 

For  hydrocarbons,  the  values  of  a0  are  usually  very  small,  typically  of  the  order  10'13 
mho/m.  The  mobility  of  electron  in  Octoil  is  approximately3  10'8  mho/C.  In  this  experiment, 
typical  charge  densities  for  the  droplets  are  between  10  and  105  C/m3.  Thus,  we  may 
approximate  equation  (4.1 1)  by, 


o  -  P0Pe 


(4.12) 


The  fractional  charge  lost  after  time,  t,  can  be  written  as, 

f  =  (l-e-‘/t)  (4.13) 

where  x  =  e  /  o,  is  the  time  constant  for  charge  relaxation  .  If  ^  is  the  time  allowed  for  charge 
transfer,  the  fractional  charge  lost  in  time  t0  is  given  by, 

ft  =  (1  -  e  *  l0/T  )  (4.14) 

the  actual  size  of  the  droplet,  ra  can  then  be  calculated  from  the  measured  radius,  r^  using  the 
relationship49, 

ra  =  (ft)-1/3rm  (415) 

The  correction  factor  ( ft )  * 1  /  3  is  a  function  of  the  volumetric  charge  density  (and  therefore  a 
function  of  the  charge-to-mass  ratio)  of  the  droplets.  Figure  4.11  shows  the  correction  to  the 


Figure  4.1 1  Effect  of  Finite  Conductivity  on  Octoil 


Droplet  Radius  (Microns) 


measured  radius  due  to  charge  relaxation.  The  correction  is  negligible  at  charge  densities  higher 
than  103  C/m3. 


These  corrected  data  showed  excellent  correlation  to  Hendrick's  data  (Figure  4.12). 
Since  the  Hendrick's  data  are  not  affected  by  the  finite  electrical  conductivity,  they  are  used  as  a 
standard  to  check  the  feasibility  of  the  conductivity  argument.  Figure  4.12  shows  that  the 
conductivity  assumption  offers  a  pausible  explanation  for  the  deviation  observed.  It  would  also 
be  worthwhile  in  future  tests  to  enhance  the  conductivity  through  additives.  This  should 
eliminate  the  need  for  the  correction  term  discussed  in  this  section. 
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Chapter  V 
Summary 


Theoretical  modeling  indicates  that  droplet  charging  in  electrostatic  sprays  will  exhibit 
phase-transition  like  behavior.  A  test  facility  has  been  constructed  to  verify  this.  A  low 
charge-to-mass  ratio  quadrupole  mass  spectrometer/high  speed  electrometer  apparatus  is  used  to 
determine  the  droplet  size  and  charge  distribution  of  charged  droplet  sprays.  These  data  are  the 
first  of  their  kind  to  be  obtained  from  electrostatic  sprays.  Comparison  with  data  published  by 
Hendrick's  has  shown  good  corrrelation. 

Necessary  software  has  been  developed  to  enable  evaluation  of  various  spray  models 
with  the  data.  The  maximum  entropy  model  is  compared  with  the  Octoil  data  as  an  initial  test  of 
the  facility.  Several  predictions  of  the  model  are  confirmed  through  this  analysis.  The 
phase-transition  like  behavior  which  manifests  itself  as  a  multimodal  size  distribution  is 
confirmed  at  low  charge-to-mass  ratios.  Droplets  at  very  large  radius  are  shown  to  be 
describable  by  a  universal  law  (the  asymptotic  B  branch).  This  branch  has  the  special  property 
of  being  independent  of  the  fluid  properties. 

All  Lagrange  multipliers  introduced  in  the  maximum  entropy  model  are  now  numerically 
defined  through  an  extension  of  the  model  equations. 

The  effective  surface  electron  emission  field  strength  is  found  to  be  larger  on  the 
asymptotic  B  branch  compared  to  the  smaller  droplets  on  the  A  branch.  This  result  can  also  be 
independently  obtained  from  comparison  of  the  theoretical  and  measured  charging  characteristics 
of  the  spray.  This  probably  indicates  a  difference  in  surface  condition  for  the  two  spray 
branches. 

Charge  relaxtion  due  to  finite  electrical  conductivity  of  Octoil  at  low  charge  densities  was 
found  to  cause  underestimation  of  the  droplet  size.  A  correction  factor  is  derived  and  the 
corrected  data  show  very  good  agreement  with  the  previous  data  and  the  model. 
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Appendix  A:  Minimum  Energy  Model 


The  minimum  energy  model  proposed  by  Vonnegut  and  Neubauer16  considers  a  volume  of 
liquid  V  and  total  charge  Q  undergoing  breakup.  The  final  state  of  the  breakup  process  is  a 
group  of  smaller  droplets  of  equal  size  each  containing  equal  amount  of  charge  q.  The  final 
droplet  radius  and  the  total  number  of  droplets  in  the  end  state  are  then  determined  from  the 
condition  of  minimum  energy,  i.e.,  the  final  state  must  have  equal  or  less  energy  than  the  initial 
state. 

The  minimum  energy  condition  then  requires  the  setting  of  the  derivative  of  the  combined 
surface  energy  and  electrostatic  potential  energy  with  respect  to  radius  to  zero.  From  this 
condition,  an  expression  for  the  final  radius  of  the  droplet  can  be  obtained. 

The  condition  of  equal  final  droplet  size  and  charge  makes  the  derivation  very  restrictive. 
In  fact,  the  theory  applies  only  to  monodisperse  sprays.  In  Hogen's  paper,  the  derivation  is 
carried  further  to  give  the  charge-to-mass  ratio  as  a  function  of  radius, 
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M  3/2 

Pr 


(A.l) 


where  eQ  is  the  permittivity  of  free  space,  y  is  the  surface  tension  of  the  fluid  and  p  is  the 
density  of  the  fluid. 

Equation  (A.l)  shows  that  the  charge-to-mass  ratio  scales  inversely  as  the  three-half  power 
of  the  radius.  This  describes  Octoil  and ,  glycerine  spray  data  fairly  well.  However,  the 
correlation  with  data  for  Wood's  metal  is  not  very  good. 

Observations  of  electrostatic  sprays1,  have  indicated  that  the  droplets  produced  are  far  from 
uniform  in  size.  Furthermore,  the  charge-to-mass  ratios  of  the  droplets  are  not  always  equal  as 
assumed  by  Vonnegut.  The  conditions  for  the  monodisperse  spraying  are  highly  sensitive  to  flow 
rate,  charging  voltage,  and  fluid  properties.  Therefore,  the  validity  of  the  model  will  depend  on 
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these  controlling  factors.  Failure  of  the  model  to  describe  Hendrick's  data  could  be  a  violation  of 
these  basic  assumptions  of  the  model. 
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Physical  insight  is  the  motivation  for  this  model.  It's  been  observed  that  a  charged  liquid 
meniscus  forms  Taylor  cones^.13  when  acted  upon  by  an  electric  field.  At  the  apex  of  these 
cones,  liquid  is  drawn  out  into  filaments  due  to  electric  field  concentration.  These  electrified 
liquid  filaments  are  unstable  to  perturbations;  thus  they  are  prone  to  develop  instabilities  where 
they  eventually  breakup  into  droplets. 

Based  on  the  work  of  Melcher15  on  field-coupled  surface  waves,  Schneider8  developed  a 
model,  which  was  later  improved  by  Neukermans14.  The  model  relates  the  charge-to-mass  ratio 
of  a  droplet  to  the  charging  voltage  on  the  fluid  meniscus  and  the  geometry  of  the  charging 
electrode.  The  expression  for  charge-to-mass  ratio  is  given  by. 
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where  a  is  the  radius  of  the  liquid  jet  before  instability  sets  in,  b  is  the  radius  of  the  cylindrical 
charging  electrode,  q  is  the  density  of  the  fluid  and  Va  is  the  charging  voltage  with  respect  to 
ground. 

The  radius  of  the  droplet  after  breakup  can  be  expressed  as  r(j=(3a^Vj/4f)l/3,  where  f  is  the 
frequency  of  breakup  and  vj  is  the  velocity  of  the  liquid  jet.  In  this  case,  the  charge-to-mass  ratio 
is  not  strictly  a  pure  function  of  the  radius  of  the  droplet  but  is  a  function  of  the  charging  voltage, 
and  jet  velocity  as  well. 

In  this  formulation,  the  Rayleigh  limit  is  not  built  into  the  model.  One  could  have  droplets 
exceeding  the  Rayleigh  limit,  thus,  rendering  the  model  unrealistic.  This  model,  similar  to  the 
minimum  energy  model,  does  not  model  the  size  distribution  of  the  spray. 

Studies  on  the  breakup  process  can  only  yield  expressions  for  the  average  charge-to-mass 
ratio  on  the  basis  of  some  amplified  disturbances  in  the  fluid  jet.  The  problem  of  size  distribution 
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of  the  electrostatic  spray  has  not  been  addressed  by  other  authors.  To  model  the  size  distribution, 
a  different  approach  must  be  taken.  The  approach  that  has  been  used  is  based  on  the  maximum 
entropy  formalism.  Kelly2  had  used  this  formalism  to  derive  a  size  distribution  for  the 
electrostatic  spray.  A  similar  derivation  by  Swellen  and  Brzustowski21  exists  also  for  droplet  size 
distribution  in  uncharged  sprays.  These  models  are  referred  to  as  first  principle  models  because 
they  utilize  the  basic  conservation  equations  for  mass,  momentum,  energy,  and  charge. 
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Appendix  C:  Operation  Principles  of_the  Quadrupole  Mass  Si 


eter 


Designed  in  1958  by  Paul33  for  fusion  fuel  separation,  the  mass  spectrometer  selects  charged 
droplets  of  given  charge-to-mass  ratio  from  the  incoming  droplet  beam.  These  droplets  follow 
stable  trajectories  through  the  spectrometer  and  are  then  picked  up  by  a  charge  detection  assembly 
situated  at  the  quadrupole  exit.  Droplets  with  the  incorrect  charge-to-mass  ratio  will  not  be 
focused  and  their  trajectories  will  diverge,  eventually  colliding  with  the  quadrupole  rods. 

The  quadrupole  does  not  operate  with  a  magnetic  field,  therefore  it  is  lighter  in  weight. 
The  attractive  features  of  the  device  consist  of  variable  resolution  at  all  charge-to-mass  ratios,  low 
power  consumption  and  high  sensitivity.  Such  devices  can  easily  scan  charge-to-mass  ratios 
over  several  orders  of  magnitude.  The  instrument  developed  for  this  electrostatic  spray 
experiment  has  a  four  order  of  magnitude  range  from  0.01  to  10  C/kg. 

A  schematic  of  the  mass  spectrometer  is  shown  in  Figure  3.1.  It  consists  of  four 
stainless  steel  rods  each  1.43  meter  long  mounted  in  a  square  array.  The  original  device 
conceived  by  Paul  requires  the  use  of  hyperbolic  electrode  surfaces.  However,  the 
manufacturing  of  such  surfaces  are  far  more  difficult  than  circular  ones.  The  best  approximate 
field  achivable  is  obtained  when  circular  rods  of  radius  16%  smaller  than  the  field  radius34  ^S-36 
is  used.  Studies  shown  that  devices  with  these  specifications  can  provide  pole  fields  that  are 
accurate  to  within  a  percent  of  the  field  generated  by  hyperbolic  rods37-43. 

To  operate  the  quadrupole,  potentials  <t>i  and  <J>2  are  applied  to  the  rods  in  diagonal  pairs 

where. 


<j>!  =  U  +  V  Coscot  (C.l) 

<J>2  “  -(U  +  V  CosCOt)  (C .  2) 

The  <j>'s  represent  two  d.c.  biased  oscillating  voltage  signals.  Each  of  the  a.c. 
components  are  180  degrees  out  of  phase  from  each  other  and  they  are  cLc.  biased  in  the  opposite 


37 


direction  by  equal  magnitudes. 

Given  these  potentials,  a  time  varying  potential  is  set  up  in  the  region  surrounded  by  the 
rods.  The  potential  in  this  field  has  the  form, 

2  2 

0=  (U  +  VCosCOt )  - — — — —  (C  .3) 

2 

r 

Charged  droplets  entering  this  time  varying  field  will  be  set  into  oscillatory  motions,  but  only 
those  with  bounded  oscillations  will  exit  the  filter  at  the  detector. 

The  governing  equations  for  the  quadrupole  do  not  take  drag  due  to  air  friction  into 
account.  Thus,  it  is  necessary  to  operate  the  quadrupole  under  vacuum.  Gravitational  force  is 
found  to  be  small  compared  to  the  electrostatic  forces  experienced  by  the  charged  droplets  and 
can  be  neglected. 

To  derive  the  equations  of  motion  for  a  charged  droplet  in  a  quadrupole  field,  we  take  the 
gradient  of  the  potential  to  obtain  the  time  varying  electric  field.  Assuming  a  droplet  of  mass  m 
and  charge  Q,  the  equations  of  motion  can  be  expressed  as: 


mX  +  2Q(U  +  VCosCOt)  —  -  0  (C .  4 ) 

tx  2 

r 

mY  -  2Q(U  +  VCosCOt)—  -  0  (C.5) 

tt  2 

r 

mZtt-  0  (C . 6) 

Equation  (C.6)  indicates  that  motion  in  the  axial  direction  is  not  affected  by  the  applied  field. 
Equations  (C.4)  and  (C.5)  are  not  in  standard  canonical  form.  Transforming  the  equations  with 
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the  introduction  of  a  dimensionless  time  variable  5  such  that  tot=2^,  we  obtain  the  equations, 


X^+  (a  +  2q  Cos2^)  =  0  <C .  6) 

Y^-  (a  +  2qCos2£)  =  0  (C.7) 

where  a  and  q  can  be  expressed  as, 
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In  these  forms,  they  are  the  well  known  Mathieu's  Equations,  the  solution  of  which,  known  as 
the  Mathieu's  Function,  describes  the  trajectories  of  the  charged  droplets  in  the  quadrupole  fieid. 
Mathieu's  Equations  are  linear  differential  equations  with  periodic  coefficients41.  They  exhibit 
bounded  stable  solutions  only  for  limited  values  of  a  and  q  41 34.35,36  (Note:  Stable  solution  in 
this  case  refers  to  bounded  solution,  there  is  no  requirement  that  the  solutions  be  periodic). 

Figure  C.  1  shows  the  stability  region  in  the  a-q  space.  The  clear  regions  denote  the 
values  of  a  and  q  in  which  stable  solutions  exist,  whereas  the  cross-hatched  regions  correspond 
to  the  values  of  a  and  q  for  unstable  solutions.  In  the  quadrupole  mass  spectrometer,  there  are 
two  pairs  of  rods  independently  connected  to  different  signal  sources.  This  gives  two  equations 
of  motion  (equation  C.4  and  C.5)  for  a  single  charged  droplet  entering  the  quadrupole.  The 
motion  along  the  axis  of  the  quadrupole  is  not  influenced  by  the  applied  electric  field.  Equation 
C.6  indicates  that  the  charged  droplet  will  retain  its  velocity  when  it  is  injected  into  the 
quadrupole. 

Equation  C.6  and  C.7  have  coefficients  that  are  different  in  signs  causing  their  stability 
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regions  in  the  a-q  space  to  be  different.  For  a  charged  droplet  to  successfully  pass  through  the 
quadrupole  both  solutions  to  the  Mathieu's  Equations  must  be  bounded.  The  stability  diagram 
for  the  quadrupole  therefore  consists  of  the  superposition  of  the  two  individual  stability  curves. 
The  intersection  of  the  stable  regions  are  then  the  stable  regions  governing  the  operation  of  the 
quadrupole  mass  spectrometer. 

To  operate  the  quadrupole,  a  point  in  one  of  these  stable  regions  of  the  a-q  diagram  is 
chosen  as  the  operating  point.  From  either  expression  (C.8)  or  (C.9),  the  parameter  a  or  q  can 
be  related  to  the  charge-to-mass  ratio  of  the  droplet  and  other  quantities  such  as  the  frequency  of 
the  a.c.  signal  and  amplitude  of  the  applied  potential.  The  region  closest  to  the  axis  is  chosen  for 
this  device.  It  is  also  a  common  practice  among  designers  of  quadrupole  devices.  Other  "higher" 
stable  regions  on  the  stability  map  could  be  used,  but  require  higher  operating  voltages  on  the 
quadrupole  rods  and  are  usually  not  practical36. 

Figure  C.2  shows  an  enlarged  view  of  the  first  stable  region  for  the  quadrupole.  Having 
chosen  the  operating  point  (a,q),  U  can  be  calculated  if  V  is  prescribed  or  vice  versa.  Dividing 
equation  C.8  by  C.9  we  obtain,  a  /  q  *  2U/V.  It  should  be  noted  that  once  these  voltages  are 
chosen,  they  do  not  depend  on  other  parameters  of  the  quadrupole.  From  the  chosen  value  of  a 
or  q,  the  scanning  charge-to-mass  ratio  can  be  calculated  from  equation  C.8  or  C.9  as  a,  q  and 
the  field  radius  r  are  fixed.  The  charge-to-mass  ratio  is  therefore  a  function  of  the  applied 
frequency  of  the  A.C.  voltage. 

It  is  possible  to  operate  the  quadrupole  in  two  different  modes.  One  of  them,  described 
in  the  previous  paragraph  requires  fixing  the  quadrupole  rod  voltage  and  allowing  the  variation  of 
the  charge-to-mass  ratio  by  changing  the  A.C.  frequency.  The  alternate  operating  mode 
corresponds  to  fixing  the  applied  quadrupole  A.C.  frequency  while  die  scanning  charge-to-mass 
ratio  is  changed  by  varying  the  amplitude  of  the  applied  potentials.  Both  modes  are  equivalent 
but  the  actual  implementation  depends  on  the  capability  of  the  power  source  available.  In  the 
present  design,  the  first  mode  is  used. 

Resolution  of  the  instrument  is  defined  as  the  ratio  of  the  scanning  charge-to-mass  ratio  to 
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its  bandwidth  at  that  ratio.  In  principle,  one  would  like  to  have  infinite  resolution.  In  Figure 
C.2,  the  apex  of  the  stability  region  corresponds  to  the  point  of  infinite  resolution.  Operations  of 
the  quadrupole  at  this  point  would  permit  no  droplet  to  pass  through.  Drawing  a  line  from  the 
origin  to  the  chosen  operating  point,  we  obtain  the  mass  line.  Once  the  operating  point  (a,q)  is 
defined,  droplets  of  different  charge-to-mass  ratio  would  correspond  to  different  locations  on  the 
mass  line.  On  this  line  the  charge-to-mass  ratio  increases  as  one  moves  away  from  the  origin. 

Notice  that  only  a  small  portion  of  this  line  is  within  the  stable  region.  Droplets  with 
charge-to-mass  ratios  corresponding  to  the  stable  portion  of  the  mass  line  would  be  selected 
(Figure  C.2).  Therefore,  the  narrower  the  stable  portion  of  the  mass  line,  the  higher  the 
resolution  of  the  quadrupole.  In  practice,  one  would  like  to  operate  very  close  to  the  apex  to 
achieve  a  high  resolution.  On  the  other  hand,  one  should  also  consider  the  fact  that  the  higher  the 
resolution,  the  more  restrictive  the  droplet  selection  process.  It  will  reduce  the  quadrupole 
transmission  efficiency  (i.e.  the  amount  of  droplets  to  be  collected  will  be  reduced).  To  have  a 
reasonable  capture  rate  of  the  droplets  the  resolution  should  not  be  too  high.  The  operating  point 
chosen  for  this  device  in  the  a-q  plane  is  (0.223,0.69).  This  corresponds  to  a  scanning  resolution 
of  20,  which  means  that  if  the  scanning  charge-to-mass  ratio  is  10  C/kg,  droplets  with 
charge-to-mass  ratio  between  9.5  to  10.5  C/kg  will  be  selected. 

To  ensure  that  a  specified  resolution  is  attained,  the  motion  of  the  droplets  along  the 
quadrupole  axis  must  be  considered.  Although  the  applied  fields  do  not  affect  the  droplets  axial 
motion,  their  axial  velocities  are  retained  and  govern  the  resident  time  of  the  droplets  in  the  field. 
The  resident  time  of  the  droplets  determines  the  effectiveness  of  the  quadrupole  as  a  filtering 
device  because  the  quadrupole  removes  unwanted  droplets  by  causing  their  trajectories  to  diverge 
and  these  diverging  trajectories  require  finite  time  to  develop. 

Therefore,  the  higher  the  resolution,  the  more  selective  the  filtering  process  and  requires 
the  droplets  to  perform  more  oscillations  for  a  finer  selecting  process.  Unfortunately,  there  is  no 
straight  forward  rule  for  determining  the  number  of  oscillations  required.  A  general  rule  for  the 
required  number  of  particle  oscillations  can  be  obtained  from  related  ion  beam  research34.  It  is 
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entirely  empirical  and  should  be  use  with  caution.  The  number  of  oscillations,  N,  a  droplet 
needs  to  perform  during  its  passage  through  the  quadrupole  is  expressed  in  terms  of  the 
resolution  R  where, 


N  -  (3.5  R)l£  (C.10) 

For  a  given  resolution,  the  number  of  oscillations  required  for  complete  removal  of 
unwanted  droplets  can  be  determined  from  equation  (C.10).  Within  the  transit  time,  all  droplets 
must  perform  at  least  N  oscillations.  Insufficient  resident  time  would  result  in  a  decrease  in 
effective  resolution.  Since  the  extractor  electrode  is  maintained  at  ground  potential,  droplets 
formed  will  be  accelerated  by  the  atomization  potential.  Assuming  that  the  droplet  carries  a 
charge  Q  and  mass  M  then  energy  conservation  gives, 

Q  V  =  0.5  M  v2  (C.ll) 

where  V  is  the  charging  potential  with  respect  to  ground  potential  and  v  is  the  final  velocity  of  the 
droplet. 

The  above  expression  permits  the  final  velocity  to  be  expressed  as  a  function  of  the 
charge-to-mass  ratio  and  the  atomization  potential.  The  transit  time  in  the  quadrupole  is  then 
calculated  from  the  length  of  the  quadrupole  field.  The  length  of  the  quadrupole  rod  is  taken  to 
be  the  length  of  the  field  i.e.  end  effects  have  been  neglected.  Multiplying  the  transit  time  by  the 
scanning  frequency  enables  us  to  determine  the  number  of  oscillations  the  droplets  perform. 

In  the  preliminary  design  of  the  quadrupole,  high  rod  voltages  were  not  possible.  As  a 
result,  the  number  of  oscillations  were  insufficient  for  the  designed  resolution  of  20.  Attempts 
have  been  made  to  decelerate  the  droplets  using  a  deceleration  electrode  placed  at  the  quadrupole 
entrance.  This  method  did  not  work  because  the  velocity  of  die  droplet  scales  as  the  square-root 
of  the  potential  drop.  A  voltage  comparable  to  the  charging  potential  is  required  to  slow  the 
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droplet  sufficiently.  The  method  failed  because  application  of  high  potential  at  the  deceleration 
electrode  hinders  the  smooth  operation  of  the  atomizer  due  to  their  close  proximity.  This  causes 
frequent  arcings  across  the  decelerating  electrode  which  inhibit  the  smooth  operation  of  the 
sprayer. 

An  alternate  method  was  then  used  to  increase  the  number  of  oscillations  of  the  droplets 
during  transit  through  the  quadrupole.  This  is  achieved  by  increasing  the  scanning  frequency  of 
the  quadrupole  which  increases  the  number  of  oscillations  for  the  same  atomization  potential. 
However,  to  maintain  the  same  charge-to-mass  ratio,  higher  d.c.  and  a.c.  amplitude  must  be 
used  in  the  quadrupole  signals.  By  using  transformers  to  amplify  the  existing  input  voltages  the 
operating  point  can  be  maintained  at  higher  frequencies  (cf.  C.8  and  C.9).  Since  these 
transformers  were  originally  designed  to  operate  at  a  fixed  frequency  of  400  Hz,  their  effective 
gain  changes  when  other  frequencies  are  applied.  The  gain  on  these  transformers  tend  to 
decrease  with  increasing  operating  frequency.  As  a  result,  the  present  quadrupole  is  only  capable 
of  scanning  up  to  2.9  C/kg. 


.Quadrupole  Fawer  SubpIy  Circuit 

The  voltage  source  is  a  Behlman  variable  A.C.  power  supply  capable  of  a  maximum  RMS  value 
of  270  volts.  Calculations  show  that  for  our  device  an  A.C.  amplitude  of  approximately  3000 
volts  is  required  to  maintain  a  scanning  resolution  of  20  assuming  no  deceleration  electrode  is 
used.  To  obtain  the  high  voltages  two  1000:1  voltage  transformers  are  used  to  amplified  the 
quadrupole  signals.  Figure  C.3  is  the  circuit  diagram  of  the  voltage  amplifier  circuit.  The 
original  waveform  from  the  Behlman  unit  is  duplicated;  they  are  made  180°  out  of  phase  from  the 
one  another  and  are  amplified  by  the  two  amplifiers. 

The  required  biased  voltage  is  supplied  by  a  Kepco  power  supply.  In  the  experiment,  it 
is  required  to  operate  the  quadrupole  at  frequencies  between  50  to  2000  Hz.  Since  the  gain  of  the 
transformers  cannot  be  maintained  constant  at  different  frequency  and  that  the  transformers  have 
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Figure  C.3  Circuit  Diagram  for  Voltc&e*  Amplifier 


